texts: Obtaining regular raster grid from netCDF format grid data in ArcGIS Desktop or R?

texts:
I downloaded E-OBS gridded data from European climate data assessment (grid datasets website and dataset 1 and dataset 2), and I am interested in 0.25 degree grid data of land surface temperature and precipitation for Germany only, but E-OBS grid data is conducted on whole Europe and all measurement in daily level. However, this grid datasets were divided and each chunk has 15 years climate data observation, and I need 30 years data. So I downloaded two big .nc file from this site, but this file is quite big. I want to read those in ArcGIS desktop as raster and intend to pull out Germany. So I did these in ArcGIS:

from Multidimension Tool.tbx -> Make NetCDF Feature Layer, to see its attribute table of those grid data in ArcGis, but I have only coordinate pair while respective daily temperature observations' attribute table is not shown. Here is the screenshot of reading .nc file as raster in ArcGIS:

+-----+-------------+------------+------------+
| FID |  longitude  |  latitude  |     tg     |
+-----+-------------+------------+------------+
|   1 | -14.6250000 | 25.3750000 | 17.3799992 |
|   2 | -14.3750000 | 25.3750000 | 17.3799992 |
|   3 | -14.1250000 | 25.3750000 | 17.3500004 |
|   4 | -13.8750000 | 25.3750000 | 17.2899990 |
|   5 | -13.6250000 | 25.3750000 | 17.3400002 |
|   6 | -13.3750000 | 25.3750000 | 17.3299999 |
|   7 | -13.1250000 | 25.3750000 | 16.9400005 |
|   8 | -12.8750000 | 25.3750000 | 16.7600002 |
|   9 | -12.6250000 | 25.3750000 | 16.4300003 |
|  10 | -12.3750000 | 25.3750000 | 16.1099987 |
|  11 | -12.1250000 | 25.3750000 | 15.8699999 |
|  12 | -11.8750000 | 25.3750000 | 15.7500000 |
|  13 | -11.6250000 | 25.3750000 | 16.1299992 |
+-----+-------------+------------+------------+

Desired output:

I want to make raster grid from this .nc file and want to have attribute table where grid node' coordinate with respective daily temperature observations.

I am not familiar with working .nc files in R, so I want to read those netCDF file in ArcGIS and want to get raster grid first, then crop all grids of Germany where geo-coordinate must come along with respective daily temperature observation.

How can I make this happen in ArcGIS?

How can I convert grid data in .nc file to regular raster grid in ArcGIS?

Can anyone help me how to get this done easily in geospatial software tools?


texts: Calculating ratios between two shapefiles

texts:
So I have two shapefiles, one is a grid and one is a map of the world's land. I am trying to find the ratio of land in each cell of the 'grid' shapefile to show the land/water ratio but I'm stuck.
Here are the shapefiles
import geopandas as gpd
from shapely.geometry import point, Polygon, mapping, shape

grid =gpd.read_file('grid.shp')
World_Land = gpd.read_file("World_Land.shp")

for index1, row1 in World_Land.iterrows():
    for index2, row2 in grid.iterrows():
        if row1['geometry'].intersection(row2['geometry']):
            grid.loc[index2,'ratio'] = (row1['geometry'].area / row2['geometry'].area)

grid.plot(column = 'ratio')

This gives me something but I want to make sure I am making a choropleth map correctly and calculating the ratios correctly.


texts: How to make my PostGIS vector layer in QGIS visible?

texts:
I have loaded a shp file (simple polygon with one entry) into PostgreSQL 9.3 (with extention PostGIS v2.1.3).  I now want to load this layer into QGIS (v2.2.0 Valmeira) via the "load Postgis layer) command.  The connection is OK and the layer is transfered to QGIS without a problem.  I can see the table attributes in QGIS and the "Zoomin to the layer" fct works (the function zooms in with the correct center point and scale), but the layer is invisible!  Regardless of the symbology attributes that I use, I cannot see the layer.

I have also made a view of the table and loaded it into QGIS (specifying my coord syst - 4269 - and a int4 primary key) and I have the same problem.  Any idea of what the problem is?


texts: Get both the old and new zoom values in OpenLayers 4

texts:
In an OpenLayers 4 application, I want to refresh the view of a layer when the user crosses a particular zoom level, either zooming in or zooming out (the tile rendering changes at that point from less to more complex or vice versa).

I can capture change:resolution, and get the new value, but is there a way of getting the old value?


texts:
I have a polygon shapefile with holes (representing Islands and Rivers) within it as seen below:

I need an automated way to fill the polygon in QGIS as seen below: That is remove all Islands and divisions within and outside the polygon.

I know I can creating a separate layer and manually digitize the area, but there are many of such polygons.

I tried the following, non worked:

1- Vector > Geometry Tools > Multipart to Singlepart

2- Vector > Geoprocessing Tools > Dissolve

3- Vector > Geoprocessing Tools > Eliminate Silver Polygons

Graphic/image explanation above:
What I have is the first image, but what I want is the second image.

Edited:
The answer here How to smooth/generalize a polygon in QGIS over generalized the boundaries into straight line (polygon shape deformed) and it does it for separate Island. What I wanted is to eliminate the gaps (Islands) in between and keep the bigger polygon retaining the shape as close as possible.

The shapefile of the image above is here


texts: How do I get the area of a WGS84 polygon in square meters?

texts:
The following gives me the area of the polygon in square degrees:

SELECT ST_Area(ST_GeomFromText('POLYGON((871325.790874952 6105405.3261047,871418.748307692 6105359.72944624,871346.22022442 6105215.141258,871254.85408906 6105261.72007212,871325.790874952 6105405.3261047))',4326)) As sqm;

Result: 16586.2319335938

How do I get the area in square meters? I tried the following:

SELECT ST_Area(ST_Transform(ST_GeomFromText('POLYGON((871325.790874952 6105405.3261047,871418.748307692 6105359.72944624,871346.22022442 6105215.141258,871254.85408906 6105261.72007212,871325.790874952 6105405.3261047))',4326)),31467) As sqm;

But it seems like ST_Transform doesn't work. I get the following error message:

ERROR:  function st_transform(geometry) does not exist
LINE 1: SELECT ST_Area(ST_Transform(ST_GeomFromText('POLYGON((871325...
                       ^
HINT:  No function matches the given name and argument types. You might need to add explicit type casts.
********** Error **********
 ERROR: function st_transform(geometry) does not exist
SQL Status:42883
Hint:No function matches the given name and argument types. You might need to add explicit type casts.

Missing Proj support shouldn't be the problem, with SELECT PostGIS_full_version() I get

"POSTGIS="1.5.3" GEOS="3.3.3-CAPI-1.7.4" PROJ="Rel. 4.8.0, 6 March 2012" LIBXML="2.8.0" USE_STATS"


texts:
For some reason I am unable to edit my map in ArcGIS. When I click 'start editing' on the editor toolbar, I get the following error message:

I've ensured that my field names in the Excel table don't have any spaces or special characters. I've also tried adding an ObjectID field, but it won't recognize this.I'm specifically looking to draw lines on the map and write labels to identify certain urban districts.


texts:
Does anyone have a simple work-flow for drawing cross sections of earth features in QGIS?

My thing is that QGIS is being used to present the plan view, but we have a little interest in some elevations as well. All the data is one set and is simply an elevation value in the attribute table, I don't want to have to move all that to CAD.

One logical workaround is to present the z value as the y.

I am a dunce in paper space in CAD, and would be more comfortable with the QGIS Composer.


texts: Create Friendly Labels for ArcGIS Print Service Options with JS API

texts: How can I change the print options received from an ArcGIS print service to be more friendly and descriptive using the Esri sample print service?


texts: How to let user only rearrange layers in a group layer?

texts:
We have a stand-alone application were we let the users add local shape and raster layers. We would like to add functionality so that the user can rearrange the added layers within a specific group layer.

We know it's possible to enable drag and drop on the TocControl but it lets the user rearrange all layers and that's not an option for us. Is there a way to enable drag and drop only on a specific group layer or somehow mimic the same behavior?


texts: QGIS 2.8 (openSuse 13.1) pdf export with print composer not working

texts:
I'm using QGIS 2.8 with a openSuse 13.1 system. If I try to export a map (with openlayers-plugin active) with 120 dpi it works fine, but if I change the dpi setting to something higher (like 130 dpi) the map dissapears from the screen and a pdf-file with the size of 0 byte will be created. So after one hour of rendering I aborted the process.

In the terminal the following will be showed:

  [WORKER THREAD] Waiting for the async request to complete
  Speicherzugriffsfehler (segmention fault)

Are there any missing components (e.g. some python-stuff) I need to install, or other workarounds?


texts: Qfield viewing features in layer without geometry

texts:
I created Qfield project that includes layer without geometry. I can successfully create new features but couldnt find a way to edit or view created features on mobile device. I looked through documentation but could not find anything about it.

Since features in layers with geometry can be viewed by selecting on the map i presume viewing features without geometry is not (yet) possible within Qfield. Or am i missing something?


texts: Dropdown menu to select vector layer from current displayed layers in qgis

texts:
I'm trying to write a python plugin to create an input file for an urban climate model out of three shape files. One of the shape files contains the buildings (polygon), one the trees (point), and the third one the center point of the new model input file (point).

I therefore need a function to select one layer out of the currently displayed ones. I tried messing around with

qgis.utils.iface.legendInterface().layers()

but then I get some ugly output like

qgis.core.QgsVectorLayer object at 0xb4b6f5c.

Is there a possibility to get readable names out of this? Or is there an easier way to select the layers?

Any ideas or help greatly appreciated!


texts:
I have more than 1000 polygons, but they do not have the correct size and location. I have a table with the correct coordinates (centroid) "x, y" and  different scale factor for each polygon. How can I automatically resize them and move them all to their correct location? In QGIS I used use the plugin "Transformation Afine" for individualy polygons or scaling tool in advanced edit (ArcMap). My shapefile use EPSG: 5367 (local system CRTM05)

The batch processing does not recognize fields, only values


texts: Holding and using Romanised names with Diacritics

texts:
I have two True type font files (Permanent Committee on Geographical Names) PCGNA.TTF AND PCGNTN.TTF.

The romanised diacritics are well read in MS excel.

But when I convert this excel in shp or pgdb feature class, it becomes unreadable.

Does anyone know how to overcome this problem?


texts: PostGIS create topogeom from existing geometry column

texts:
I tried to create topogeom using postgis.topology extention.
SELECT topology.AddTopoGeometryColumn('roadtopo','public','roads','topogeom ','LINESTRING');

UPDATE roads SET topogeom = topology.toTopoGeom(geom,'roadtopo',1,0.00028);

But I receive an error;

ERROR:  column "topogeom" of relation "roads" does not exist
LINE 1: UPDATE roads SET topogeom = topology.toTopoGeom(geom,'roadto...

What I am doing wrong?


texts: Looking to color code census block groups by different categories in different columns in QGIS

texts:
I'm very new to QGIS and I am looking for the best way to color code census block groups by race. I have a CSV with columns for each race with the values in percents. What I'm looking to do is color code each block group as the race with the largest percent. For example, if a block group is %60 white, 20% Black and %20 Asian, that would be the color code for white.
Is something like this possible?


texts: I was trying to convert a very simple Erdas Imagine AOI file (its a single polygon) into a shapefile and was mostly finding documentation on the inverse (SHP to AOI). I was trying to follow the steps listed here, but keep getting the following "There are No Selected Elements". Below is a workflow that worked for me in case it is helpful for anyone else.


texts:
I am new in Arcmap and I am trying to calculate the distance between some points I have in an excel file that have an FID shape and latitude and longitude. However whichever tool I am trying to open eg. near, proximity I see only an empty grey box and no drop down menu like the tutorials.

What am I doing wrong?

Is there any other way to calculate the distances?


texts: Maintaining default value for fields in ArcGIS online map / feature service

texts:
I have an enterprise geodatabase within which there are several feature classes that contain subtypes and domains.  A web map has been created from this geodatabase on ArcGIS Online and it is used for data input by many other users.

Within the enterprise GDB, I set the default value for one of the fields in my feature classes and for all subtypes within the feature class.  The web map was then updated with changes.  When I add features directly into the enterprise GDB, the default value is honored within this field.  However, whenever I try to add new features to the web map at ArcGIS Online or on my Collector App, the default value does not populate the field for which it is set.

Do ArcGIS Online maps honor default values for fields?  If so, is there something that I need to be doing to get the default value to automatically fill the field?


texts: openlayers draw feautre(line, polygon) callback

texts:
I need to enable feature drawing functionality. Code:

var drawingControl = new OpenLayers.Control.DrawFeature(someLayer, OpenLayers.Handler.Path);
map.addControl(drawingControl);

But I also need to do some actions when drawing is done, on double click.
Therefore, how can I define a dblclick callback (for example, add there console.log("1"))?


texts:
I'm doing tests to query layers in Leaflet with WMS.

I'm seeing this

https://github.com/mmsmdali/ll1.0.3/blob/master/GetFeatureInfoWMS.htm and

http://jsfiddle.net/mmsmdali/ybzkg1tL/.

I'm trying to take WMS data from my Geoserver and replaced the WMS URL, the name of the layer, and the name of the fields in the popup window but I have a problem:
I see the layers in the map, but it does not make me select, and I do not see the fields in the popup.

Do I need to change some settings in my Geoserver?


texts:
I'm attempting to determine the 2D area beneath around 3,000 polyline features as in the diagram below:

(Polyline feature stretches from upflow to downflow end)

I have a DEM representing a flat surface (the base of the feature) and another containing the relief of the feature, the part I'm interested in is the two-dimensional area beneath the curve. I'm not too sure how to determine this parameter, does anybody have any ideas on this?

I'm using ArcGIS 10.4.1 and all extensions are available.


texts: how to fix this ERROR: Start vertex was not found

texts:
I am attempting to make an SQL selection but receive an error. My statement:

SELECT *
 FROM network
 JOIN
 (SELECT * FROM shortest_path('
 SELECT gid AS id,
 start_id::int4 AS source,
 end_id::int4 AS target,
 "length"::float8 AS cost
 FROM network',
 10,
 12,
 false,
 false)) AS route
 ON
 network.gid = route.edge_id;

The error:

ERROR:  Start vertex was not found.
********** Error **********

ERROR: Start vertex was not found.
SQL state: XX000


texts:
I've encountered an issue implementing a code into ArcMap. The original code was:

for fc in fclist:
... arcpy.Clip_analysis(fc, r’C:\GIS_LabManual\Lab10_Manual\Lab10_data\Arb_clip_poly.shp’, r’C:\GIS_LabManual\Lab10_Manual\Lab10_results\’+fc)

The error it shows is "Parsing error SyntaxError: invalid syntax (line 2)."

I then tried to adjust it to the following code (which also had the same error):

for fc in fclist:
... arcpy.Clip_analysis(fc,r’F:GIS_Joseph\Lab10_Joseph\Data_Joseph\Lab_10_data\Lab_10_data\Arb_clip_poly.shp’,r’F:GIS_Joseph\Lab10_Joseph\Results_Joseph\’+fc)

Could you tell what the error might be?


texts:
I'm having trouble with the Alter Field  Data Management tool -- it will not allow me to create a field name longer than 31 characters. The table I'm using is a File Geodatabase table, where the character limit on field names is 64 characters.

When importing and creating tables, I'm allowed to have field names (and aliases) that are longer than 31 characters, but when I try to alter the field name, this restriction is placed.

Can't seem to find information on this anywhere.


texts:
I have a small geoprocessing service that reads data-driven pages from an MXD, then exports a PDF.

 mxd = arcpy.mapping.MapDocument(mxd_path)
 mxd.dataDrivenPages.exportToPDF(r"\\test_report.pdf", "RANGE", page_range, "PDF_SINGLE_FILE", 300, "BEST", "RGB", True, "ADAPTIVE",
                            "RASTERIZE_BITMAP", False, True, "NONE", True, 80, False)

Where in the ArcGIS Server directory do I need to place the MXD I am trying to reference?

E.g. C:\arcgisserver\directories\arcgisjobs ?


texts: What is the FME Equivalent to the ArcGIS Union tool?

texts: I'm attempting to union 2 large featureclasses in ArcCatalog, QGIS, and FME.  What is the transformer that would be equivalent to the ArcGIS Union tool?


texts: How to resolve QGIS 3 syntax error from processing.algorithmHelp?

texts:
I am using QGIS 3.4.3-Madeira and OSX Mojave 10.14.3.
I am trying to use QGIS 3 through the Python Console.
I am following the QGIS 3 documentation.

I am trying to determine the syntax of the Split Vector Layer function.
I entered this code into the QGIS Python Console:

import processing
processing.algorithmHelp(qgis:splitvectorlayer)

I got this this syntax error:

    Traceback (most recent call last):
  File "/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/code.py", line 64, in runsource
    code = self.compile(source, filename, symbol)
  File "/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/codeop.py", line 168, in __call__
    return _maybe_compile(self.compiler, source, filename, symbol)
  File "/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/codeop.py", line 99, in _maybe_compile
    raise err1
  File "/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/codeop.py", line 87, in _maybe_compile
    code1 = compiler(source + "\n", filename, symbol)
  File "/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/codeop.py", line 133, in __call__
    codeob = compile(source, filename, symbol, self.flags, 1)
  File "<input>", line 1
    processing.algorithmHelp(qgis:splitvectorlayer)
                                 ^
SyntaxError: invalid syntax

I also tried the following but I got the same syntax error:

import processing
processing.algorithmHelp(native.splitvectorlayer)

How do I find the syntax and parameters for this function?


texts: Creating observation points(vertices) around a line feature in ArcGIS 10.1

texts: Is there any way to create observation points(vertices) or lets just say points on a line feature. What I have is a line feature and I wanna put observation points around that feature, each one with a distance of 10 Meters to the other one. Does anybody know a way to do that? Maybe a tool or something.


texts: Saving changed positions in SQL Server table X, Y columns

texts:
I have a SQL Server table with X,Y columns linked to QGIS 3.6, and I made a virtual layer of it using "MakePoint(X,Y,CRS_id) as geometry".

Now I want to shift a point in the map and have the change stored in the X,Y columns of the SQL Server table.
What I see is that I can not edit the virtual layer, so I can not shift its points. So, how to shift the points then? And if shifting points is possible, how to store the X,Y in the SQL Server table?
Our SQL Server version does not support spatial, still version 2005 (!).


texts: Publishing a geoprocessing service that imports modules besides arcpy?

texts:
I'm trying to publish a geoprocessing service (10.2.1) that will be used to generate a PDF from a webmap using an MXD layout as a template.

The problem is that if I try to import any module besides arcpy, the publishing just runs forever until I kill it.

For example, I tried to import os, and it just hangs while it's publishing. When I remove import os, it publishes just fine. Same thing happens with any other module I try to import.

Any idea what's wrong here?

I found this thread that has the exact same problem, but there wasn't really a solution:

https://geonet.esri.com/thread/86236

I tried @KHibma's suggestion, and created a "script" with only import os, and no parameters at all. Of course it runs with no errors, but the same thing happens when I try to publish the geoprocessing tool. It just hangs until I decide to kill it.

I also did the same thing, but with only import arcpy, and it publishes fine.


texts: Something going wrong trying to transform longitude, latitude to NAD 83 using pyproj

texts:
I'd like to convert a latitude longitude pair to NAD 83.  I've read NAD 83 is equivalent to epsg:4269 and the shapefile I am working with is a US Census TIGER shapefile, all of which are supposed to be in NAD 83 and geopandas reports my GeoDataFrame as being in epsg 4269. All of this makes sense.  The issue I'm running into occurs when I want to plot a point described as a lat/long pair onto the shapefile.  My understanding is that, in order to do this, I need to transform the point using the relevant map projection (NAD 83/epsg:4269), however, when I try to do this using pyproj, the point returned is way off:

import geopandas as gp
import matplotlib.pyplot as plt
from shapely.geometry import Point
import pyproj

gdf = gp.read_file('./state_shapefiles/WI/tl_2016_55_bg.shp')
projection = pyproj.Proj(init='epsg:4269')

lat = 43.087806
lon = -89.430121

# this returns tuple(-1.5608500619069237, 0.7520240821605678) which is
# nowhere near the point I am trying to indicate on my map (Madison
proj_coords = projection(lon, lat)

point = Point(proj_coords[1], proj_coords[0]).buffer(1.0)

layer = gp.GeoDataFrame({'geometry': [point]})

layer_plot = layer.plot(color='red')

gdf.plot(ax=layer_plot)
plt.show()

print('end')

Am I misinterpreting something or have implemented this wrong?

Edit: Fixed latitude and longitude mix up, issue persists.


texts: Iterating polygon creation in ArcGIS ModelBuilder?

texts:
I have contour lines in a feature class. I covert the closed lines into polygons using the 'Feature to Polygon' tool. With this method I get concentric stripes, but I need the area of each area enclosed by a polygon.

For better understanding see the picture below. So I have the two contours - A and B. If if convert them by using the 'Feature to Polygon' tool on the whole featureclass I get the areas in option one: area enclosed by contour A, and a stripe between A and B. But I need option 2: the area enclosed by A and the area enclosed by B.

To do this I tried to run the 'Feature to Polygon' tool in model builder using on the output of a Feature Selection iterator, but I got an empty feature class.

What am I missing or the iterator tool is a wrong choice?


texts:
I am working on a map mixing OpenLayers and leaflet libraries (it is possible, correct  ?).
I am such a beginner that, after lots of readings, I thought creating pop ups was less difficult for me with leaflet than OpenLayers.

After hours spent to display pop ups, I can't find the right manner. I get a

TypeError:a.div is undefined

HTML :

<title>OpenLayers Example</title>

<link rel="stylesheet" href="./leaflet-0.7.3/leaflet.css" />


texts: How can the pixel values of a raster image be extracted under a polygon fishnet?

texts: I have a single band raster image of 30m by 30m resolution. I have already created a polygon fishnet (grid) of 300m by 300m with point labels over that. How can I extract the mean pixel value of the raster image overlaying each feature in the fishnet grid? Is it possible with Arcgis?  Please inform me.


texts:
I have to lat/lon Points. I want to draw a line segment between the 2 (in kml), but I don't see functionality in Google Earth to do this. LineString is not it. Surely I am missing something.

<?xml version='1.0' encoding='UTF-8'?>
<kml xmlns:ns="http://earth.google.com/kml/2.0">
  <Folder>
    <open>1</open>
    <Placemark>
      <Point>
        <coordinates>-107.980362141, 31.872060174</coordinates>
      </Point>
    </Placemark>
    <Placemark>
      <Point>
        <coordinates>-107.963280345, 31.80230482</coordinates>
      </Point>
    </Placemark>
    <Placemark>
      <LineString>
        <extrude>1</extrude>
        <tessellate>1</tessellate>
        <coordinates>-107.980362141, 31.872060174, 0  -107.963280345, 31.80230482, 0</coordinates>
      </LineString>
    </Placemark>
    <GroundOverlay>
      <Icon>
        <href>redBorder.png</href>
      </Icon>
      <LatLonBox>
        <north>31.82103476</north>
        <south>31.78357488</south>
        <east>-107.91400752</east>
        <west>-108.01255317</west>
      </LatLonBox>
    </GroundOverlay>
  </Folder>
</kml>


texts: ST_Intersect and ST_Equals return false positive?

texts:
I try to find equal lines between two table using the following script

SELECT tableA.ID, tableB.ID tableA.geom from tableA, tableB
where ST_equals(tableA.geom, tableB.geom)

But my problem is that it seem to return all rows (~17000 rows).
I get the same result if I am using ST_Intersect.

I know that the overlapping geometries are exactly the same, they come from the same source. Not sure what is going on here? I am expecting about 1400 equal/intersecting records.

EDIT: Deleted sample data as it did not really contributed to the question.


texts: Vegetation indices extraction from RGB images

texts:
I have been working on the project to classify crops as health or unhealthy. For this task I have taken the open source dataset from the the website https://www.sensefly.com/education/datasets . It has the RGB crop field images and orthomosaic is also available.
I have been reading research paper regarding this and most of the approaches are classification using vegetation indices (NDVI).
I am facing problem in extracting this detail. The method to extract NDVI are being applied on satellite imageries.
How can I extract this data from my aerial images?


texts: Generate LiDAR point cloud understory division through histogram and height threshold

texts:
This is somehow a good post that gives me a bit of direction;
Retrieving LiDAR metrics from relative vertical layers of point cloud

I'm using FUSION and rLIDAR most of the time. I'm looking for a way to detect understory crown for Individual tree. So the steps done so far, I've had done generated CHM, applied watershed delineation and extracted the raw LiDAR returns from each segmented area. Meaning i had saved all the point clouds into separate las files represent individual tree. So there's thousands of trees for my study area. Right now, i'm looking for a way to plot a histogram and get the height threshold and then separate the point clouds based on the threshold (because later i need to process the lower layer CHM for the understory).

From @Andre Silva answer in that post this is something that i can interpolate, but can anyone advice me in terms of giving me example step using R if possible on how to do this histogram profiling, thresholding thing and separate the point clouds based on the threshold.


texts: How to get Point feature from Mouse interaction in openlayers 3

texts:
I am unable to get point features from event coordinate like below

 $(map.getViewport()).on('click', function (evt) {
            //var pixel = map.getEventPixel(evt.originalEvent);
            var coordinate = map.getEventCoordinate(evt.originalEvent);
layer.getSource().getFeaturesAtCoordinate(coordinate);
}

It's returning empty array. But if I use map.forEachFeatureAtPixel(pixel,function) its working but it is not dependable.
Is there any solution?


texts: QGIS raster to vector with shapefile feed procedure

texts:
Does a procedure like Polygonize raster to vector exist in QGIS that can do the following:

I have a raster image, a shapefile (.shp) with all the regions of interest and I need a vector output with those polygons only.


texts: Installing QGIS from source in Suse Enterprise Desktop 11

texts:
I have installed QGIS from source in Ubuntu using the this link. In ubuntu I compiled qgis from source without any problem.

but now I moved to Suse Enterprise Desktop 11. In internet every one pointing to this link. But For Suse 11 there is no repositories available.

So In Suse,

how to compile Qgis from sorce?
Where can I get all qgis dependency list for suse 11?


texts:
In my project, I use 3 layers with a line geometry.
Every line in the layer is a contour with the same depth "valdco".

For operational ease, I color these lines blue, red en orange (as seen in the screenshot).

To color these lines I use rule-based symbology.
When I need to change the value of the blue line I open the properties of the layer and change the value for the blue line in the rule, i.e. "valdco"= 'value'.

I would rather like to have a toolbar where I can set the value for the colored line.

My idea is to change back to a simple line and use the data override function on the color of the line:

CASE
 WHEN "valdco" = 'value_set_in_dialogbox1_of_toolbar' THEN BLUE
 WHEN "valdco" = 'value_set_in_dialogbox2_of_toolbar' THEN ORANGE
 WHEN "valdco" = 'value_set_in_dialogbox3_of_toolbar' THEN RED
 ELSE BLACK
END

Question: How do I make a toolbar with 3 dialogboxes where I can set the desired value for the colored lines. Is there a plugin available or a decent tutorial written?

I've started to build the mentioned dialogboxes in Qt5 and made the layout as I would like it to have in QGIS. But this is my absolute first time.

What is my next step?

Can I make a simple database (with no geom) and then have the values from my plugin write-out to a field in the database?
Can I directly reference, the values in the plugin, to the symbology of the chart, so to manipulate the colors?
I am drawn to the signal/slot editor but not sure this is the right method.


texts: Hiding particular field in parameter combo box using ArcObjects and C#?

texts:
I am currently fetching the all fields/attributes of feature class and display in Combo box parameter using below code:

paramEdit = new GPParameterClass();
paramEdit.DataType = new DEFeatureClassType() as IGPDataType;
paramEdit.Value = new DEFeatureClass() as IGPValue;
paramEdit.Direction = esriGPParameterDirection.esriGPParameterDirectionInput;
paramEdit.DisplayName = "Feature Class 1";
paramEdit.DisplayOrder = inputTable1Index;
paramEdit.Enabled = true;
paramEdit.Name = "in_line_featureclass1";
paramEdit.ParameterType = esriGPParameterType.esriGPParameterTypeRequired;
paramEdit.Domain = lineFeatureClassDomain as IGPDomain;
paramArray.Add(paramEdit);

paramEdit = new GPParameterClass();
paramEdit.DataType = (IGPDataType)new FieldType();
paramEdit.Value = (IGPValue)new Field();
paramEdit.Direction = esriGPParameterDirection.esriGPParameterDirectionInput;
paramEdit.DisplayName = "Select attributes from Feature Class 1";
paramEdit.DisplayOrder = inputDynamicList1Index;
paramEdit.Enabled = true;
paramEdit.Name = "in_line_featureclass1_attrib";
paramEdit.ParameterType = esriGPParameterType.esriGPParameterTypeOptional;
paramEdit.AddDependency("in_line_featureclass1");
IGPFeatureSchema outputSchema = new GPFeatureSchemaClass();
IGPSchema schema = (IGPSchema)outputSchema;
schema.CloneDependency = true;
paramEdit.Schema = outputSchema as IGPSchema;
paramArray.Add(paramEdit);

I want to hide some fields which display in this combo box parameter.

Is it possible to hide the fields or display selected fields from feature class ?

As shown in above image when I selecting feature class it display the fields of that feature class in combo box but I want to hide some fields in that combo box so how can I do this?


texts:
I have point data along lines shown in picture by black color. I made a grid by inverse distance method in QGIS. Maybe the result is not that appealing to eye or something (especially when eye is used to Oasis montaj minimum curvature grids).

Do you have any hints how to make my grid in QGIS looking more like that in Oasis.

I am not an expert on grid algorithms so maybe inverse distance is not even the best choice here but QGIS does not offer minimum curvature.

EDIT NOTE: Added new example dataset to demonstrate difference between QGIS and Oasis montaj grids.

Picture 1. QGIS grid (inverse distance)
Picture 2. Oasis montaj grid (minimum curvature)
Picture 3. Original data

Data is supposedly gravimetric data from a mine site.


texts: Limiting number of duplicate labels using QGIS

texts:
I am using QGIS 3.10.

I am trying to draw a map of the US interstate highway.
How could I limit the number of labels to just a handful per each interestate highway?
I have tried all the various settings but haven't got the required result.


texts: Why OSM data has different attribute values from the same element and same time frame?

texts:
Here is my question: Does anyone know why OSM has different attribute values from the same element and same frame time?

The example explains itself, when I try the Shapefile from different mirrors like Planet.osm.org , bbbike.org and geofabrik.de; I find that they have a mistake with the Direction of Travel on its attribute value.

When I check in the website: here the link , I see on its attribute value that there is no mistake at all and has its right value.

So, I thought that those mirrors extract the exact information from osm servers, but I think that there is a lack of info. Is that right?

Thanks for your time,


texts: How can I prevent ogr2ogr from trying to get the SRID from PostGIS?

texts:
I'm running ogr2ogr for a PostGIS to GeoJSON export like this:

ogr2ogr -s_srs EPSG:4326 -f GeoJSON unique_geoms.json "PG:host=localhost dbname=local_data user=local_user password=local_pass" -sql "SELECT MAX(geom)::geometry FROM geom_table GROUP BY ST_AsBinary(geom)"

ogr2ogr isn't very bright, and tries this query on my database:

SELECT ST_SRID("max") FROM(SELECT MAX(geom)::geometry FROM geom_table GROUP BY ST_AsBinary(geom)) AS ogrpggetsrid LIMIT 1

So the database has to run the entire subquery to even get a response! This is fine for small tables, but for really big tables, the operation takes forever.

As you can see, I tried using -s_srs EPSG:4326 to override the built-in SRID detection, but it doesn't work.

My query is designed to de-duplicate my point geometry. I could create a new table based on this query, but I'm not sure that would help.


texts: Google Maps: Decoded polylines showing up incorrectly

texts:
I'm trying to draw an encoded polyline via the Google Maps API (v3), and it seems to draw 'sort of' correctly (as in, roughly in the right area of the world ;) ) but otherwise is very, very wrong. (Pictures to follow)

Essentially, I'm using the Strava API (v3) to get road 'segments', which include encoded polylines (http://strava.github.io/api/#polylines). The short version is, I'm decoding these as per the Google documentation and then creating a new Polyline object with the new path and adding it to the map.

I've been very brief in summing up that part, simply because my method doesn't seem to be the issue here... If I take the encoded polyline string directly from the response of the API call and paste it into Google's Interactive Polyline Encoder Utility (https://developers.google.com/maps/documentation/utilities/polylineutility) I see the same 'wrong' path.

So, that leaves two possibilities... Either the Strava API is passing out incorrect data, or I'm not correctly understanding what I'm supposed to do with it.

I've attached the relevant images + polyline below. Am I doing something obviously wrong?

Example 'encoded' polyline:

{ae{HntiQtCcDzG_I|^uc@rFgHhC{CxAiA~AaA~BkAvB}A|F_G|AgBbBkCtAwCd@sA|BoIVw@Pc@|@gBt@}@|@y@lCwBvA_B`@k@~@aBt@iBlAaE~@oEp@sDX{BP_BJaDAcEIeCe@gHo@yMUaEk@uDm@iD]mCAwBNsDXyDL}@nByIZyCt@cLr@gNB_ABoEAkFJmDTkBVeAZ_Af@gAnDwF|@gBbAoChHgUPWlAT`@B|@GbE_@dAW`Cu@vBe@tDs@xD{@`Bg@bBq@hBaAtB}@dCi@bF}@jBg@pBeAj@SNE\\C^@\\DbAZ`Ah@~C`A\\H|ALzAFLA^Gl@UdBgAjBaBZSh@Qz@MjD_@`FoAtCa@j@Ez@DxE|@xF\\nBP~@TxHvBf@Tb@\\pBvC\\^`@XxAf@fBT|BDfAIr@MfBe@rBa@rBMvBYxBg@xA_@^Ir@@NF|@l@nBfAjAj@dBV`Bb@lBbAbB~ALPhC`FV`@n@z@^VNBX?LGZa@d@eAp@qAt@Sx@Cz@G\\IZOhCcBb@c@T]jA_CrE_HfEiFz@}@p@k@|@o@`C{A`A{@rBwBx@oAbByCp@wArAoDLWxA}BhAcBjAqAlAiB~AaDr@sBp@{CD[TkC^}FZyD^oCx@gF`@qAh@kAz@yAtAgBpD_E|JoKdDuEjBcCfC{ExCqGdAgBlBuBrAyBpEkIpEsI\\]^YbAg@|GaBzKeEfBe@lCW`AQr@U|A_AtAkAhDyCpAeA|Aq@`EeCrDgBvA{@tD}C`BmAzBm@t@QvAQxBOl@Q~Ai@~BsAlCcB

Strava map (correct):

Decoded polyline (wrong):


texts: Can I create a virtual mosaic using other virtual mosaics?

texts:
I am following this example - http://linfiniti.com/2009/09/image-mosaicking-with-gdal/

I have 6TB of aerial imagery in 6 folders 1TB each and each folder has many subdirectories including the geotiffs.
I am using gdal_translate to compress the data first. Then I use gdalbuildvrt to create the mosaics, and finally I build the pyramids/overviews for the mosaics using gdaladdo.

My target is to use geoserver to serve all the aerial data connected together. My concern is that a subdirectory may be 100GB which takes quite a while to process (either compress or building pyramids). This is why I did not use gdal on all 6TB at once and preferred handling them in different directories. Thus, my question is whether or not I will be able to use something like gdal_merge to connect all the pyramitised mosaics together to a bigger mosaic and serve it using geoserver or similar OS software.

Apologies for the huge post! Thanks in advance for any hints or tips.


texts: Getting QGIS python modules to work from PyScripter in Windows 10?

texts:
Been trying lots of "solutions" from net research but nothing seems to work.

I have the current versions of QGIS from OSGEO and PyScripter 2.5

Based on one of the solutions I tried...

import sys
sys.path.extend([r"C:\OSGeo4W64\apps\qgis\python",r"C:\OSGeo4W64\apps\qgis\bin"])
sys.path.reverse()
import qgis.core

gives me

  Message   File Name   Line    Position     Traceback
           4
        C:\OSGeo4W64\apps\qgis\python\qgis__init__.py  26       ImportError: No module named sip

import processing

gives me

  ImportError: No module named processing

and path is

PATH=C:\ProgramData\Oracle\Java\javapath;C:\Python27\;C:\Python27\Scripts;C:\Program Files (x86)\Intel\iCLS Client\;C:\Program Files\Intel\iCLS Client\;C:\Windows\system32;C:\Windows;C:\Windows\System32\Wbem;C:\Windows\System32\WindowsPowerShell\v1.0\;C:\Program Files\Intel\Intel(R) Management Engine Components\DAL;C:\Program Files\Intel\Intel(R) Management Engine Components\IPT;C:\Program Files (x86)\Intel\Intel(R) Management Engine Components\DAL;C:\Program Files (x86)\Intel\Intel(R) Management Engine Components\IPT;C:\Program Files\Intel\WiFi\bin\;C:\Program Files\Common Files\Intel\WirelessCommon\;C:\Program Files (x86)\NVIDIA Corporation\PhysX\Common;C:\Program Files (x86)\Windows Live\Shared;C:\Program Files (x86)\Skype\Phone\;C:\WINDOWS\system32;C:\WINDOWS;C:\WINDOWS\System32\Wbem;C:\WINDOWS\System32\WindowsPowerShell\v1.0\;C:\Program Files (x86)\QuickTime\QTSystem\;C:\Program Files\Intel\WiFi\bin\;C:\Program Files\Common Files\Intel\WirelessCommon\; C:\OSGeo4W64\bin

pythonpath includes

C:\OSGeo4W64\apps\qgis\bin
C:\OSGeo4W64\apps\qgis\python

Tried "solutions"

https://mapoholic.wordpress.com/2012/06/28/configure-pyscripter-qgis/
How can I Setup Python to recognize QGIS
(and links from it)
https://stackoverflow.com/questions/15297391/problems-while-importing-pyqgis-modules-in-python
Import error for qgis.core when running OSGeo4w shell script

and

Problem with import qgis.core when writing a stand-alone PyQGIS script

after running the .bat in above I get "ImportError: no module name qgis.core" in python within osgeo shell as well.


texts: How to add Audio into QGIS recorded and geotagged using OSMand app

texts:
Is there a plugin or easy way of adding audio recordings, recorded and tagged using OSMand app, into QGIS.
I'm using a PC and have no coding knowledge!
I was able to  load the photographs using the "load photos" plugin then pop them up using the eVis plugin.
I understand that I could also use eVis for audio but seems to require a database set up which I don't have.


texts:
I have a small .OSM map which has some rules set up based of whether it's deemed a primary road, secondary road, path etc. What I would like to do is branch each one of these out to its own layer for exporting to a DXF. I've tried to duplicate the layer and removed the conditional formatting but when I export the DXF it still contains all of the data.

Is there any way to do this?


texts: Leaflet Draw Popup text unselectable in Chrome

texts:
I have a map here: http://maps.gcc.tas.gov.au/measure.html

I have some simple popups happening, and it's working well. But, I can't select the text for markers when using Chrome or Firefox.

This does work in IE, funnily enough. What is going on?


texts: Import ARCGIS web dataset contours into Sketchup

texts:
My local council has online ARCGIS data sources including 1m elevation contours.  I am trying to import these into a sandbox terrain model in sketchup to do some landscape planning.

Conceptually I know what I need to do, and can generate the necessary REST API queries to get the data, but am just not able to string it all together.

The contour dataset is here:

http://gisservices.scc.qld.gov.au/arcgis/rest/services/Elevation

An example bounding box for the query is {17008011.278814,-3090615.13623876,17008111.278814,-3090515.13623876}

I can dump the contours in KMZ, and they look great in google earth, however Sketchup (15.3) is unable to import the KML, it just gives a generic 'Import error' dialog.

Also, looking at the KMZ data it doesn't actually have the elevation data, it's pure 2D contours which Google Earth is just clamping to ground level.

If I do a JSON export of the contours, and include the 'CONTOUR' field in the output dataset then I get lots of contour polylines (and elevations) but then I am stuck trying to get that into Sketchup.

I should be clear, I don't have the ARCGIS software at all, I'm trying to do this purely by scraping the REST API of the online data repository.  This is just a private/hobby project.


texts: Derive (default) value in ValueTable based on values in ValueTable of Python Toolbox using ArcPy?

texts:
I am writing a Python Toolbox which needs Shapefiles (or Featureclasses) as input and calculates an output for each of the input files. I would like to have the possibility to rename the files in order to avoid an overwriting or an exception.
My approach is to use a toolbox parameter of type GPValueTable with 2 columns. The input files is given in the first column (by the user) and I want the 2. column to be automatically filled by the base name without extension of the input file.

Unfortunatelly, I am not able to fill the second column based on the first.
I tried various ways already but here is a snipped of sample code:

    def getParameterInfo(self):
    """Define parameter definitions"""

    input_paths = arcpy.Parameter(
        displayName="Input files",
        name="in_paths",
        datatype="GPValueTable",
        parameterType="Required",
        direction="Input",
        multiValue=True)
    input_paths.columns = [['GPFeatureLayer','Eingabedatensatz'], ['GPString','Ausgabename']]

    output_path = arcpy.Parameter(
        displayName="Output path",
        name="out_path",
        datatype="DEFolder",
        parameterType="Required",
        direction="Input")

    params = [input_paths,output_path]
    return params

def isLicensed(self):
    """Set whether tool is licensed to execute."""
    return True

def updateParameters(self, parameters):
    """Modify the values and properties of parameters before internal
    validation is performed.  This method is called whenever a parameter
    has been changed."""
    if parameters[0].value:
        new_vt = []
        gp_value_table = parameters[0].value
        vt_length = len(gp_value_table)
        if  gp_value_table[vt_length-1][0] and not not gp_value_table[vt_length-1][1]:
                    new_vt = [[v[0],v[2]] for v in gp_value_table]
                    input_path = gp_value_table[vt_length-1][0]
                    new_vt.append([vt_length-1,arcpy.Describe(input_path).baseName])
                    parameters[0].value = new_vt
    return

In this case nothing happens and the 2. column remains emtpy.
How can I - row by row - derive the default output name based on the input?


texts: How do I load offline google map/earth vrt files as a backdrop in qgis

texts:
I was trying to follow the process : saving-google-earth-location

in order to cache google maps/earth offline for digitizing, but I am fairly new to qgis. When I try to add the vrt file as a raster in qgis I get nothing but a blank image. I also noticed that pyramiding gave me an error which I suspect is because of the resolution. In addition if I wanted to load the vrt as a backdrop/base map does it not require a bunch of tiles loaded to create an image mosaic? I am working in Windows so I would prefer to not have to run some gdal commands in the command prompt.

I'm also not so sure how to use tile caching and I would prefer the easiest way possible without having to standup any wms type of server.

*I know its against the terms of service.


texts: I don't think this exists but just wanted to double check with the GIS community. Is there a tool or plugin that connects Google Earth and ArcMap such that when you are in Street View in Google Earth there is dynamic point in ArcMap that moves around relative to your position in Google Earth?


texts: drop list from field to execute SQL request

texts:
I'm trying to create a python script that allows a selection attribute based on the name of a field in the attribute table of a layer.

I would like to configure my script such a way that permits the user to select the SOURCE NAME from a drop-down list.
Here's the code I have written:

import arcpy
arcpy.env.workspace = 'in_memory'
arcpy.env.overwriteOutput = True
mxd = arcpy.mapping.MapDocument("CURRENT")
df = arcpy.mapping.ListDataFrames(mxd)[0]
PosteLayer=arcpy.GetParameterAsText(0)
field="NAME_SOURCE"
rows=arcpy.da.SearchCursor(PosteLayer,["NAME_SOURCE"])
list=[]
for row in rows:
   if str(row[O]) not in list:
       liste.append(row[0])
POSTES=liste
POSTES=arcpy.GetParameterAsText(1)

I'm using the Advanced Version Of Arcgis 10.2


texts: PyQt resources file to display images on QGIS custom feature form

texts:
I am having issues creating a PyQt4 resources file to display images in my QGIS 2 custom feature form using Windows XP. EDIT This still gives me problems with windows 7 on 32 and 64 bit machines, so I'm pretty sure this isn't system related. Any suggestions to check if pyqt or pyrcc is installed correctly..?

These are the steps I followed:

-Added QLabel and clicked "choose resource…"

-Created my .qrc in PyQt4 designer, linked .png image file. Checked my file with notepad, and following lead of some tutorials made a few modifications (because it didn't work with the untouched resource file direct from PyQt4) My .qrc looks like this:

<!DOCTYPE RCC><RCC version ="1.0">
    <qresource prefix="/">
        <file>img/image.png</file>
    </qresource>
</RCC>

-Saved the ui and compiled the resource file in the windows command line using:

pyrcc4 -o image_rc.py image.qrc

-I even tried compiling the ui to py and moved the import space_rc to the top. I'm pretty sure this has nothing to do with the qrc file though, as it is for compiling the ui file.

pyuic4 -o form.py form.ui

then:

import image_rc

the image is called like so:

self.label_7.setPixmap(QtGui.QPixmap(_fromUtf8(":/img/image.png")))

my file system structure is like so:

ui/resource.qrc
ui/form.ui
ui/img/image.png

Under my layer's properties in QGIS, I chose to add custom ui and linked to the ui created in PyQt designer. When opening the feature form, I see only my frame, but no image. After reading many tutorials on what steps to follow, I can't see what I am missing.

I had believed it was a environment variable PATH issue (conflicting with GDAL), but now I switched computers and my issue is back. The windows environmental path for PyQt is the very first one, and I still don't see my image.

Path     C:\Python27\Lib\site-packages\PyQt4;...


texts: How to extract data from a ESRI ArcGIS MapServer?

texts:
I have connected my QGIS to the Sao Paulo Hydropower potential MapServer with a WMTS connection using the following URL:

http://tps.geosur.info/arcgis/rest/services/maps/GeoSUR_SPHydropower/MapServer/WMTS/1.0.0/WMTSCapabilities.xml

Visually this is great but how can I access the data source behind this specific layer: Hydropower Potential MW (ID: 7)?. In other words, how can I click on a coloured line in QGIS to see what the potential value is?.

I noticed there is a link to a json version but not sure how to load this into QGIS (or if it is indeed the best method).


texts: I have a QGIS (2.14.3) viewshed analysis that is taking too long to run on my desktop computer.  Is it possible to run QGIS on the Amazon cloud computing (or some other similar service) to speed up the process?


texts: Import vector data from PostGIS to an R session

texts:
So now I know that it has been asked before many times but I wanted to know what are the main problems that one encounters when using rgdal. To be more specific I have 2 lines of code:

dcV="PG:dname='remote' host=localhost user='postgres' password='pass' port=5432 scehma='public' table='vector2' mode=2"
poly=readOGR(dcV)

but I have I get an error that says:

  Error in ogrListLayers(dsn = dsn) : Cannot open data source

and I have no idea what is causing it because I loaded a raster layer the same way but instead of "vector2" I wrote "raster2".And yes both of them exist and can be fetched from QGIS
P.S I am working with Linux and I am calling the Rscript from the terminal

EDIT:
When I typed drv = ogrDrivers();drv[grepl("PostGIS",drv$long_name),] I got the following output:


texts: How to create a web app that will display data in unique urls for the end-user?

texts:
Can anyone provide insight on how a novice can construct a fairly rudimentary web app? I would like to send out unique URLs that contain my finished maps instead of exporting them into PDFs and sending them by email. I am using ESRI software and making my final maps in ArcMap 10.1.

The maps each contain 1 polygon, 4 buffers, and between 1 to 500 vector points. I would like the end user to be able to see the points I've mapped for them, zoom in and out of the map area, and click on points to see more data about them. They should not be able to edit any data or add any additional points to the map. I have no prior experience in creating mapping applications.


texts: Draw Coordinates problem in QGIS map composer

texts: I'm new to QGIS and am having trouble with the Draw Coordinates feature in Map Composer. Extensive Googling has not provided any solutions. With reference to the attached image, the coordinates around the page perimeter vary between UTM and degrees, which I find odd. I'm suspecting something to do with the project CRS.


texts: Performing multiple "Closest Facilities" analyses for multiple feature classes

texts:
I am trying to calculate network distance from 10 points within 896 polygons "SA1s" (a total of 8960 facilities), each to ten points within 186 polygons "DZNs" (a total of 1860 incidents). I would like to iterate over each of the 896 facility groups (each has it's own feature class, containing 10 points (I would undertake some post-processing of these feature classes, deleting the main outputs with each one.

The iteration is not particularly problematic, but I'm having a problem with the actual python script running over a single iteration.

I set up a script in Model Builder on a test set of data (10 facilities from one SA1 to 10 incidents from one DZN) and ran it - it ran perfectly. I then deleted the outputs and exported the model to a python script, shown below:

# Import arcpy module
import arcpy

# Check out any necessary licenses
arcpy.CheckOutExtension("Network")

# Local variables:
Incidents = "in_memory\\{DC034A5A-77FD-4B9E-BD36-ECDE5ABD576B}"
Facilities = "in_memory\\{0AA3C7A5-47AA-494D-B97E-6C5A25C7F522}"
in_memory = "in_memory"
Point_Barriers = "in_memory\\{DDEAEEC2-AAD9-4B99-9D69-1020AEEF8B6F}"
Line_Barriers = "in_memory\\{C7F6D7A5-B9FB-447B-8373-E98398408F2C}"
Polygon_Barriers = "in_memory\\{40ABDBEB-B48B-4C95-8354-9C3DFB60EDB9}"
Attribute_Parameter_Values = "in_memory\\{D84163A6-7A47-419B-B8FF-79420ABC0651}"
DZN_Test_Points =    "C:\\Users\\xgoldie\\Desktop\\ModeShare_Study\\Mode_Share_Study.gdb\\DZN_Test_Points"
SA1_Test_Points = "C:\\Users\\xgoldie\\Desktop\\ModeShare_Study\\Mode_Share_Study.gdb\\SA1_Test_Points"
ACT_Road_Network_ND_nd = "C:\\Users\\xgoldie\\Desktop\\ModeShare_Study\\ACT_Road_Network_ND.nd"
Mode_Share_Study_gdb = "C:\\Users\\xgoldie\\Desktop\\ModeShare_Study\\Mode_Share_Study.gdb"

# Process: Find Closest Facilities
arcpy.FindClosestFacilities_na(DZN_Test_Points, SA1_Test_Points, "Meters", ACT_Road_Network_ND_nd, Mode_Share_Study_gdb, "Routes_Test", "Directions_Test", "ClosestFacilities_Test", "10", "", "TRAVEL_TO", "", "NOT_USED", "GEO_LOCAL", "ALLOW_UTURNS", Point_Barriers, Line_Barriers, Polygon_Barriers, "Minutes", "Minutes", "Meters", "Meters", "NO_HIERARCHY", "", Attribute_Parameter_Values, "", "20 Kilometers", "\"ACT_Road_Network\" #;\"ACT_Road_Network_ND_Junctions\" #", "TRUE_LINES_WITH_MEASURES", "10 Meters", "NO_DIRECTIONS", "", "Miles", "", "", "", "", "", "", "", "", "NO_SAVE_OUTPUT_LAYER")

I then attempted to run it again from IDLE, and it failed. Giving me the following error message:

Traceback (most recent call last):
  File "C:\Users\xgoldie\Desktop\SCRIPTS\networkanalyst.py", line 30, in <module>
    arcpy.FindClosestFacilities_na(DZN_Test_Points, SA1_Test_Points, "Meters", ACT_Road_Network_ND_nd, Mode_Share_Study_gdb, "Routes_Test", "Directions_Test", "ClosestFacilities_Test", "10", "", "TRAVEL_TO", "", "NOT_USED", "GEO_LOCAL", "ALLOW_UTURNS", Point_Barriers, Line_Barriers, Polygon_Barriers, "Minutes", "Minutes", "Meters", "Meters", "NO_HIERARCHY", "", Attribute_Parameter_Values, "", "20 Kilometers", "\"ACT_Road_Network\" #;\"ACT_Road_Network_ND_Junctions\" #", "TRUE_LINES_WITH_MEASURES", "10 Meters", "NO_DIRECTIONS", "", "Miles", "", "", "", "", "", "", "", "", "NO_SAVE_OUTPUT_LAYER")
  File "C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy\arcpy\na.py", line 2595, in FindClosestFacilities
    raise e
ExecuteError: A geoprocessing error occurred.
Failed to execute. Parameters are not valid.
ERROR 000735: Input Rows: Value is required
Failed to execute (GetCount).
Failed to execute (FindClosestFacilities).

Does anyone why this is occuring, and how to fix it?


texts: How to configure Proxy.cgi file for Geoserver On Tomcat?

texts:
I have a web map application in my local machine(windows xp and dont have IIS) and I have postgis data in Geoserver on Tomcat (Ubuntu). To get the WMSgetfeatureInfo, I have configured the proxy.cgi file in Geoserver  (usr/lip/cgi-bin/proxy.cgi) and added my local ip address in

allowedHosts = ['xxx.xxx.x.xxx', 'www.openlayers.org', 'openlayers.org', 'labs.metacarta.com', 'world.freemap.in', 'prototype.openmnnd.org', 'geo.openplans.org', 'sigma.openplans.org', 'demo.opengeo.org', 'www.openstreetmap.org', 'sample.azavea.com', 'v2.suite.opengeo.org', 'v-swe.uni-muenster.de:8080', 'vmap0.tiles.osgeo.org','www.openrouteservice.org']

and added the below line in my script.

OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url=";

but it throws an error in chrome console:

GET localhost:3890/usr/lib/cgi-bin/….…T%3Dapplication%252Fvnd.ogc.gml%26SRS%3DEPSG%253A27700%26X%3D682%26Y%3D284 404 (Not Found)

I dont know whether its right or wrong. Proxy.cgi should be copied in Geoserver or my machine? in which path? added my IP into allowedhosts is correct?

Please anyone suggest me to fix this issue


texts:
I am trying to preprocess a netcdf dataset with arcgis such that it has the same extent, cellsize and format as an ascii dataset I am using. Both are visualized in the image below:

The big one in grey is the netcdf dataset, it has a longitude range of 0-360. The small dataset in Black is ascii, it has a range of -180-+180 decimal degrees. Problem is, the netcdf dataset is split right at the greenwich meridian, which splits africa, my region of interest. Now, i can change the prime meridian in the dataset to -180°, which makes the image look right, and both layers are congruent. Problem is, I want to clip the large dataset with the extent of the small dataset (rectangular), but since the minimal X is -30, the clipped image is cut exactly at 0 degrees. I tried to shift the whole dataset to the left (X-180°) but then again, it cuts off at 0°

Thanks for the help


texts: Has anyone use AcuGIS, a web hosting service for PostGIS and GeoServe?  I am curious to see how it works with QGIS.  I am looking at some cloud hosting services for PostGIS and I would like to be able to access the database in QGIS and pgAdmin.  Has anyone had any success or issues with AcuGIS.  Or can you recommend any other cloud hosting services.


texts: Using NASA's GIBS servers with GRASS GIS

texts:
I am trying to download WMS tiles into a grass mapset. I'd like to use NASA's GIBS servers which replaced the JPL OnEarth services. I'm not sure how much remains the same between the two services, but it seems that the JPL OnEarth server was supported. Here is an example command that I run in the Grass shell

r.in.wms url=https://gibs.earthdata.nasa.gov/wmts/epsg4326/best/wmts.cgi layer=BlueMarble_ShadedRelief output=testing format=png

This does connect to and download from the server and it is reported that the raster is imported. However, the entire image is populated with no data. Am I missing something really basic here?


texts: Unable to get QGIS 3.4 to work correctly with a clean install - Windows Registry Remnants

texts:
I installed QGIS 3.4 and I keep getting errors attempting to use certain tools. GRASS r.slope.aspect via QGIS 2.18 vs QGIS 3.4 - not working
I uninstalled but it is not removing all traces of the software to allow a completely clean re-install. I am left with remnants in the Windows 10 Registry for various plugins. Is this what causes my problems? These entries should not be there.

It seems to be that on Windows 10 systems:
A. 3.4.2-Madeira (released 23/11/18) has issues interacting with the GRASS tools.
B. 3.4.3-Madeira (released 21/12/18) has issues interacting with the GRASS, SAGA tools.


texts: How to "select by location" certain points

texts: I have a small problem - how to select the points on the line but not select the points at the intersections of lines using "select by location"? Illustration below. I am using ArcGIS 10.0


texts: Georeferencing map with known Central Lat-long and distances?

texts:
I have the scanned copy of an Old Map which was created using a Plane Table survey, and based on my reading of literature around it, is supposed to be quite Accurate.

The central Lat-long is known, as is the fact that this map represents an Area Of 1000 feet(North-South) by 1400 feet (East-West)

How can I Georeference this Map using just these facts?

I have access to QGIS & GDAL.


texts:
I am not yet used to QGIS (I often get lost in menus for simple things), and I would like to make shapefiles using historical maps. Here is an example: (the real map is in higher resolution).

I would like to put it in QGIS, align it to the actual geographical map, and use it to trace some of the elements and export it into a shapefile.
How can I do it?


texts: Why does arcpy.ListFields list GUID field as SmallInteger?

texts:
I want to get a list of all GUID fields in a feature class, but it seems that arcpy.ListFields is listing my GUID field as a SmallInteger. What am I doing wrong?

global_id_fieldnames = [f.name for f in arcpy.ListFields(layer) if f.type.lower() == 'guid']

after this list was empty, I called arcpy.AddMessage() for each of the fields:

for f in arcpy.ListFields(layer):
    arcpy.AddMessage('{},{}'.format(f.name, f.type))


texts: Making GEE code for Pansharpening Sentinel-2

texts:
The following link contain information about Pansharpening Sentinel-2
https://leclab.wixsite.com/spatial/post/pansharpening-sentinel-2-imagery-in-google-earth-engine
The code in picture form and long . How it can be converted into GEE JavaScript.
Please suggest some alternate code link .
Regards


texts: Error fetching feature types when adding WFS to ArcMap using Interoperability Extension

texts:
I'm trying to add a WFS layer to ArcMap using the Interoperability Extension. I get the following error when clicking the '...' next to the 'Feature Types' option.

I'm unable to add the data because of this.

Any ideas on what could cause this error? As far as I know the Data Interoperability extension and license are correctly installed and activated. The WFS link should be correct. I tried other URLs of different WFS layers and they all give the same error.


texts: Set attribute values while creating an ArcGIS Server feature in esri-leaflet?

texts:
The Esri-Leaflet sample map Editing feature layers shows how the geometry of an ArcGIS Server feature layer can be created in Leaflet, by digitising a new polygon which is inserted into the feature layer.

Is it possible to set the attributes of the newly created feature?

This code inserts the new feature into the feature layer: pedestrianDistricts.addFeature(e.layer.toGeoJSON())

Examining the value of e.layer.toGeoJSON() in the debugger I can see that there are no properties (ie attributes) on the editable layer. Is it possible to pass the desired attributes while creating the geometry, during the addFeature?


texts: Missing v.stream.order of GRASS in QGIS processing tool

texts:
I want to use v.stream.order of GRASS in QGIS for stream but it is not available in QGIS processing tool.
How to add that in processing tool

QGIS version
2.14.17-Essen


texts: "No rows were georeferenced" error in CartoDB

texts:
I am trying to geolocate my data by the city name and country, but receive an error:

  No rows were georeferenced It seems that some of your rows didn't
  finish successfully. Perhaps these rows contained empty values or
  perhaps we just didn't know what the values meant. We encourage you to
  take a look and try again.

Basically, I have an .xls dataset with columns: Data - that contain data of patent codes; City - the city where this patent registered; and Country - Italy. I've been trying to map the location for patents, but I did not succeed with even one address. What might I be doing wrong?


texts: Finding source code of "fill sinks (wang and liu)"

texts:
I would like to use the source code of "fill sinks (wang and liu)".

I checked through QGIS3 folders, but could not find much related to it.

Anyone know where the code could be?


texts: Checking if point lies with GeoJSON Polygon in PostgreSQL?

texts:
I have already parsed my GeoJSON in PostgreSQL and my table looks exactly like the answer to Storing GeoJSON FeatureCollection to PostgreSQL with PostGIS?

Now I want to query PostgreSQL with lat/lon and check if a point lies within a polygon or not and if the point lies within a polygon then print the city name.
So far I have tried this :

SELECT * FROM locations WHERE ST_DWithin(geom, 'POINT(71.013 11)', 10.0);

But this is obviously wrong because I don't want to specify a distance within which I want to check. I just want if a point lies in a certain polygon in database then print the row


texts: How to handle a UnicodeEncodeError with sampling tool?

texts:
I'm using qgis and I keep getting an error when using the point sampling tool, its keeping me from getting the sampling done. Does anyone has an idea what is it about?

Thanks!


texts:
I am using two sets of data, 'Part' and 'Total'. The part set is what you would think, a fraction of the total. Part = 13 , Total = 29.
I want to overlap this data on a map using Graduated Symbols but my problem is the 'Total' quantitiy is much greater than the 'Part' quantity. This makes my graduated symbols different scales and gives a false visual representation.
Even though both tables have different quantities, how can I make it so, visually, the symbols are still accurate to true quantity?


texts: Dynamic rendering on a layer with a long SQL query

texts:
I have a question on dynamic rendering with the CartoDB APIs. These APIs are new for me.

Data is a big layer stored in CartoDB (about 30 000 polygons).

In my application the default rendering of this layer must be replaced by a rendering defined by the user.  For each subset of selected polygons per the user, the color change. But in a subset the number of polygons can be about 10 000 polygons for example.

In my code I use the subLayerOptions by setting SQL and cartocss properties. The SQL is like "sql:"SELECT * FROM polygons where ID in ('10000','40202'...".

It works fine with a few values but it doesn't work when the SQL request contains a lot value because the request is too long.

Is there an another method to do this?


texts: Calculating Average NDVI value for research Plot?

texts:
I am working on a Field trial with small research plots. I've created a GeoTiff orthomosaic with 4 bands from a multispectral camera (Red, Green, NIR, Red Edge) and want to be able to extract the mean pixel value for each band in each plot using ArcGIS. I have tried creating a fishnet grid to cover each plot, but am having trouble getting it to align. Furthermore, once I've created my fishnet, I'm not entirely sure what to do next.

Would I use the Zonal Statistics toolbox to do this?

Also, is there a simpler way of doing this, for instance creating a simple polygon around a plot and then computing the values for each band?


texts: How to move an Overlay with animation in Openlayers 3?

texts:
I would like to animate a Overlay marker from one coordinate to another with a smooth movement.

In the last Overlay documentation:
http://openlayers.org/en/v3.9.0/apidoc/ol.Overlay.html
there are autoPan & autoPanAnimation options, but it is not working for me as you can see here:
http://jsfiddle.net/q2x4vf6x/2/

If I define the marker without those options the marker changes its position as expected, but without animation:

marker_pos = new ol.Overlay({
        element: $('<img id="marker_pos" src="http://1.bp.blogspot.com/_GZzKwf6g1o8/S6xwK6CSghI/AAAAAAAAA98/_iA3r4Ehclk/s1600/marker-green.png">')
        .css({marginTop: '-10px', marginLeft: '-10px', cursor: 'pointer'})
});

Am I missing something?


texts: Cell value based on positions of neighbouring cells

texts:
I want to know if there is a way to assign values to a raster cell based on the values and positions of neighbouring cells. To be more precise, I want to look for:

the occurance of neighbours with matching values,
the spatial relationship of matching and non-matching neighbours.

For example, these situations should result in different values (o = matching values):
1:
x x x
o o x
x x x

2:
x o x
x o x
x x x

3:
x o o
x o x
o x x

Only the studied cell (the cell in the middle) should be assigned the new value, but all cells should be studied in the same way. To preserve the original raster, I guess the ideal would be to store the new values in a different raster.
To my understanding, focal statistics doesn't deal with spatial relationships in quite this manner, but maybe there is another way?
I use ArcMap 10.4.1.


texts: How to clip polygon by line feature QGIS?

texts: The clip tool will only allow for clipping by another polygon, I would like to clip polygon features by a line feature. Had tried installing the digitizing feature however the cut by lines button is no longer there.


texts: SQL Geometry type and Geography type conversions/projections

texts:
I am working with a replica database managed in Sql Server Management Studio 2016. My geographic data of points, lines, and polygons are projected in WGS_1984_UTM_Zone_17N Meter and SRID of 32617. Its geographic system is GCS_WGS_1984 Degree and SRID 4326. Everyone of my feature classes contains a latitude and longitude field; however, all of my feature classes are geometry types. I'm in the process of developing triggers to auto-populate several fields but the lat/longs are giving me trouble.

Does anyone know of a way to project/transform geometry data into a geography variable where a lat/long can be extracted? Are there other ways to do get this information that are not accomplished within the ArcGIS applications?

The trigger I'm developing on my A-delta table shows that I am stumped after calling the x and y since I don't know how to get the lat long.

USE [DATABASE]
GO
/****** Object:  Trigger [dbo].[AreaCalc_Insert]    Script Date: 2/13/2018     12:52:59 PM ******/
SET ANSI_NULLS ON
GO
SET QUOTED_IDENTIFIER ON
GO
ALTER TRIGGER [dbo].[AreaCalc_Insert]
ON [dbo].[a101]
FOR INSERT
AS

BEGIN
DECLARE @areaCalc decimal(10, 2);
SET @areaCalc = 10.7639;
DECLARE @Lftfrommeter decimal(38, 8);
SET @Lftfrommeter = 0.3048;
DECLARE @AreaUOM nvarchar(20);
SET @AreaUOM = 'squareFoot';
DECLARE @perimUOM nvarchar(25);
SET @perimUOM = 'USsurveyFoot';
DECLARE @strt nvarchar(20);
SET @strt = '00010000000';
DECLARE @fstr nvarchar(20);
DECLARE @oid int, @IDPK int, @total int;
DECLARE @X float, @Y float;

SET NOCOUNT ON;

UPDATE dbo.a101
    SET @X = shape.STCentroid().STX,
        @Y = shape.STCentroid().STY,
        @oid = OBJECTID,
        @IDPK = DATALENGTH(@oid),
        @total = 20 - @IDPK,
        @fstr = LEFT(@strt, @total) + CAST(@oid as nvarchar(20)),
        buildingIDPK = @fstr,
        onInstallation = 'yes',
        installationID = '%%%%',
        installationName = '&&&&',
        siteID = '&&&&0001',
        realPropertySiteUniqueID = '4949',
        majorCommand = '!!!!',
        wacInnrCode = '76',
        country = 'US',
        owner = '****',
        areaSizeUOM = @AreaUOM,
        perimeterSizeUOM = @perimUOM,
        perimeterSize = Round((Shape.STLength() * @Lftfrommeter), 0),
        areaSize = Round((Shape.STArea() * @areaCalc), 0)
    WHERE OBJECTID IN (SELECT OBJECTID FROM INSERTED)

update dbo.BUILDING_A
    SET @oid = OBJECTID,
        @IDPK = DATALENGTH(@oid),
        @total = 20 - @IDPK,
        @fstr = LEFT(@strt, @total) + CAST(@oid as nvarchar(20)),
        buildingIDPK = @fstr,
        onInstallation = 'yes',
        installationID = '%%%%',
        installationName = '&&&&',
        siteID = '&&&&0001',
        realPropertySiteUniqueID = '4949',
        majorCommand = '!!!!',
        wacInnrCode = '76',
        country = 'US',
        owner = '****',
        areaSizeUOM = @AreaUOM,
        perimeterSizeUOM = @perimUOM,
        perimeterSize = Round((Shape.STLength() * @Lftfrommeter), 0),
        areaSize = Round((Shape.STArea() * @areaCalc), 0)
    WHERE OBJECTID IN (SELECT OBJECTID FROM INSERTED)

END


texts: How does my phone get an accurate reference time for calculating distance from GPS satellites?

texts:
So if I understand this correctly, my phone uses trilateration to figure out a position in 3D space using the distance from various satellites.

In order to figure out the distance from those satellites, you need to figure out the time difference between when the signal was sent and when it was received.

But how does my phone have an accurate enough time for that? Or battery powered GPS? Do those get the correct time from cell phone towers or something!


texts: How to style by date monthly ranges QGIS?

texts:
So I have mypointtable with columns id, mydate and geom and I need to style it according to the months in the date column

For example:

id | mydate
1  | 2019-11-01
2  | 2019-11-15
3  | 2019-11-28
4  | 2019-12-03
5  | 2020-01-01

I've tried using month(age(mydate, now()))  <=1

But besides being a bit slow, it just takes all the features from now and 30 days from now. So wrongly capturing ids 2,3 and 4

Intended result is to capture all features in November, so ids 1, 2 and 3 but not 4

I tried a solution for this month only which is month(mydate) - month(now()) = 0 and year(mydate) - year(now()) = 0. But the problem with this approach is that while it does work for this month, it does not work for December.

How to create a solution to style current month and the next 5 months?

By next 5 months I mean

style for current month (november at the moment)
style for next month (december at the moment)
style for in 2 months (january at the moment)
style for in 3 months (feb at the moment)
style for in 4 months (mar at the moment)
style for in 5 months (apr at the moment)
style for anything after 5 months months (may onwards)


texts:
I am trying to create a basemap in tilemill.  Very new to this software.  When I add the label styling for the roads, it adds the name several times.  Is there a way to have it only place the label on the map once?

My code

#MajorRoads{
  [STREET_NAM='NC-191']{
   text-name:[STREET_NAM];
   text-face-name:'Arial Regular';
   line-width:3;
   line-color:#606060;
  }


texts: Acceptable Reduce() amounts for google maps zoom levels?

texts: I'm using SQL Server Spatial to return large polygons to Google Maps for display.  Given that google maps has a zoom scale of 1 - 20, I would like to reduce my polygons (reduce function uses the Douglas Peucker formula) appropriately for the zoom level being displayed.  Does anyone have a table of visually acceptable levels of reduction per zoom level, or have a formula that would work well?


texts: I am using ArcGis 10 Template windows application(Desktop) in VisualStudio 2010. i have added an ArcGis Engine Map Control and engine TOC Control in that. when I am Adding my digitized map to the map control it is only showing my map not the digitized features and layers in map control. whereas in the TOC Control it is showing the name of all the layers. please help me.


texts: I added the coordinates in the map but values in the middle of the map overlap. How to fix that?


texts: Having difficulty in arcmap 10.5 changing projections from GCS_North_American_1983 WKID: 4269 Authority: EPSG to NAD_1983_Ontario_MNR_Lambert WKID: 3161 Authority: EPSG. Every attempt constantly displays the newly projected data outside the view extent of the current Arcmap layer even with the new projection. To view the new projected data I have to use "zoom to layer"  then zoom to layer back to the data in my layer.


texts: Converting integer value to unicode for labeling in QGIS

texts:
It is possible to stylize a layer with integer numbers stored in a shapefile and then choose a symbol out of a font.

I would like to use the label engine instead, because i'm more flexible, for example i can rotate the label.

Is it possible to label a point layer with unicode numbers stored as integer within a shapefile?


texts:
I am a young developer in Africa.

I want to use OSM on the web and mobile (Android and IOS).

Is there a free cross-platform SDK at 100% and open source for this?
Or separate free SDKs for Android and IOS?


texts: ArcGIS Online - Perform analysis based on different feature services

texts:
I am creating a COVID response map and we want to show number of infected/recovered/deaths/testing per % population in each country.

The issue is that the key authoritative source is a feature layer - https://services1.arcgis.com/0MSEUqKaxRlEPj5g/arcgis/rest/services/Coronavirus_2019_nCoV_Cases/FeatureServer/2

I want to use this dynamic service and create a table for mapping that uses https://hub.arcgis.com/datasets/a21fdb46d23e4ef896f31475217cbb08_1/data - a polygon of all countries and population (from https://www.worldometers.info/world-population/population-by-country/)

Is this possible in ArcGIS Online?

The population by country can be created as a static file, the key question here is can I then divide population by the covid values dynamically or do I need to do something funky like a import into Google Sheets to do this?


texts: Pyqgis Qgis2 override Processing framework settings doesn't work

texts:
Inside my Python plugin (qgis2.18) I use some Processing algorithms.

I'm trying to override the Processing setting that allows using only the selected features. I' want to set it to False and following this answer I use that code inside my plugin before running Processing algs:

QSettings().setValue('/processing/configuration/use_selected', False)

Unfortunately after running my plugin Processing algorithms still working only with selected features.
But if I check Qgis settings I correctly see that :

Some hints about where I'm wrong?

Edit:
I noticed that inside Processing own properties the checkbox  is still checked (sorry for Italian menu) and I need to restart Qgis for applying changes.


texts: How to rename PostgreSQL/PostGIS columns while importing from WFS, using OGR?

texts:
so I'm trying to rename PostgreSQL/PostGIS columns when importing from WFS using ogr2ogr.

The post here covers how to achieve this when importing from a Shapefile, using the -sql option like this:

-sql "SELECT col_1 AS BetterName, col_2 AS ImprovedName FROM YourShapefile"

However, as I'm importing from WFS, I'm trying to do something like this:

ogr2ogr -overwrite -f PostgreSQL PG:"dbname=YourDB host=YourHost user=User password=Password active_schema=Schema" -nlt POLYGON -nln TableName "WFS:[WFS_URL_HERE]" -a_srs EPSG:4230 -sql "SELECT shape.area AS shape_area, shape.len AS shape_len FROM [What goes here?]"

How does the -sql option work in this particular use case? I have tried adding the WFS URL after the FROM clause, but this did not work.

Does anyone have any experience of this, is it even possible?


texts:
I'm running QGIS 2.6.0 and trying to process a 3x3 SAGA Grid Rank filter over multiple rasters. I've not found a lot of documentation on it, but the filter seems to allow you to select a percentile (n, where 0 < n < 100) and assigns to the central pixel of each moving window the estimated nth percentile of values in that 3x3 grid.

My rasters include null pixels. With a little experimenting I've found the filter works ok except where a raster includes a 3x3 grid with a single non-null pixel with 8 null pixels, at which point the process fails. I understand you can't estimate  percentiles based on a sample size of one, but I assumed the there would be a provision in the code for cases like this that would result in a null value for the pixel rather than crashing the whole process.

My question is whether I am missing something here so that there is a way to make this tool work, or do I need to look for an alternative solution?

If it helps, below is the log from a processing failure.

_____________________________________________

library path:   C:/PROGRA~1/QGISBR~1/apps\saga\modules
library name:   grid_filter
module name :   Rank Filter
author      :   O.Conrad (c) 2010
_____________________________________________
go...

Load grid: C:\Users\beutelt\AppData\Local\Temp\processing\a026098f4a23479ab0aa9e79a3db0f26\Test30.sgrd...

ready

okay

Parameters

Grid system: 30; 3x 3y; 2014860x -3046100y
Grid: Test30
Filtered Grid: Filtered Grid
Search Mode: Square
Radius: 1
Rank [Percent]: 20.000000

Module caused access violation!

Rank Filter: Module caused access violation!

error: executing module [Rank Filter]

C:\PROGRA~1\QGISBR~1\bin>saga_cmd io_gdal 1 -GRIDS "C:\Users\beutelt\AppData\Local\Temp\processing\54874f152788496eaa5298042ff3de9d\RESULT.tif.sgrd" -FORMAT 1 -TYPE 0 -FILE "C:\Users\beutelt\AppData\Local\Temp\processing\54874f152788496eaa5298042ff3de9d\RESULT.tif"
_____________________________________________
#####   ##   #####    #####     ###  ########
_____________________________________________

library path:   C:/PROGRA~1/QGISBR~1/apps\saga\modules
library name:   io_gdal
module name :   GDAL: Export Raster
author      :   O.Conrad (c) 2007
_____________________________________________
go...

Load grid: C:\Users\beutelt\AppData\Local\Temp\processing\54874f152788496eaa5298042ff3de9d\RESULT.tif.sgrd...

failed

error: Grid file could not be opened.

error: input file [C:\Users\beutelt\AppData\Local\Temp\processing\54874f152788496eaa5298042ff3de9d\RESULT.tif.sgrd]

error: empty input list [GRIDS]
Usage: 1 -GRIDS <str> [-FILE <str>] [-FORMAT <str>] [-TYPE <str>]
-GRIDS:<str>    Grid(s)
Grid list (input)
-FILE:<str>     File
File path
-FORMAT:<str>   Format
Choice
Available Choices:
[0] Virtual Raster
[1] GeoTIFF
[2] National Imagery Transmission Format
[3] Erdas Imagine Images (.img)
[4] ELAS
[5] In Memory Raster
[6] MS Windows Device Independent Bitmap
[7] PCIDSK Database File
[8] ILWIS Raster Map
[9] SGI Image File Format 1.0
[10] Leveller heightfield
[11] Terragen heightfield
[12] ERMapper .ers Labelled
[13] Raster Matrix Format
[14] Idrisi Raster A.1
[15] Intergraph Raster
[16] Golden Software Binary Grid (.grd)
[17] Portable Pixmap Format (netpbm)
[18] ENVI .hdr Labelled
[19] ESRI .hdr Labelled
[20] PCI .aux Labelled
[21] Vexcel MFF Raster
[22] Vexcel MFF2 (HKV) Raster
[23] VTP .bt (Binary Terrain) 1.3 Format
[24] Image Data and Analysis
[25] NOAA Vertical Datum .GTX
[26] NTv2 Datum Grid Shift
[27] ARC Digitized Raster Graphics
[28] SAGA GIS Binary Grid (.sdat)
-TYPE:<str>     Data Type
Choice
Available Choices:
[0] match input data
[1] 8 bit unsigned integer
[2] 16 bit unsigned integer
[3] 16 bit signed integer
[4] 32 bit unsigned integer
[5] 32 bit signed integer
[6] 32 bit floating point
[7] 64 bit floating point

error: executing module [GDAL: Export Raster]

C:\PROGRA~1\QGISBR~1\bin>exit


texts: How to know the required scene required to cover an area?

texts: I want to know how to find the nos. of scene of LISS 4 data required to cover an area(say 84000 sq km)? any suggestion/solution/calculation?


texts:
I would like to install sen2cor in Anaconda 4.4.0 virtual env.
How to know the list of dependencies needed? Is there environment.yaml file for sen2cor?


texts: ArcGIS Rest service is returning multiple identical records when getting KMZ

texts:
This REST URI returns 299 identical records:

http://dec.alaska.gov/arcgis/rest/services/EH/Solid_Waste_Sites/MapServer/1/query?where=1%3D1&text=&objectIds=&time=&geometry=&geometryType=esriGeometryPolygon&inSR=&spatialRel=esriSpatialRelIntersects&relationParam=&outFields=*&returnGeometry=true&returnTrueCurves=false&maxAllowableOffset=&geometryPrecision=&outSR=&returnIdsOnly=false&returnCountOnly=false&orderByFields=&groupByFieldsForStatistics=&outStatistics=&returnZ=false&returnM=false&gdbVersion=&returnDistinctValues=false&resultOffset=&resultRecordCount=&f=kmz

However if you have it return JSON or HTML, you get the distinct records:

http://dec.alaska.gov/arcgis/rest/services/EH/Solid_Waste_Sites/MapServer/1/query?where=1%3D1&text=&objectIds=&time=&geometry=&geometryType=esriGeometryPolygon&inSR=&spatialRel=esriSpatialRelIntersects&relationParam=&outFields=*&returnGeometry=true&returnTrueCurves=false&maxAllowableOffset=&geometryPrecision=&outSR=&returnIdsOnly=false&returnCountOnly=false&orderByFields=&groupByFieldsForStatistics=&outStatistics=&returnZ=false&returnM=false&gdbVersion=&returnDistinctValues=false&resultOffset=&resultRecordCount=&f=html

This was also happening with layer 0 and I somehow got it to stop by recreating the URI from the query form. I am now using the identical URI with just a different layer id and geometry type - layer 0 is point and layer 1 is polygon. I have also tried recreating the URI from the query form.

It appears that the server is using ArcGIS 10.31.

The server is back up but problem is still there. I should clarify that each record has a unique geometry, but the data fields (ObjectID, Site_ID, Comments, Site Name) are the same on every record in the KMZ. When getting the HTML or JSON the data fields are correct.

Update: I just found a second server, this one at the EPA, that is doing the same thing. The geometries are unique but the data, OBJECTID and PGM_SYS_ID is the same for every record. Here is the URI to get the KMZ:

https://geodata.epa.gov/arcgis/rest/services/OEI/FRS_INTERESTS/MapServer/27/query?where=1%3D1&text=&objectIds=&time=&geometry=&geometryType=esriGeometryEnvelope&inSR=&spatialRel=esriSpatialRelIntersects&relationParam=&outFields=PGM_SYS_ID&returnGeometry=true&maxAllowableOffset=&geometryPrecision=&outSR=&returnIdsOnly=false&returnCountOnly=false&orderByFields=&groupByFieldsForStatistics=&outStatistics=&returnZ=false&returnM=false&gdbVersion=&returnDistinctValues=false&returnTrueCurves=false&resultOffset=&resultRecordCount=&f=kmz

The EPA server also says it is v10.31.


texts: Add link to function in InfoWindow without overriding 'Zoom to' - ArcGIS

texts:
I am trying to see if there is a way to add a link in the InfoWindow that fires a function in JavaScript. I'm working with a map that shows rail lines in the US and the idea is to store selected lines in an array to collect information.

I know that you can override the 'Zoom To' link with dojo.disconnect(map.infoWindow._eventConnections[4]) to replace the zoom action with your own, but I do not want to sacrifice this function. I would use .setContent() and add the link to the function using HTML, but because I'm using .setFeatures(), this would remove information about the feature currently selected.


texts: Merging different .asc rasters located in different folders

texts:
I'm using ArcMap 10.6 and am trying to create a model with model builder to mosaic 3 rasters with the same name located in 3 different folders and repeat the process for all the remaining rasters in the three folders.If you look at the picture, you will see that the rasters are stored in 3 different folders ( named bio_a1, bio_a2 and bio_a3) and I need to mosaic the raster named bio_1.asc together, then bio_2.asc together etc...

I have built a sub model that iterates throught the folders and used collect values to create a list and set the list as a parameter (picture 1). Then I used the submodel in my main model to create a mosaic.
I have tried to run the model like that and nothing was mosaiced. So I modified the Iterator within the submodel so it would only retrieve bio_1.asc files (picture 2) and I still don't have a mosaic (picture 3).

Why doesn't this work?


texts: Assigning values to vector grids using summary statistics QGIS

texts:
In QGIS, I have a 500 sqft grid overlaid onto my map. I want to perform summary statistics (sum, mean, median, max/min) on certain vector features under the grid and assign the grid cell that value. I did this before (purple grid image below), but I can NOT remember what I did.

Here's my map with the grid overlaid. Let's say I want to calculate the total linear feet of streets in each cell. Or even the number of intersections.

Here's the example I completed a few weeks ago. BUT I CAN'T REMEMBER HOW.


texts: Can this Be Recoded to Be More Efficient with Respect to Processing Time?

texts:
I'm very new to GIS but have been asked to parse through any given ESRI GDB file-folder without using any ESRI product.  I've been trying to focus on using Fiona and Shapely as these seem to be mentioned the most as alternatives.

I am trying to determine if one particular feature in one layer of the GDB is WITHIN another feature of another layer of the SAME GDB.

As it stands right now the below sample code takes over 5.5 hours to parse through an 80MB GDB file-folder.  I was told by the GDB owner it takes ~ 5 minutes via arcpy.

Is there anything that looks obvious to a GIS veteran that I should / should not be doing that is causing execution to be unnecessarily longer?

I tried using the next() command to try and get the next layer of the GDB as I traversed the GDB, but this would crash my script.  Using Fiona's listlayers worked around this problem.

with fiona.drivers():
    #The GDB full path is:  sourceDatabaseFile
    numBLDGinPerennialWaters = 0

    for i, featureName in enumerate(fiona.listlayers(sourceDatabaseFile)):
        for j, bldgName in enumerate(fiona.listlayers(sourceDatabaseFile)):
            with fiona.open(sourceDatabaseFile, 'r', layer=i) as gdb:
                with fiona.open(sourceDatabaseFile, 'r', layer=j) as gdbStructures:

                    for feature in gdb:
                        if ((featureName == "HydrographySrf") and ('ZI024_HYP' in feature['properties']) and (feature['properties']['ZI024_HYP'] > 0)):

                            for buildings in gdbStructures:
                                if (bldgName == "StructureSrf"):

                                    if (('WID' in buildings['properties'])  and (buildings['properties']['WID'] > 0)):
                                        if shape(buildings['geometry']).within(shape(featureName['geometry'])):
                                            numBLDGinPerennialWaters = numBLDGinPerennialWaters + 1

    print "Finished parsing GDB file."


texts: What are the most effective methods for monitoring changes in polygon layers over time

texts:
I work with many different thematic areas with GIS but there is one underlying task which as yet has never been satisfactorily methodized. This is how to effectively monitor changes over time. Particularly in the areas of Landuse planning and Environmental programs there is the need to analyse how layers have changed since the last update. Non-GIS users want a simple understandable method to viusalise and document change but is often very complex in nature. Usually my method involves geoprocessing tools such as "union"  or "identify" with ArcGIS to save old and new features in the same feature class and start the analysis from there. This method is however sometimes difficult for non-GIS users to understand and then to query.
What are the most effective methods for monitoring changes using GIS?


texts: I would like start/run FME workbench from one external python script. Is there some simple syntax/command to do this?


texts: PostGIS: Ordering by the sum of two distances

texts:
I have an SQL statement that selects possible bus routes between two points, ordered by the walking distance required from the start point to the departure bus stop (distance1) and the arrival bus stop to destination (distance2).

SELECT * FROM
  (SELECT route_id,
          st_distance(journey_points,ST_SetSRID(ST_MakePoint('-4.1425658', '50.3712659'),4326)) AS distance1,
          st_distance(journey_points,ST_SetSRID(ST_MakePoint('-2.2447445', '53.4791466'),4326)) AS distance2
   FROM route
   WHERE ST_DWithin(journey_points, ST_SetSRID(ST_MakePoint('-4.1425658', '50.3712659'),4326), '3.638751848588139')
     AND ST_DWithin(journey_points, ST_SetSRID(ST_MakePoint('-2.2447445', '53.4791466'),4326), '3.638751848588139')) r
ORDER BY distance1+distance2 LIMIT 300

This works great for short distances as the ST_DWithin limits it to only a few possible routes but over long distances it is very slow. I would like to be able to use <-> to order by but it only does bounding boxes and I need it to be the sum of both distance1 and distance2. Is there any way to speed up this SQL, maybe taking advantage of indexes?


texts: OpenLayers v6.5.0: How to create a non-intersecting polygon passing through all given points

texts:
Is there any function/library that I can use to create a non-intersecting polygon?
Below is the code that I used to draw Polygon from the specified points but it is not working for the non-intersecting part.
  public drawPolygon(points: LonLat[]): Feature {
   let coordinates =  [];
   points.forEach(function (point) {
    coordinates.push(fromLonLat([point.lon, point.lat]))
   });
   let feature = new Feature(new Polygon([coordinates]));
   return feature;
  }


texts: Does Google Maps API provide statistical tools or even raster to vector conversion?

texts:
I am working with a web developer who is working with Google maps API, trying to determine what pre-processing steps I should perform (e.g. in R and GRASS), and what can be done within the google maps API. This leads me to the following questions:

Is there a way to convert raster to vector data in the google maps API?
Are there any statistical tools associated with the API?


texts: Increasing developer productivity with ArcGIS platform?

texts:
We are a small team of .NET developers. We have ample GIS experience, and none of us are new to software/database development or system administration. We have technical degrees and many years of industry experience. We have attended Esri Developer Summits.

Esri's technology - primarily ArcGIS Server, ArcSDE, and ArcObjects - plays a small but necessary role in all the software we develop. Despite ESRI's minority status in our technology stack, we spend an inordinate amount of time troubleshooting elusive bugs, crafting workarounds, deciphering its vague error messages, tracking down performance problems, and recycling processes.

Typically our problems stem from genuine bugs, poor exception handling, limiting design/architectural decisions, lack of documentation, instability, or some combination thereof. (I'm speaking of the ESRI stack here.)

From a project manager's perspective, I'm very concerned about team productivity. This costs us a lot of time. We do not have time to learn every idiosyncrasy of the ESRI stack, but we still need to get things done. (Can't live with it, can't live without it.)

What pragmatic suggestions do you have for increasing developer productivity with ESRI in the mix?

I'm not looking for suggestions on alternative technology stacks.


texts: Using QGIS Python interpreter outside QGIS

texts:
I'm writing a simple Python Geoprocessing script which leverages two specific modules:

from osgeo import gdal
import pyproj

The script must be shared with colleagues in my organisation.To avoid hooking them up with a new Python installation, I intend to leverage the QGIS python interpreter which has these two modules by default (using QGIS 3.4.3 Madeira on W10). See image below:

I'm having difficulties finding and using QGIS's Python interpreter from which my colleagues could run my script. Firstly, the normal method for identifying the interpreter doesn't point to a Python interpeter but to qgis-bin.exe instead:

Secondly, the most obvious candidate for the Python interpeter (located at c:\program files\qgis3.4\apps\Python37\python.exe) doesn't recognize the same modules:

To summarize: using Python from within QGIS provides the right modules, but using the (presumably) same Python interpreter outside QGIS fails to import the same modules. How can this be?


texts: Extract overlapped parts of the routes - Network Analysis

texts:
I am on the project to identify the specific roads that are more frequently used than others. With the network shapefiles (OS ITN), origins and destinations (OD) points, route was created for each set of OD (Closest facility tool). What I want to do next is to identify the roads that are used for more than one routes; for instance, assuming that there are 5 people, each of them departs from different point (origin) to get to the point A (destination), but there are some roads that all of them need to use along the way. In other words, I want to extract the lines that "overlapped". As an experiment, I created 3 routes separately and tried to extract overlapped features using overlay tool - intersect - but it didn't work. I tried converting it to Raster and did "Overlay sum" but didn't work.

If anyone can suggest the way to solve my problem, that would be really helpful.


texts:
I am trying to do some pre-processing operations. However I use some tools and usually I cannot execute the commands. For example, I use Fill nodata from GDAL and I got the following error:
INFO GDAL execution console output

     gdal_fillnodata [-q] [-md max_distance] [-si smooth_iterations]

      [-o name=value] [-b band]

      srcfile [-nomask] [-mask filename] [-of format] [-co name=value]* [dstfile]

Moreover, I use Fill Sinks (Wang & Liu) and again I cannot run the tool having again error.

I use the last version of QGIS 3.10.1 'A Coruna', actually I downloaded the OSGeo4W Network.

I downloaded a DEM from: https://earthexplorer.usgs.gov/
My DEM is n39_e021_1arc_v3

What is happening and how can I fix the tools?


texts: Rotate symbols perpendicular to line angle (QGIS)

texts:
I have lines representing curving streams and want to angle symbols at various points on the streams so that they are perpendicular to the angle in their particular locations.
What QGIS seems to do is give each point an angle that possibly represents the angle between the top and bottom of the stream.  I've looked into writing an expression, but there doesn't seem to be a way to include degrees or angles as components.
I'm using the LRS plugin for these layers.


texts:
I need to find the county name for a large number of addresses. I found the following API which does that.
https://geocoding.geo.census.gov/geocoder/geographies/address?form
Is it possible to load the same data from the Census API into a SQL database and run the queries locally?
Or any other data/platforms available to install locally to get county from a US address?


texts:
The only good way of creating circles in ArcMap is using true arc, (compound curves) but I need to convert those to vertex (densifying) because of a compatibility problem. I need info for geodatabase and SDO_geometry (arcsde)

and yes there is a way. I use:

UPDATE layer1 a SET arctype = 'compound' WHERE has_compound_curves(a.shape) <> 0;

and then using ET Geotools to densify  layer by layer.

But there should be easier way of doing this straightforward in Oracle or ArcMap, and I just want to densify the features that have true arc. I don't want any unnecessary vertices on straight lines and so on.

How can I identify true arcs in ArcMap, and how can I densify them as fast and easy as possible?

Same for Oracle SQL.


texts: Why Geoprocessing service show only one result?

texts:
I've searched a lot before posting this question.

I've created with model builder on ArcGIS 10.2.2 a model for clipping a point shapefile.

It works correctly on desktop, but when I publish it as a geoprocessing service and ask it to run, the result is only one point instead of the total number of features that should be clipped.

Any help or advice?