texts: Focal statistics on a mosaic dataset in ArcMap?

texts:
I frequently generate contours from lidar-derived elevation image tiles.
My workflow is:

Load the elevation image tiles into a mosaic dataset in a GDB
Run focal statistics on on the mosaic dataset
Contour the resulting focal stats raster so as to yield smoother contour lines rather than stairstep ones.

As a result, we have the elevation imagery twice and it takes extra steps.

Can we do Focal statistics on the mosaics dataset directly without generating a new intermediate image output?


texts: Export Polygon generated with Geometry Generator as dissolved polygons

texts:
Absolute QGIS newbie (3.10.8) here, so please bear with me.
I got a QGIS document where some polygons (countries) are virtually dissolved in regions (like South America), through a rule-based symbology and styled as one polygon with the geometry generator.
buffer (collect($geometry), ' ',(array_contains(array("Reg1" ,"Reg2" , "Reg3"),'Region_name'))),0)
I want to export a file (preferably .shp) of these regional objects into dissolved polygons based on the attribute "region_name" when it occurs in one of the fields "Reg1", "Reg2", "Reg3", some countries being part of several political or economical region (South America as a geographical ensemble where some of the countries belong to a specific economical institution and some other not)
How do I do that?
Edit: Value and field names, explanations


texts: Pansharpen Landsat 7 RGB using gdal_pansharpen.py

texts:
Is it possible to pansharpen Landsat 7 RGB band combination (without NIR) to have approximately the same intensities as the original RGB image with gdal_pansharpen.py (weighted Brovey algorithm)?
What is the best way to determine the coefficients?

For example, for Landsat 8 I can achieve such a result after setting weight coefficients  from here:

gdal_pansharpen.py -w 0.52 -w 0.25 -w 0.23 panchromatic.tif rgb.tif result.tif


texts:
I have a data.frame consisting of a points with an associated continuos value and their X,Y coordinates. Altough it's not stored anywhere it's projection is a Lambert conformal conic. I think this is why when I try to make it a raster object:

rast <- raster(dataframe)

I get the following error:

Error en .local(x, ...) : data are not on a regular grid

How can I tackle this issue? Thank you.


texts: GPS Pathfinder (ssf) to ESRI ArcMap (shp) Display Issue in ArcMap

texts:
I am working on a project in the Coordinate System NAD_1983_UTM_17N within ArcMap. I used a Yuma Trimble to collect data in Terrasync and made sure the projection was the same. I transferred the files to GPS Pathfinder in order to export the ssf file to a shp file in ArcMap.

I fist noticed that most of my points were being filtered out, so that the imported shp file did not match the actual area. So I exported the ssf file once more, but made sure to check off that uncorrected points were included, and the imagery matched finally.

However, even though the new shp file has the same projection, and I made sure using the Projection tools in ArcMap, the shp file appears in what must be the Arctic Circle. I am working in Florida. I traced the area myself in ArcMap using aerial imagery, but obviously the Trimble data would be more accurate.

I am at a loss of what to do, I have constantly checked the projection, made sure the same projection was listed in the layer in ArcMap, and even double checked the prj file that GPS Pathfinder references to export into ESRI.

Has anyone experience anything similar?


texts: Grouping A Single Layer Into Multiple Groups

texts:
I have a single layer that is divided up into "sub-layers" (there is only 1 layer). The "sub-layers" are shown in the Layers Panel with different names and colors based on location. I would like to group these "sub-layers" together based on their name (the way I am grouping them doesn't really matters) what matters is how I can group a single layer into groups. I would like to group these Lines below by "AK" or "AL" (again the specifics are not important). I cannot select any of these options (AK - North Slope), for lack of a better word, drag and drop does not work.

Any suggestions would be great.

I tried following information found here and here but these resources are specifically talking about multiple layers while I am working with a single layer.

I also found similar questions and this one but again they are dealing with multiple layers.

I am trying to do this using the GUI first and then through a Python script so any information you can provide will be really helpful

I am using QGIS 2.4, Postgresql 9.6, and Python 2.x


texts: Adding layer titles to summary table using batch processing and ModelBuilder

texts:
I have 180+ Landsat Raster Layers open. Part of the layer names in the TOC include the date the image was taken.

I'm trying to get weighted average NDVI for each layer in one summary table. This is the output of my model currently:
.
I can't figure out how to get the names of the rasters into a new field. Here's my process so far...

Add Landsat NDVI layers from ArcGIS online
Select by location to only include rasters in my polygon
Add selected rasters to map, creating what you see there in TOC
Run batch extract by mask tool. Here I can't customize names of output. The batch tool is awful. However, in ModelBuilder I can't iterate through these rasters because they aren't in a Raster Catalog, yet. There doesn't seem to be a way to export these directly to a Catalog without first running a tool, probably because they are huge files and from ArcGIS online?
Save each extracted raster in a Raster Catalog
Run the model which iterates through this new Raster Catalog, runs a series of calculations and gives me the output table you see.

I can very painstakingly manually add dates to the summary table. Obviously not scalable and that's why I'm here. Seems that I can't use ModelBuilder to fix the issue due to what I described above.
What can I do to get the raster names reflecting what you see in the TOC when I run the Extract by Mask tool?
I don't know any Python yet.


texts: How to include GeoPDF support for GDAL on OSX?

texts:
I have installed GDAL 2.2.1 and Poppler using homebrew

brew install gdal
brew install poppler

Everything seems to be working correctly but I cannot translate a GeoPDF to a GeoTiffs.

gdal_translate -of GTiff source.pdf dest.tiff
ERROR 4: `source.pdf' not recognized as a supported file format.

And there doesn't seem to be a driver available to it

gdalinfo --formats | grep PDF

Any idea how I would go about including GeoPDF support?


texts: Selecting current layer only using QGIS vertex tool

texts:
I am using the vertex / node tool in QGIS 3.0.2 to move vertices in a line layer. Part of the new behaviour is confusing / annoying and I don't know if anyone can tell me how to change the relevant options. I have looked in the manual but it does not help - the behaviour is not described:
https://docs.qgis.org/testing/en/docs/user_manual/working_with_vector/editing_geometry_attributes.html#id38
When I switch on this tool the cursor automatically highlights all lines that I pass. In all layers. Even if they are not visible. So if I have lines that overlap or a are close to the line of interest, the cursor switches to another line (even if it can't be seen on screen) and I need to move the cursor back to my target line in order for it to 'get focus'. This is REALLY annoying, especially as I have a vector contour layer in my project, so each hidden line gets focus even though the vertex tool doesn't actually change the vertex of a layer if it is not in edit mode (so that in itself is counter intuitive)...
There might be a good reason for this behaviour, but I'd like to switch it off and select features from only the layer I am currently trying to edit, as happened with the previous version of QGIS (2.18 for me). If it's a bug then I could report it, but I thought I would ask here first. Ah, another search of QGIS Issues has just revealed something - https://issues.qgis.org/issues/18434
It appears to be a 'feature' but maybe not all of us understand its benefits.


texts: I found a google maps api that allows me to overlay my map with a kmz/kml file however it must be hosted remotely. I have my kmz/kml files locally and would like to overlay them on my google maps, how can i get my local files to overlay my google map?


texts: Set page origin in QGIS Print Composer OR ArcMap Layout

texts:
Working on QGIS 3.6.1 I am reproducing an ArcMap Layout in Print Composer. This is basically pretty straightforward, but I need to recalculate the y-position of each map element, because in QGIS Print Composer page origin is top-left corner whereas in ArcMap Layout page origin is bottom-left (thus y_qgis = page_height - y_arcmap using the same anchor point).

Since manually doing this calculus for each map element is quite painful, my question: Is there a way to define/set page origin in one of both applications QGIS Print Composer OR ArcMap Layout?


texts: Permanently Change ArcMap's "Maximum Sample Size" Setting?

texts:
I've lately been working on feature classes with a huge number of polygons (e.g. all of the Census blocks within the Chicago metro area). I've been needing to symbolize them by quantities a fair amount to display demographic information, but because there are so many polygons I receive the Maximum sample size reached. Not all records are being used. Use this sample or change maximum sample size. message every single time.

I know that I can can manually change the sample size by pressing the Sampling... button in the Classification window, but I'm wondering if it's possible to change the default sample size somewhere, so that I will no longer have to go through this whole rigmarole. I haven't been able to find anything in Esri's documentation.


texts: ArcGIS Runtime: How to identify the topmost feature across all feature layers?

texts:
From ArcGIS Runtime Java API documentation, there is an identifyLayersAsync() method.
From documentation:

Asynchronously identifies the topmost visible geo-element in each
identifiable layer attached to the map, near given screen point.

Is there a way to for ArcGIS Runtime to return a single geo-element belonging to any layer, such that it is the topmost element?
P.S. Also asked on Stackoverflow.
Edit
I'm still having this problem. To rephrase it, I need to find out the feature (that belongs to any layer within a list of layers) that the user has clicked on. The feature returned must be feature that the user thought he/she has clicked on (based on visual judgement).
Additionally, the method must work for both maps and scenes.
I'm going to offer bounty on both sites. At this point I'm quite sure that the API doesn't do this out-of-the-box, so I'm looking for work-around solutions.
Edit 2
It looks like my question stating about "topmost feature" has caused some confusion. The "topmost" condition is only true when it is a MapView (2D). What I really need is to identify the feature that the user is actually seeing on the MapView/SceneView when the user does the clicking.
Additionally, features at the bottom layer may have a very high "z" (height) value and may appear nearest to the camera viewport.
If you call identifyLayerAsync() for a single layer, the API would be able to figure out the feature that is nearest to the viewport. However, the API does not have any mean to do the same analysis taking in all features from all layers (or selected list of layers).


texts: Graphhopper Map-Matching looses time attribute

texts:
I try to use the Graphhopper Map-Matching. My GPX source file looks something like:

<?xml version="1.0" encoding="UTF-8"?>
<gpx version="1.0" creator="GPSBabel - http://www.gpsbabel.org" xmlns="http://www.topografix.com/GPX/1/0">
  <time>2015-10-14T05:42:02.857Z</time>
  <bounds minlat="49.021236667" minlon="8.430790000" maxlat="49.023695000" maxlon="8.433425000"/>
  <trk>
    <trkseg>
      <trkpt lat="49.021240000" lon="8.431605000">
        <time>2015-08-24T05:48:08.916Z</time>
      </trkpt>
      <trkpt lat="49.021240000" lon="8.431605000">
        <time>2015-08-24T05:48:08.927Z</time>
      </trkpt>
      <trkpt lat="49.021240000" lon="8.431605000">
        <time>2015-08-24T05:48:08.938Z</time>
      </trkpt>

      <!-- much much more points -->

      <trkpt lat="49.021415000" lon="8.431866667">
        <time>2015-08-24T05:50:29.128Z</time>
      </trkpt>
    </trkseg>
  </trk>
</gpx>

Then I run:

./map-matching.sh action=import datasource=mypath/karlsruhe.osm.pbf       vehicle=car

./map-matching.sh action=match gpx=mypath/myfile.gpx

But the resulting gpx file lost all it's times.

<?xml version="1.0" encoding="UTF-8" standalone="no" ?><gpx xmlns="http://www.topografix.com/GPX/1/1" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" creator="Graphhopper" version="1.1" xmlns:gh="https://graphhopper.com/public/schema/gpx/1.1">
<metadata><copyright author="OpenStreetMap contributors"/><link href="http://graphhopper.com"><text>GraphHopper GPX</text></link><time>1970-01-01T01:00:00+01:00</time></metadata>
<trk><name>GraphHopper MapMatching</name><trkseg>
<trkpt lat="49.020944" lon="8.430902"><time>1970-01-01T01:00:00+01:00</time></trkpt>
<trkpt lat="49.021033" lon="8.431043"><time>1970-01-01T01:00:00+01:00</time></trkpt>

<!-- some more points -->

<trkpt lat="49.020944" lon="8.430902"><time>1970-01-01T01:00:00+01:00</time></trkpt></trkseg></trk></gpx>

I tried some different time formats but nothing worked. What am I doing wrong?


texts: Can ArcGIS Python API automate functionality provided in Web AppBuilder?

texts:
I was feeling adventurous the other day and I decided to dive into ArcGIS Pro and the new ArcGIS Python API.  I watched some training videos and since I've been working with ArcGIS Online a lot, I thought it would be a good idea to learn the new API.  The admin functionality seems like it's worth the effort to learn but they state it's very helpful for automating the production of apps.  I was able to automate the basic web map production but nowadays everybody is using the Web AppBuilder or they're doing all the coding through the new Server software.  It hardly seems worth the effort of automating the map production if you still have to configure the app with all the widgets.

I guess what I'm asking is can you automate workflows needed to build the sort of apps you get using the Online Web AppBuilder?


texts: Calculating number of degrees of longitude and latitude between points at same latitude on map?

texts:
I'm working on a project that use Google maps and i want to know:
Giving 2 points with same latitude but with different longitudes, how can i calculate the distance between them in terms of longitude degrees?

For example, how can we calculate how far point 1 from point 2 in longitude degrees, the latitude degrees will be 0 since they are at the same latitude line.

I know how to calculate distance in kilometers, meters.


texts: Quickly assign coordinates to region names using R?

texts:
I'm wondering if there is any quick way anyone knows to assign long / lat coordinates to a list of regions without doing so manually.

I need coordinates which are simply contained by the region so it doesn't matter where the specific coordinates fall as long as they fall into the boundaries of the area.
So for example, given the list immediately below of Austrian regions, is there any way to do this? FYI, the reason I am looking for shortcuts is that overall I have something like 400 regions across Europe so manually assigning points could take time.

Burgenland
AT12    Niederösterreich
AT13    Wien
AT21    Kärnten
AT22    Steiermark
AT31    Oberösterreich
AT32    Salzburg
AT33    Tirol
AT34    Vorarlberg

There may be no simple answer to this which I'm prepared to accept, but if anyone has any ideas please let me know?

I know about spatial joining but in this scenario I need coordinates first to be able to perform other sorts of analysis. At this point, I just need the coordinates. I use R Studio normally.


texts: Rasterize shapefile to a new gtiff file with gdal_rasterize utility,however,I cannot open the tif file?

texts:
Now I'm tring to rasterize a shapefile with a utility "gdal_rasterize" and want to generate a new gtiff file,but when I open the rasterized gtiff file,it prompts "Windows Photo Viewer can't open this picture because the file appears to be damaged,corrupted,or is too large.

.........................................................................................
.........................................................................................
Command Line:gdal_rasterize -burn 255 -burn 0 -burn 0 -ts 800 800 -l Countries02 "D:\Map Data\Countries02.shp" "D:\Result\Countries02.tif"

Map Data:
https://skydrive.live.com/redir?resid=7EF4FFBCD523D8E0!431&authkey=!AKiLt1a0bpOGi7w

The output gdalinfo returns is as follow:
Driver: GTiff/GeoTIFF
Files: d:\Result\World.tif
Size is 800, 800
Coordinate System is `'
Origin = (-180.000000000000000,83.623596191406250)
Pixel Size = (0.450000000000000,-0.217029495239258)
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (-180.0000000,  83.6235962)
Lower Left  (-180.0000000, -90.0000000)
Upper Right ( 180.0000000,  83.6235962)
Lower Right ( 180.0000000, -90.0000000)
Center      (   0.0000000,  -3.1882019)
Band 1 Block=800x1 Type=Float64, ColorInterp=Gray
  Minimum=0.000, Maximum=255.000, Mean=87.551, StdDev=121.080
  Metadata:
    STATISTICS_MAXIMUM=255
    STATISTICS_MEAN=87.5514609375
    STATISTICS_MINIMUM=0
    STATISTICS_STDDEV=121.07999102565
Band 2 Block=800x1 Type=Float64, ColorInterp=Undefined
  Minimum=0.000, Maximum=0.000, Mean=0.000, StdDev=0.000
  Metadata:
    STATISTICS_MAXIMUM=0
    STATISTICS_MEAN=0
    STATISTICS_MINIMUM=0
    STATISTICS_STDDEV=0
Band 3 Block=800x1 Type=Float64, ColorInterp=Undefined
  Minimum=0.000, Maximum=0.000, Mean=0.000, StdDev=0.000
  Metadata:
    STATISTICS_MAXIMUM=0
    STATISTICS_MEAN=0
    STATISTICS_MINIMUM=0
    STATISTICS_STDDEV=0


texts:
I must be overlooking something really simple. I want to use a QComboxBox dropdown list to control the form input interpretation, rather than provide a choice of input values. Maybe this example would help clarify:

The dropdown list would specify what weight units will be entered by the user. In other words, if the 'kilograms' list value were chosen, a subsequent entry of '10' by the user in the textedit area would be interpreted as 10 kilograms. As far as I can tell, there is no way to bind a QComboBox in this way. Suggestions?


texts: Looking for Web service standard to disseminate water level data

texts:
I’m about to lead a project that has to disseminate observed, forecasted and predicted water level data coming from tide gauges. I did a lot of research to identify which dissemination standards to support (Web Service). Many of the recommendations are pointing for the OGC SOS but it seems that it only does the observation dissemination part of our needs. No way to embed forecast and prediction and the transactional aspects are not adapted to support QA/QC advanced process.

So, I’m wondering if there are standards that exist or that are on their way to be adopted by OGC (such as SOS) but concerning the Water levels time series (observation, forecast and prediction). We have already looked at the SensorThingAPI. It could be an option but it is still not adopted by the community (at the time that I write this post).

Does anybody is aware of a Web service standard that can meet our needs?


texts: Corporate logo added to openlayer web mab created with Qgis2web

texts:
I have used the advice from a similar post:

Megacorp Industries

This works to put in a title at the top but it doesn't tell me how to float a corporate logo in the bottom right corner (for example) with some control over position and size.
It seems simple but I can't seem to find the answer.
I want to be able to stick the logo png in the "image" directory created by Qgis2web and insert it into my html from there.


texts: How to locate where shapefiles do NOT intersect? Two road datasets

texts:
I am trying to match a Canadian road dataset with a US road dataset in ArcGIS 10.1. One way to assess how complete the "matching process" is is by seeing where the roads intersect and where they do NOT intersect (see picture; US roads in red, Canada roads in blue, US-Canada Border in black).

In the picture above, I have three areas where the lines intersect and cross the border and one area where the US road reaches the border but does not have a connection on the other side. Presumably, the one US road (in red) does not just stop at the US-Canada border (in black), but should continue to a Canadian road (in blue). To see where the roads intersect, I use the Intersect tool to create points where the shapefiles intersect. However, I do not have a method to assessing where the US and Canadian roads do NOT intersect.

Is there a way to locate where the US roads and US-Canada border intersect, but where Canadian roads do NOT intersect? Or vice versa with Canadian roads?

This would locate where the matching of the two datasets fails, an important thing to assess for my goals.


texts:
I am using Boundless Server 1.0.2(2.12-server-1.0.2) on Cent OS 7

My Geowebcache is located in following directory & My data directory located in EFS storage

/var/opt/boundless/server/geoserver/tilecache

I want to move tilecache folder into EFS storage & edited web.xml file as per this link then restarted tomcat server but it unable to move its EFS storage

-Clustering is used for this Geoserver

-Tried to edit 'geowebcache.xml' file but unable to update its storage

-Tried same on another standalone Geoserver in my PC (on windows platform),it works perfectly

-Wondering why its not working for Boundless Geoserver.

-Researching on Web


texts: Lines from Oracle Spatial Layer in QGIS are not completly visible

texts:
I connected a spatial layer from the Oracle database which contains lines (OracleDB objectclass typ is polyline) to be visible in QGIS.

An example geometry is as follows:

sdo_geometry(2002,2056,null,
    sdo_elem_info_array(1,4,3,1,2,1,7,2,2,11,2,1),
    sdo_ordinate_array(2610857.31245004,1266613.32122543,2610857.70740207,1266610.23735626,2610858.72627747,1266602.3016907,2610859.88613344,1266593.26306758,2610859.84612155,1266592.56309404,2610859.5561095,1266591.92511322,2610859.27210097,1266591.49512447))

In QGIS only a part of the above line is visible. I know that this line is not only straight. So I was wondering, if there is kind of an "arc-problem" behind this and if there is any known solution to the problem?

I didn't find anything in the official QGIS documentation.


texts:
cnn.Open "Provider=Microsoft.Jet.OLEDB.4.0;" & _
         "Data Source=C:\Workspace\Data.gdb"
rst.Open "SELECT DISTINCT [Adi],[KimlikNo]FROM MahalleSiniri ORDER BY [Adi];", _
         cnn, adOpenStatic
rst.MoveFirst
    ComboBox4.Clear
    With Me.ComboBox1
    .Clear
   ComboBox1.Text = "Mahalle Seçiniz..."
    Do
        ComboBox1.AddItem rst![Adi]
        ComboBox4.AddItem rst![KimlikNo]
        rst.MoveNext
    Loop Until rst.EOF
End With
'UserForm_Initialize_Exit:'
On Error Resume Next
rst.Close
cnn.Close
Exit Sub
UserForm_Initialize_Err:
MsgBox Err.Number & vbCrLf & Err.Description, vbCritical, "Error!"
Resume UserForm_Initialize_Exit
End Sub

In this code block, I can add *.mdb database to my interface and it works fine. But my data contains raster data and I want it to connect to a *.gdb database file.

How can I correct it? I tried several solutions I found on the internet but they didn't solve my problem.


texts:
I'd like to develop a simple spatial system with database for plotting and managing geographic features.I thought I should install on an ubuntu server:

qgis server
postgreSQL
postGIS

and write python qt modules to access the data from android devices and desktop PCs.

Could helpful people let me know if I'm on the right track?
Thanks
Mini


texts: Import arcpy using PyScripter or IDLE gives No module named RuntimeError?

texts:
I'm trying to figure out this problem. Yesterday I installed PyScripter and since then, scripts doesn't work. When I run any script (in PyScripter or IDLE) which uses "import arcpy", I get this error message:

import arcpy

Traceback (most recent call last):
  File "<pyshell#0>", line 1, in <module>
    import arcpy
  File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\__init__.py", line 17, in     <module>
    from geoprocessing import gp
  File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing\__init__.py", line 14, in <module>
from _base import *
  File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing\_base.py", line 14, in <module>
    import arcgisscripting
  File "C:\Python26\ArcGIS10.0\lib\arcgisscripting.py", line 81, in <module>
    from RuntimeError import RuntimeError
ImportError: No module named RuntimeError

Are there any suggestions on a potential to fix?


texts:
Is it possible to modify the default fonts to be used from "MS Shell Dlg 2" to another?

If yes, where do I have to set it?


texts: View database connections in QGIS with Python

texts:
I am designing a GUI for a QGIS plugin I am developing, which is intended to plot information from a MySQL database. In the GUI, I would like to have a combo box of all the database connections saved in QGIS.

These can be viewed by printing allKeys() for QSettings() in PyQt4.QtCore. When this is done, the keys I am interested in are:

...
MySQL/connections/adsb/database
MySQL/connections/adsb/host
MySQL/connections/adsb/password
MySQL/connections/adsb/port
MySQL/connections/adsb/save
MySQL/connections/adsb/username
MySQL/connections/database
MySQL/connections/host
MySQL/connections/password
MySQL/connections/port
MySQL/connections/save
MySQL/connections/selected
MySQL/connections/strath_adsb/database
MySQL/connections/strath_adsb/host
MySQL/connections/strath_adsb/password
MySQL/connections/strath_adsb/port
MySQL/connections/strath_adsb/save
MySQL/connections/strath_adsb/username
...

Where 'adsb' and 'strath_adsb' are databases. How would it be possible to read the names of all connected databases and store them in an array, so that they could then be displayed as options in the combo box?


texts: Why do each of my buffer-polygons have different areas

texts:
I`m using Qgis 2.2 and 2.4 and loaded a wms-map using the coordinate system ERTS89 UTM 32N. I set 8 points (all in on shp-fileand in UTM 32N) on various locations on the map. After that i created a buffer for each of the points with 72 and later 99 cricle-segments, because i thought this is more accurate according to a real circle. Each buffer has a radius of 50m.Then i calculated the area of this buffer-polygons using the fieldcalculator.

Now the problem.

Not only that the calculated areas differ from a mathematic calculated circle (with a radius of 50m, the area should be around 7853m² --> fieldcalc gives areas around 7830m²). The areas of each of 8 polygons differ between each other. And the more circle-segments i use, the more deviation between the results, alldough the difference between calculated area and mathematical area gets smaller. WHY does that deviation exist, alldough each polygon should have the same radius and therefore the same area.

I understand the difference between qgis-calculated area and mathematic area, but not the difference between the calculated ones. I tried to use different coordinate-systems like WGS84, Nad83 (NSRS2007), Lambert, DHDN GK4,......
All with the same problem.

Is it a qgis-discrepancy i have simply to deal with or did i something wrong?


texts: Reading empty feature table with FME / Python?

texts:
I would like to convert a feature table Geodatabase (empty) to MapInfo, but my feature table is empty.

I would like to write the schema of feature table in MapInfo file.

Is it possible to force the writer MapInfo ?


texts: leaflet map crs is 3857 but coordinates 4326

texts:
I am developing my first leaflet applications. So this is my trial application demo.

       var osmUrl = 'http://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png',
            osm = L.tileLayer(osmUrl, {
                maxZoom: 18,

            });

        var map = L.map('map').setView([51.505, -0.09], 12).addLayer(osm);

I am getting the CRS of the map, and CRS is EPSG:3857.

 console.log("map crs: " + map.options.crs.code)

I can add a point data with EPSG:4326 coordinates.

L.marker([51.5, -0.09]).addTo(map)
.bindPopup('A popup.')
.openPopup();

But I can not add EPSG:3857 coordinates.

    L.marker([5732953.77585,-10018.7582914]).addTo(map)
    .bindPopup('A popup.')
    .openPopup();

This goes wrong place.


texts:
I have multiple wind data at varying heights in .asc-file format.

I want to interpolate them in QGSI but it seems I can't even display the points even if I click add raster. It's still on a table.

How can I achieve points out of .asc-file?


texts: Multipatch georeferencing in ArcGIS Pro?

texts:
I am currently using ArcGIS Pro to create a 3D campus model and I have floor plans extruded by each individual floor to represent buildings. Another group is working on the modeling side of things and have multipatch features (.shp) that they have created, but there are a lot of issues that arise when trying to bring them into ArcGIS Pro. First they are not georeferenced, so we have to find them wherever they end up (usually at 0,0 in the center of the globe) and then when we get them to fit (this is done manually).

Does anyone know of a more automatic way to get these multipatches georeferenced with a more automatic method?


texts: Delete a Feature from a Different Version While Editing

texts:
We have one feature class that always loads from the default version (named 'Work Request'), even when you are editing inside a different version.  We use this feature class to identify areas where work in ongoing and we pull it from default so it doesn't get outdated when it displays while working in a version.

My question is: Is there a way to delete features from Work Request while editing in a different version (that isn't default)?

Here's a screenshot showing the version I am editing inside of (DN_1857) and the 'Work Request' feature class that is from Default:

The problem I am trying to solve is that it's part of my job to open a version, verify that it's correct, then reconcile and post it to default.  However, this process does not delete the work request feature that was identifying the work that I've now completed and posted.  So I have to start a new editing session just to delete the single polygon.  I would like to be able to delete that polygon while I have the version open and save myself the extra steps.


texts: Display InSAR calculation results on satellite images

texts:
I have used PSInSAR technology to process and got the result of ground subsidence velocity. I saw in [a paper]1 that the results of ground subsidence velocity (a layer) can be displayed on satellite imagery（the figure as following）.

I plan to use Sentinel-2 remote sensing images, and the software is QGIS, but I don’t know how to do it, so some processing tutorials or ideas are needed here.


texts: Exporting only selected items in My Places to another KML File in Google Earth Pro?

texts:
It seems that when I right click and export from my places in google earth pro, it will export the entire kml instead of the selected or checked items within it.

How can I just simply export the checked items to another kml?

Even when I import the KML into arcmap, there is no indication of what was selected etc. IF that was the case, I could do a query.


texts:
I installed ArcGIS 10 SP1 last week.  My graphics card meets the spec and I'm up to date on drivers.  I have a layer of parcels that are tint-banded in the MXD.
At some scales, straight lines that are the same color as the tint band shoot across the map.  After this happens, the next, second, or third extent change crashes ArcMap (usually it's the 1st or 2nd).  When I successfully change extent, I'll get a different set of straight lines.
In the image below, only two of the polygons are (or should be) tint-banded:

Has anyone else experienced this and know the cause? And does anyone have a solution (besides not using tint bands)?

I did a check/repair geometry, but the issue continues.
I was able to replicate the issue with US Counties: here's a layer package of the US Counties with a green tint band.  As I zoom in/out, I experience the same issue.

Perhaps this is a result of the way I create the tint bands.


texts: I tried to download and install an older version of QGIS, but I don't understand how it works... I downloaded and opened this file: qgis-2.0.0.tar.bz2, but I don't see any program file here. I also tried to download the two .exe-files for this version, but I couldn't open them on my mac. What do I do wrong? I already have QGIS 2.18 installed.


texts: Generating anchor points using table in ArcGIS for Desktop?

texts:
I have geocoded multiple points based on GPS data over 7 days. Each point has a duration, speed (see table and arcmap image). I would like to create anchor points from this file, that is merge all points that are stationary in the same location, across all points e.g. Home as anchor (which have geocoded in another file), shops, school etc.

Could anyone suggest how I might achieve this?

I would like to also measure the frequency of trips to the same destination and the durations of these trips.


texts: Creating an address locator/geocoding addresses in ArcMap without using City or Placename

texts: I am creating an address locator, and I have it working, but if I want to search for any address, I have to enter the entire address into the search. I want to be able to enter everything besides the City or Placename, but I can't seem to make it work. 95% of my addresses are rural addresses located within a Canadian municipality, so the City or Placename field is superfluous, since we use a grid system for our road network.


texts: Creating multiple lines from one Array using ArcPy in ArcGIS Pro

texts:
I am trying to create multiple lines from a single array. What I want to do is have a single destination point and have all the other point make lines to this destinations point. I am using arrays to create polylines, however I am unsure how to create a loop which makes polylines using one point and then going through the array and creating a new line for each point in the array.
I am using python 3 in ArcGIS Pro.
Io illustrate:
Points Array : <Array [<Point (1,2), <Point (3,4)>,<Point (5,6)>, <Point (7,8)>]>
I want to make the lines using these points:
(1,2),(3,4)
(1,2), (5,6)
(1,2), (7,8)


texts: Is there a statistical test to calculate difference/similarity between two sets of lat/lng data points?

texts:
I have two sets of data points where each data point is a latitude-longitude combination. The two sets represent two different time frames, let's say 2010-2014 and 2015-2018. The number of data points within the first time frame does not equal the number of data points in the second time frame. The data points are also independent, meaning that there exists no mapping from time frame 1 to time frame 2.

I would like to know if there exists a statistical test to compare the two sets of data points and see if the data points have shifted in geographical location, or stayed the same?

For instance, let's say the bulk of data points from 2010-2014 lies within the US, and within time frame 2015-2018 the majority of points is still in the US but also in Canada. Can we statistically say something about the shifting spatial characteristics of the two time frames.

I have tried figuring out a way to do this (PySAL for Python) but I seem to be lost.


texts: Populate field using a combination matrix

texts:
I have a polygon shapefile that represents combinations of soil hydrologic type (A,B,C,D) and soil cover (53 different possible covers, coded from 0 to 53). Each of the 53 soil covers can be either soil hydrologic type "A", "B", "C" or "D". Depending on this, each intersection will have assigned a number, i.e. field "CNII". I would like to populate CNII field automatically based on a table that looks like this:

My attribute table containing fields "Cod_Cob" (Integer) and "Clase_Hid" (Text) would need to have "CNII" (Integer) populated as a function of the former two:

I am new to Python parse in the field calculator (I am using Arcmap 10.4.1 with Standard License) but would like help finding a solution using this tool or any other kind you think is appropriate. I checked how to write if else statements and return the values input manually in the code but this doesn't seem adequate because there are more than 200.000 possible combinations for the soil type-soil cover complexes. I have 85 features I need to populate and I can do it manually by now, but I will be needing this solution for bigger projects.


texts: WPS not accepting https in process input

texts:
The following WPS request successfully styles a png when the png is served over http.

<wps:Execute version="1.0.0" service="WPS" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns="http://www.opengis.net/wps/1.0.0" xmlns:wfs="http://www.opengis.net/wfs" xmlns:wps="http://www.opengis.net/wps/1.0.0" xmlns:ows="http://www.opengis.net/ows/1.1" xmlns:gml="http://www.opengis.net/gml" xmlns:ogc="http://www.opengis.net/ogc" xmlns:wcs="http://www.opengis.net/wcs/2.0" xmlns:xlink="http://www.w3.org/1999/xlink" xsi:schemaLocation="http://www.opengis.net/wps/1.0.0 http://schemas.opengis.net/wps/1.0.0/wpsAll.xsd">
    <ows:Identifier>ras:StyleCoverage</ows:Identifier>
    <wps:DataInputs>
        <wps:Input>
            <ows:Identifier>coverage</ows:Identifier>
            <wps:Reference mimeType="image/${fileFormat}" xlink:href="${process.env.PROTOCOL}://${process.env.SERVICES_HOST}:${process.env.PORT}/${filename}" method="GET"/>
        </wps:Input>
        <wps:Input>
            <ows:Identifier>style</ows:Identifier>
            <wps:Reference mimeType="text/xml; subtype=sld/1.1.1" xlink:href="${process.env.PROTOCOL}://${process.env.GEOSERVER_HOST}/geoserver/wms?request=GetStyles&amp;layers=${layerName}&amp;service=wms&amp;version=1.1.1" method="GET"/>
        </wps:Input>
    </wps:DataInputs>
    <wps:ResponseForm>
        <wps:RawDataOutput mimeType="image/${fileFormat}">
            <ows:Identifier>result</ows:Identifier>
        </wps:RawDataOutput>
    </wps:ResponseForm>
</wps:Execute>

However when I change the png to be served over https (${process.env.PROTOCOL} === https) I get the following response from geoserver:

<ows:ExceptionText>Failed to retrieve value for input coverage
    sun.security.validator.ValidatorException: PKIX path building failed: sun.security.provider.certpath.SunCertPathBuilderException: unable to find valid certification path to requested target
    PKIX path building failed: sun.security.provider.certpath.SunCertPathBuilderException: unable to find valid certification path to requested target
    unable to find valid certification path to requested target</ows:ExceptionText></ows:Exception></ows:ExceptionReport></wps:ProcessFailed></wps:Status></wps:ExecuteResponse>

Any ideas if https is supported as a process input? I don't see it listed here.


texts: ArcGIS API For Javascript - Custom DynamicMapServiceLayer Image Being Distorted In Map

texts:
I am having some trouble with a custom DynamicMapServiceLayer that I hope someone can either help me out with, or explain what is happening...

On the application I am currently working on, there is a requirement for a 'HeatMap' style layer.  Since I have to support all the way back to IE 7 - using HTML5 canvas solutions are not possible.  For that reason, I have created a Server Object Extension that renders that heat map at server, and returns an image that should work in any browser.

The SOE is working great.  However, the image that gets added to map is a bit 'smooshed'.

But the raw image (same request as above - just grabbed this out of Google Developer Tools Network tab) getting returned by the SOE does not show this distortion.

Why is Javascript API doing this - I thought it would just display the raw image?


texts:
In one of my plugins I would like to use the distance geoprocessing algorithm. When running this algorithm from the menu Vector->Analysis Tools->Distance Matrix I get a different output than running the same function from the Geoprocessing Tools. (I am interested in the k*n matrix)

Unfortunately the output from the Geoprocessing Toolbox has different headers in the output and I need the result from the process run with the vector-analysis tools.

The first row of the processing output is

ID,DIST_1,DIST_2, ...

The Output of the menu-version is:

ID,4711,0815, ...

4711 and 0815 are my idnetifiers choosen in the dialog.

Is there a possibility to run this tool instead of

processing.runalg("qgis:distancematrix" ...

in my python plugin?


texts: Finding difference between two rasters using ArcGIS Spatial Analyst?

texts: I know how to find difference between two rasters with Spatial Analyst tools, but how to display attributes, for example differences between two Corine Land (2006 - 2012) rasters, how to see codes of Corine Land that change from year 2006 to 2012.


texts: drawing on AGSSketchGraphicsLayer leads to memory issue

texts:
In my app I have  a functionality to the log the GPS. I am using AGSSketch GraphicLayer to do the drawing with the points received from GPS/location manager.

I see that the memory allocations get mounted as the number of locations get increased. I used the instruments tool (comes along with Xcode) to know which objects leads to it. The tool shows that sketch layer is causing the jump in memory allocations.

Is there a way to avoid this. Few blogs suggested to update the ArcGIS sdk but that didn't help.


texts:
How do I change or add the CSS style for panZoomBar?

I.e. to change the panup, pandown, panleft, etc hover in CSS by attaching the class in CSS for those...  and which class should I create it in...


texts: Mapping NetCDF CF Mercator parameters to GDAL

texts:
I have coordinates on a mercator grid projection from a NetCDF (CF) file. The metadata attributes in the file include : standard_parallel, longitude_of_projection_origin, earth_shape (spherical), earth_radius.

I want to use the gdal python API to transform the coordinates into lat-long. I need to set up the mercator coordinate system which I presume is best accomplished using SetMercator2SP. It takes arguments StdP1, CenterLat, CenterLong, FalseEasting and FalseNorthing.

How do I map these to the data I have?

I assume StdP1 is the standard parallel.

Is the longitude of projection origin the CenterLong?

What about CenterLat ... is that typically assumed to be zero if not listed?

Do I ignore false easting and northing?


texts: Leaflet ImageOverlay Image does not appear

texts:
I'm trying to stitch the picture taken from my drone onto the map, and I'm using Django 1.11.14 and latest version leaflet 1.3.4. I tried to use ImageOverlay, there is no error but the image does not appear on my map

var imageUrl = "/media/photos/2018/10/18/DJI_0001_2.JPG";
var imageBounds = [[101.6999678,2.9244267], [101.611, 2.983226]];
L.imageOverlay(imageUrl, imageBounds).addTo(map);
map.fitBounds(imageBounds);


texts: select intersection of 2 geometries, polygon and polyline

texts:
I'm searching how to have a result of 2 geometries (villages and wetland) which don't intersect.

When I write:

select distinct v.name
from villages v , wetland w
where ST_intersects(w.geom,v.geom) ='f' ;

it selects all the names. I want just the names which don't intersect.

When I write:

select distinct v.name
from villages v , wetland w
where ST_intersects(w.geom,v.geom) ='t' ;

it selects 231 lines.


texts: Tile server that outputs .svg tiles instead of raster?

texts:
Is there one? I found all major suppliers work with raster output (png) not vector (svg).

I wonder why.


texts: How GIS and SDI (spatial data infrastructure) are related?

texts:
I read the definitions of both on wikipedia but they didn't mention any relationship between the two terms although there must be one because they are dealing with the same thing: spatial data.
Can we consider that regional GISs are building blocks of a national SDI?


texts: Leaflet routing machine waypoints with numbers

texts:
Is it possible to display numbers near waypoints on map? I create multiple routes in loop and passing waypoints via array. Routes get painted with Leaflet Routing Machine. I want to show number of waypoint near it on map. Here's my code

for(var i = 0; i < data.length; i++){
   var waypoints = [];
   for(var j = 0; j < data[i].length; j++){
       fs.push(L.latLng(coords[data[i][j]][0], coords[data[i][j]][1]));
   }
   if (fs.length > 2){
   L.Routing.control({
           waypoints:waypoints,
           lineOptions:{styles:[{color: colorArray[i]}]},
           show: false,
       }).addTo(map);
   }
}

data is array that contain orders of points to visit, coords contain coordinates of these points. By iterating through these 2 loops I create several routes. Is it possible to add something to L.Routing.control and display number of the point near it's marker? What can I use for it if it's not possible just with leaflet routing machine


texts: Does anyone know of a project using a raspberry pi and gis that would work for a research project for a computer science major? I am looking to do a research project involving a raspberry pi and incorporating gis into it. I am also looking to use python but all ideas are welcome. Very broad but looking for any ideas. Possibly a mapping application.


texts: SLD polygon styling + labeling in Geoserver

texts:
I'm trying to label a polygon alongside its boundary. I found out how to do it but the stroke color of the polygon looks wear because of a black line above it. I think that the step I use to label the boundary of the polygon instead of the polygon itself generate also a small polyline which I cannot manipulate. I'd like to make this black polyline disappear or at least to give it a color like my polygon stroke.
    <Rule>
      <Name>Zone 100m</Name>
      <Title>Zone 100m</Title>
     <ogc:Filter>
        <ogc:PropertyIsEqualTo>
          <ogc:PropertyName>ToBufDist</ogc:PropertyName>
          <ogc:Literal>100</ogc:Literal>
        </ogc:PropertyIsEqualTo>
      </ogc:Filter>
      <PolygonSymbolizer>
        <Stroke>
          <CssParameter name="stroke">#FF0000</CssParameter>
          <CssParameter name="stroke-width">4</CssParameter>
        </Stroke>
      </PolygonSymbolizer>
      <TextSymbolizer>
      <Geometry>
        <ogc:Function name="boundary">
        <ogc:PropertyName>the_geom</ogc:PropertyName>
        </ogc:Function>
      </Geometry>
        <Label>
        <ogc:PropertyName>Project</ogc:PropertyName>
        </Label>
        <LabelPlacement>
            <LinePlacement/>
        </LabelPlacement>
         <Halo>
            <Radius>3</Radius>
            <Fill>
                <CssParameter name="fill">#FFFFFF</CssParameter>
            </Fill>
        </Halo>
        <Fill>
            <CssParameter name="fill">#000000</CssParameter>
        </Fill>
        <VendorOption name="followLine">true</VendorOption>
      </TextSymbolizer>
    </Rule>

Sorry, forgot to post a picture


texts: Calculating cost, in time, of walking on sand with ArcMap?

texts:
Using ArcMap 10.2, I'm trying to calculate the path distance around some archaeological sites and later the least cost path to the nearest water source. I'm interested to see how far the site's inhabitants could have walked in a given time. In the path-distance tool, there is an option to add a horizontal and/or vertical factor for a more accurate cost of travel time on different terrains and slopes.Most of the area is comfortable to walk so I don't want a complicated surface factor, only to add one area of increased time-cost because of sand dunes.  I managed to add the cost of the slope as the vertical factor using Tobler's (1993) method, but can't find how to determine the relative increase in time-cost of walking on sand, as the horizontal factor. I only found a paper discussing the increase in energetic cost: http://jeb.biologists.org/content/jexbio/201/13/2071.full.pdf
According to the authors, the cost of walking on sand is 2.1-2.7 times more energy consuming than walking on a hard surface (at the same speed).
So far I haven't tried to add a cost for the sandy areas because I don't know what the factor should be. or put another way, I don't know how much more time it would realistically take to cross a sandy cell, compared to a non-sandy cell.

Does anyone have an idea how to calculate the difference in time, if say, the energy was the same?

Does anyone know of previous publications where this is mentioned?


texts: Errors while building grass-7.0.1 with Msys/MinGW

texts:
I need to build grass-7.0.1 for qgis on Windows 10 with Msys/MinGW. I succefully installed all packages from general and several packages from optional requirements. Here is my configure options (minimal):

./configure \
--prefix=/usr/local \
--enable-shared \
--disable-static \
--enable-largefile \
--disable-x11 \
--with-cxx \
--with-proj-share=/usr/local/share/proj \
--without-gdal \
--without-x \
--without-sqlite \
--without-postgres \
--disable-doc \
--without-freetype \
--without-cairo \
--without-tiff \
--without-png \
--without-opengl \
--without-fftw \
--without-regex

Then I wait about 15 minutes while compiling and get this message:

In case of errors please change into the directory with error and run 'make'.
If you get multiple errors, you need to deal with them in the order they
appear in the error log. If you get an error building a library, you will
also get errors from anything which uses the library.

Grass Makefiles contains complex referencing and it's too expensive to try fix the problems one by one. Maybe any ideas?

UPD

Ok, after some experimentation, I realized that this is a problem of missing demolocation. But that ticket is outdated, so I still need help. This error occurs during the compilation of all these packages:

access: No such file or directory
ERROR: LOCATION </dist.i686-pc-mingw32/demolocation> not available


texts: Issue On Having HomeButton and Overview Together on ArcGIS Sever JavaScript API

texts:
Using ArcGIS 10.1 and ArcGIS Server Javascript Api I am having issue on having Home Button and Map Overview on the map. for what ever reason the home button is not showing up the map
Here is the Script which I have is:

<script>
  require([
    "esri/map",
    "esri/dijit/OverviewMap", "dojo/parser","dijit/layout/BorderContainer", "dijit/layout/ContentPane",
    "esri/dijit/HomeButton", "esri/layers/FeatureLayer", "dojo/dom-construct", "dojo/domReady!"
  ], function(
      Map, OverviewMap,
  parser, HomeButton, FeatureLayer
     ) {
    parser.parse();
    var map = new Map("mapDiv", {
      basemap: "topo",
      center: [-126.416, 55.781],
      zoom: 6
    });
  //===================================================== Overview
    var overviewMapDijit = new OverviewMap({
      map: map,
      visible: true
    });
    overviewMapDijit.startup();

 //===================================================== Shapefile
    //add a layer to the map
    var featureLayer = new FeatureLayer("http://somewhere/1", {
      mode: FeatureLayer.MODE_ONDEMAND,
    });
 map.addLayer(featureLayer);

 //===================================================== Home Button
var home = new HomeButton({
map: map
}, "HomeButton");
 home.startup();

  });
 </script>

and CSS for Home Button is:

#HomeButton { position: absolute;  top: 120px; left: 50px;  z-index: 2; }

Can you please let me know why this is happening?


texts: Do you know any geoportals which are fully (or almost) compatible with INSPIRE directive?


texts:
I've got a shape file that I loaded into PostGIS that has both road and parcel multipolygons and I need to remove the roads. I tried loading the shp file in qgis and it loads them both as a single layer. I need to remove roads as best I can without removing parcels. I can remove most of them by testing whether the address is equal to null. However this deletes parcels where that data is unavailable.

My idea is to create another filter that tests the complexity of the geometry since most roads are complex polygons and parcels are generally simple (not all). Then I think I will find a "close enough" solution where most roads are removed, and little to no parcels are removed.

Here is what I've looked at to no avial:

st_isSimple : tests whether a polygon is anomalous (self intersection or tangent)

st_numPoints : Does not work because it only looks at the first linestring and there are no linestrings

Some others...

Anyway, what ideas do you have for classifying a polygon as a road? Is there a query I can do? A feature I can extract offline?

-N


texts:
I'm trying to write a processing script using grass r.neighbors algorithm for rasters. When you use the Gui for this algorithm, you can leave the GRASS layer extent field blanc to use the minimal extent. On the contrary, when used in the python console, I can't find how to input this kind of "default" value.

I saw on https://docs.qgis.org/2.2/en/docs/user_manual/processing/console.html that it seems to be mandatory to input the four coordinates for an extent parameter.
 I tried using iface.mapCanvas().extent() but the object returned is not accepted as a parameter and I can't find how to get its coordinates..

Is there a way to get the xmin, ymin, xmax, ymax of the layer in PyQGIS to input them in the algorithm? Am I doing it totally wrong ?


texts: Extracting/processing raster information from .lyr or .kmz file

texts:
I have a raster which is stored as a .lyr file, which I got from the following source:

http://forestrymaps.alaska.gov/arcgis/rest/services/Rasters/AK30mSlope/MapServer

When I open it in ArcMap I cannot seem to perform any Raster based analysis, as I get an invalid GP data type error - and I can't even open it in QGIS.
I am essentially trying to convert it to a raster format in order to perform raster statistics (zonal statistics in ArcMap or raster <> polygon statistics in QGIS), however neither program recognises it as a valid raster, perhaps because it is lyr file.
The link above also allows you to generate it as a KMZ file, however I haven't had any luck trying to utilise this or the contained KML file either. The link does let you export an JPEG or PNG file directly, however after playing around with the parameters I could not get the output to the desired level of quality/resolution.

How might I use this data for raster analysis?

In it's current format it seems I can open it in ArcMap to view it, but not much else.

I also tried to find the data source the .lyr file is pointing to using a python script I found in my searching, however the file does not seem to be a feature class, apparently, so I could not get the source.


texts: How to extract various Thresholds of Chl-a pixels in ArcGIS in bulk?

texts: I have 17 years worth of Chl-a data from January to December in raster format and now I just want to just look at amounts from 10-30 mg/m3 of Chl-a data. Can I extract just this data from the image, keeping the shapefile and can I do it in bulk with file names being differientiated as well?


texts: Hiding border of dot density display in ArcGIS Pro?

texts:
I have a dot density symbology but I do not want to see the polygon borders. Using dot density layers makes it seem that ArcMap provided the ability to turn off the borders but I can't figure out how to in ArcGIS Pro.

Any ideas how to just show the dots, without the polygon borders?


texts: How do I get the feature from a single-feature OGR layer in Python without looping?

texts:
A while ago I asked why the 0 index in a OGR layer was not a feature and it turned out that OGR interprets the index as Feature ID and looks that up.

Now I filtered a layer to have only one feature in it. I would like to use that feature straight-away but cannot figure out how.

The setup:

from osgeo import ogr

driver = ogr.GetDriverByName("WFS")
wfs = driver.Open("WFS:https://geodienste.hamburg.de/HH_WFS_Statistische_Gebiete?REQUEST=GetCapabilities&SERVICE=WFS")  # This is a public WFS with open data
layer = wfs.GetLayerByName("Stadtteile")  # Some polygons, possibly multipolygons

# A random point inside one of them
point_in_hamburg = ogr.CreateGeometryFromWkt("POINT (566795 5935774)")

# Checking which features on layer intersect with our point
layer.SetSpatialFilter(point_in_hamburg)

# Due to the nature of Stadtteile (=districts) we
# will get only one intersecting feature:
print("One single feature: {}".format(len(layer) == 1))  # Prints "True"

Now I want to use that single feature. The only way I managed to is by using a loop over that length 1 layer:

for feature in layer:
    print(feature.GetFID())  # Prints "94"

I find that incredibly ugly and illogical. How can I directly access the feature instead without knowing anything about it?

layer[0] is not it. Strangely enough I can use layer[94] before the Spatial Filter and I get the feature (even though the WFS server returned a "Generic WFS service error") but if I try to use it afterwards, I get an IndexError from ogr. So even if the FID is 94 and is still 94 after filtering, layer[94] will fail at that point.

Is there some alternative?


texts: Go to XY tool not working when clicked in ArcMap?

texts:
I am using ArcGIS 10.3. When I click on 'Go to XY tool', located on toolbar, it doesn't open.

What is the possible solution to fix this problem, so I could open 'Go to XY tool' in ArcGIS 10.3?


texts:
Joining shapefiles is a relatively straightforward proess in QGIS using 'add vector join' in the layer properties menu. I can add a csv file as a vector layer, as well as a shapefile, and join the two using a common field.

I would like to do a similar join using pyQGIS, however I can't seem to find any documentation on how to go about this. Is this possible to do programmatically using pyQGIS with 'add vector join' or the older 'Join by Attributes'? Or, alternately, might there be any other Python modules that might be able to merge a csv with a shapefile before loading as a map layer with pyQGIS?


texts: I have a finished project but I have to make another with any information first (mostly) and I want to make a copy of the first draft and I just want to rename the mxd in my arcgis 10.1. Is there any way to do this without altering the project?


texts: SpatiaLite database error: no such table, but it is present in the SpatiaLite database

texts:
I am new to databases in QGIS, but cannot find any answer in all the documentation and online help in a weeks worth of desperate searching.  I think this might be a newbie error, but just cannot solve this issue.
I am trying to join two attribute tables from two point shapefiles [j_pipeline_inc (multipoint) to the j_pipeline_midpts (point)] using a SQL expression (image attached):
SELECT j_pipeline_midpts.*, j_pipeline_inc.*
FROM j_pipeline_midpts
INNER JOIN j_pipeline_midpts
on j_pipeline_midpts.LSD = j_pipeline_inc.LOCATION and j_pipeline_midpts.Licence_No=j_pipeline_inc.LICENCE_NO

The join will be based on two keys (LSD & Location, Licence_No & LICENCE_NO).  I have created a SpatiaLite database (incident_move.sqlite), which has these two point shapefiles in it, but when I execute my expression it states: "no such table: j_pipeline_inc".  I don't understand why it doesn't recognize one of the files.  When I type their names in the expression box it shows both shapefile names in the autocomplete.
I have tried saving one file as a SpatiaLite file type then importing the other to that database, but it still gives me the same error.
Also, they both have the same projection (EPSG:3400 - NAD83).


texts: Showing image (photo) stored as PostGIS base64 in ArcGIS Desktop

texts:
I have a PostGIS table which stores point data with images as attribute (in base64 data type).

I would like to show the pictures from the column "bild" of the attribute table which look like this in ArcMap:

What I want to do is to click a point object in the map view and then the attribute data of the point will be shown together with the corresponding picture. I have found this webpage and followed their tips as follows in ArcMap but it still does not work:

I have verified the text for the images by pasting it to a web browser and it shows the pictures without problem.

Does anyone have any idea what's wrong and how I should solve this in ArcMap.

I have ArcMap 10.3


texts: GeoServer publish raster mbtiles with png8

texts:
I had a raster mbtles, the raster was png8 in blob field. When I published the mbtiles, the level 0 (only one picture) was correct,

but other levels were wrong.

It seemed different color palette was loaded. How could I do?


texts: Installing Location Lab plugin in QGIS on Mac OSX?

texts:
I am having problems installing "Location Lab" plugin in QGIS 3.2 (Mac OSX High Sierra).

Looking for an option to run gravity models in QGIS, I've found the Location Lab plugin for QGIS, but the problem is that I could not been able to install it, every time I got this message:

  ModuleNotFoundError: No module named 'PyQt4'  Traceback (most recent
  call last):   File
  "/Applications/QGIS3.app/Contents/MacOS/../Resources/python/qgis/utils.py",
  line 337, in startPlugin
      plugins[packageName] = package.classFactory(iface)   File "/Users/davlira82/Library/Application
  Support/QGIS/QGIS3/profiles/default/python/plugins/location_lab/init.py",
  line 26, in classFactory
      from .location_lab import LocationLab   File "/Applications/QGIS3.app/Contents/MacOS/../Resources/python/qgis/utils.py",
  line 665, in _import
      mod = _builtin_import(name, globals, locals, fromlist, level)   File "/Users/davlira82/Library/Application
  Support/QGIS/QGIS3/profiles/default/python/plugins/location_lab/location_lab.py",
  line 24, in
      from PyQt4.QtCore import QSettings, QTranslator, qVersion, QCoreApplication   File
  "/Applications/QGIS3.app/Contents/MacOS/../Resources/python/qgis/utils.py",
  line 665, in _import
      mod = _builtin_import(name, globals, locals, fromlist, level) ModuleNotFoundError: No module named 'PyQt4'

I've already installed PyQt and upgraded QGIS from 3.0 to 3.2, but I receive the same message again and again. I'm actually an urban planner not a coder, I do not know what else to do.


texts: How to stick polygons next to each other in PostGIS or QGIS?

texts:
I've polygonized a territory, and then removed the polygons that were too narrow. Statistic values calculated for them made no sense, so I filtered polygons by shrinking by X metres and checking if anything was left.

Now, there are gaps between the polygons, and they are a minor problem.

I wonder if there's an algorythm to take those remaining, and if two of them are closer than X metres, stick them to each other. The margins after this operation should look about like this:

Is there any available solution, preferrably in PostGIS or QGIS?


texts: Styling plot items by category in QGIS print layout

texts:
I am making a sensitivity atlas which gives me gridded sensitivity values (A to E) for a given area. I am using the atlas functionality in QGIS 3.10 to create several maps showing these sensitivity attributts for different regions. The categories within a region is presented in several plots. I would like to link the layer styling to the styling of the plot layout.

In the below map you will see that I have a rule based styling where I give the values from A to E an increasing redness. I would like the statistics shown in the plot grids to have similar styling. Baing able to control the order of the objects in the plot (A to E) would also be a pre.

Is there a way to do this?


texts: Distance between multiple parallel lines?

texts:
I have a set of parallel lines where the default distance between them is 1.4m, I need to identify where that distance is greater or less than the default:

I have a sample file where it was done but I do not know how it was calculated nor how the points were generated:

The atribute table os points shp:

Searching in GIS, I found this topic with an answer that seems to be the solution, but I do not know how to reproduce for a set of lines. I quote:

  "You bet. There are a bunch of ways you can do it, but here's one potential method @Rashomon. For the set of points along each river bank, add new fields to the attribute table for Xcoord and Ycoord. Then open field calculator and use the $x and $y functions to populate the fields. You can then use anything you want to compute a simple average of the two sets of coordinates. Even Excel will work. Then just plot your new x and y coordinates that correspond to the center of the river. – lambertj Jan 25 at 16:44"

Create a center line between 2 polyline - QGIS

Solutions I thought:

Create a buffer in the 70cm lines and then identify where they overlap or do not touch, but I do not know how to calculate this.
Create centerlines between parallel lines, extract nodes every 1mt for example, and measure the distance to parallel lines.

There are many lines, thousands, it would need to be an automated solution as much as possible, it could be matlab, python, r or some other, or qgis or arcgis.

The sample files: https://drive.google.com/drive/folders/1XKOUECG3bY2claSqwQ-S9-C-yjYjE0r7?usp=sharing


texts:
https://openlayers.org/en/latest/examples/rotation.html is an example of a rotated map. If I take the extents of that, they do not reference the actual corners of the displayed map area but something prior(?) to rotation.

How do I get the actual extents or corner coordinates of a rotated map?

To see/reproduce the issue, you can copy and paste this into the browser console after loading above example and then zoom out.

var extents = map.getView().calculateExtent(map.getSize());
var poly_extents = ol.geom.Polygon.fromExtent(extents);

var p = new ol.layer.Vector({
  source: new ol.source.Vector({
      features: [new ol.Feature({
          geometry: poly_extents,
      })]
  }),
});

map.addLayer(p);


texts: How to cut contour lines at the edge of a data source

texts:
I have calculated some contour-lines using gdal-contour for two elevation model into a Postgis database. One is high res and the other lower res. Now I am seeing problems at the edge of my high resolution data-source where lines are creating a black border. See the attached image.

It seems as if the lines of the high resolution model continue over the edge, as if there has been a massive drop in altitude and, therefore, there are so many lines stacked on one of another.

My question now is: Can I somehow cut those lines? I thought I might be able to use the data part of the GeoTiff that I used for my high resolution contour-lines to generate a shape. Then shrink the shape by let's say 10m or so and use that shape to cut my contour line?!?
Is that somehow possible using PostGIS and gdal?


texts: QGIS 2.14.1 LAStools install error: "wrong value for parameter 'Msys folder'..."

texts:
I went to step 4 from How to install LAStools toolbox in QGIS guide and extracted LAStools to my C drive and just made the folder there.

I can't do step 8 because it seems to already be in "advanced interface". Then, step 9 mine only says "Options".

Step 10 is where I get the error message and it won't go past that step. I keep getting this error message. Can anyone help me get this installed?


texts: How to define Subroutines and functions?

texts:
I am working in VBA. My error is sub and function are not defined.  Here is my code:

Private Sub CreateMap(cboDistrict As String)
    Dim mapDocument As New mapDocument
    If mapDocument.IsPresent("K:\TASS\4_MAPPING_DATA_SUPPORT\Traffic_Mapping\District_Maps\2012\Abilene\Abilene_Base_Map.mxd") Then
        mapDocument.Open ("K:\TASS\4_MAPPING_DATA_SUPPORT\Traffic_Mapping\District_Maps\2012\Abilene\Abilene_Base_Map.mxd")
    End If
End Sub

I am trying to open an mxd when the Map button is pushed.


texts: Iterate Row Selection Using Current Selection ArcMap 10.2 ModelBuilder

texts:
I am trying to make a model using Iterate Row Selection. I have a buffer shapefile and a line shapefile. Ideally I would like for the model to loop through the buffer shapefile and select each row one by one.

When one of the buffer polygons is selection, I need to run a select by location using the current buffer selection to assign the FID of the buffer to a field called "FAF_FID" for the line shapefile

For example, if there are 3 line features with their centroid in Buffer FID 2, the three line features would be given a FAF_FID of 2.

I am looking for advice on how to accomplish this.


texts: Collector shuts down when attempting to use offline editing

texts: We have a service with a feature class and two related tables. All three have attachments. It's complicated, but it works fine when used online in Collector. When we download it for disconnected editing to a iPad 2 or iPhone 5s, it will open offline. However, when we try to add a new entry collector shuts down with no error message. This happens every time. I've tried republishing the service, uninstalling and reinstalling collector. Then I tested it on a fairly old and dumb Samsung tablet, and it worked fine. The tablet has an older version of collector; I'm not sure if that has anything to do with the problem. Has anyone else experienced this problem?


texts: ArcGIS 10 Network Analyst Service Area Problem/Peculiarity with Solving

texts:
I am having a problem when solving some of my service areas. I am loading the locations from a point feature class, and they are all located. The majority of them return a service area that is appropriate (7 minute breaks, double checked with solving for routes), but a handful of them return something strange. The problem is that there are some areas of the specific service areas that are wildly inaccurate, and there is no way that it is a 7 minute distance, which I also checked. I have attached a picture of the problem, and I'm assuming that it has something to do with my network, but I have no idea what. Any help/ideas would be greatly appreciated.

The network is a walking network that was manually created, and I am using ArcGIS 10.

Thanks,

Eric


texts: How can I use GDAL to batch define a projection?

texts: I have a fodler with 1500 ECW rasters, no spatial reference information defined.  Can I use GDAL to batch define a projection for all 1500 in the folder?


texts: Are there any free online GIS courses based around using QGIS?

texts:
Does anyone know of any MOOCs for learning GIS which use QGIS for any practical exercises?

I've been finding QGIS a fantastic tool for work, but would like to understand it and GIS principles on a deeper level.

However, many of the online courses or MOOCs I've found when looking for something like this are based around ArcGIS.


texts: Can't seem to move labels manually in QGIS, 1.80

texts: Trying to move labels manually in QGIS.  I am having the same problem when adding x,y to the attribute table.  I picked (Type Decimal Number (real)), it's the only one out of the 3 choices that has 'double'.  After enabling the Edit tool, the 'Move' and 'Rotation' tools remained greyed out. (1.80 Lisbo, using PC)  Help?


texts: I have been asked to check for pseudonodes - none are allowed for streams but they are for roads.  I can't figure out why so I assume I am not understanding what a pseudonode is.  I thought it was a node with no reason to be there. I am using ArcGIS 10.0.


texts:
I have two columns in my QGIS attribute table. The first column contains values that aren't contained the second column and vice versa. I would like to add a new field that contains both the values from the first and second column. I thought it might be as simple as "Value 1 + Value 2) but this just gives me Null results. All values are strings.

Value 1| Value 2 | New Column
-------------------------
 Bacon |         | Bacon
 Eggs  |         | Eggs
       | Cheese  | Cheese
       | Ham     | Ham


texts: Custom SRS for openlayers 3 (ol3) without an EPSG code

texts:
How can openlayers 3 be told to use and work with a custom SRS-definition which has no EPSG-code?

May i don't understand it but it seems like ol3 needs an EPSG- or ESRI-code
to work.

I found some examples about how to use WKT- and proj4- definitions unsing Prof4js.defs, but it didn't work for me:

Proj4: https://stackoverflow.com/questions/29027295/how-to-use-openlayers-3-with-proj4js

OpenLayers 3 transform coordinates to custom projection

https://stackoverflow.com/questions/29302770/parsing-projection-wkt-in-openlayers-3

I had several errors caused by the line with the definition (proj4.defs) saying the projection is undefined. I used a WKT and a proj4 declaration in there - both were creating the same Error.

I have a WKT- and a simple proj4-definition for the needed reference body, but i don't understand how to combine it to make it work in ol3.

The simplified js I use:

var baseLayers = new ol.layer.Group({
    title: 'Base Layers',
        openInLayerSwitcher: true,
        layers:
        [   new ol.layer.Tile({
                title: "testLayer1",
                baseLayer: true,
                source: new ol.source.OSM({
                    url: 'tiled_layer1/{z}/{x}/{-y}.png'
                })
            }),
            new ol.layer.Tile({
                title: "mos_final_14_02_13.map",
                baseLayer: true,
                visible: false,
                source: new ol.source.OSM({
                    url: 'tiled_layer2/{z}/{x}/{-y}.png'
                })
            })
        ]
    });
var overlays = new ol.layer.Group({
    title: 'moasaic single frames',
        openInLayerSwitcher: true,
        layers:
        [   new ol.layer.Tile({
                title: "overlay1",
                baseLayer: false,
                source: new ol.source.OSM({
                    url: 'tiled_overlay/{z}/{x}/{-y}.png'
                })
            }),
            new ol.layer.Tile({
                title: "overlay2,
                baseLayer: false,
                visible: false,
                source: new ol.source.OSM({
                    url: 'tiled_overlay2/{z}/{x}/{-y}.png'
                })
            })
        ]
    });
    var map = new ol.Map({
        target: 'map',
        view: new ol.View({
            center: ol.proj.transform([0, 0], 'EPSG:4326', 'EPSG:3857'),
            zoom: 3,
            minZoom: 0,
            maxZoom: 9
        }),
        controls: ol.control.defaults().extend([
            new ol.control.LayerSwitcher()
        ]),
        layers: [baseLayers, overlays]
    });

OpenLayers 3 Examples showing this: http://openlayers.org/en/v3.3.0/examples/wms-image-custom-proj.html

They use a JS-file telling the Definition for the Projection - but i didn't find the place where they use it - its just standing in thier souce-code.


texts:
I'm trying to create a view of two heatmaps in the same time with QGIS.  I tried following the suggestion on this thread:
Can I combine two heat-maps in QGIS?
But it didn't work.  The virtual layer only showed one of the heatmaps and not the other, as in this screenshot.  Is there another workaround for this?


texts: is it possible to estimate the complex coherence between two RADARSAT-2 SLC images with SNAP or SarScape? If I try to estimate the Coherence I get only float data but I need it in a complex form.


texts: Can I specify a minimum scale for atlas features?

texts: Is there a way of specifying a minimum scale for all the features of the atlas? (without using the fixed scale option) and for all scales above that "predefined scale (best fit)"? I want to have a minimum scale for the features in the atlas which are too zoomed in and best fit for any scales over that.


texts:
I am new to python and breaking my head on how to use gdal_merge in a (spyder) python script. I use Windows 10, Python 3.6 and have the gdal tool installed from osgeo4w. I realize many other posts describe this problem but none could help me resolving this issue.

When i call the gdal module from cmd, it works like a charm:

 python "C:\OSGeo4W64\bin\gdal_merge.py" -o merged.tif input1.tif input2.tif

However, I cannot get it working properly in a python script (spyder).

A first approach produces an output but not with the right name: it produces an 'out.tif' file and not the 'merged.tif' file as I request:

import sys
sys.path.append('C:\\OSGeo4W64\\bin')
import gdal_merge as gm
gm.main(['-o', 'merged.tif', 'N00W078.tif', 'N00W079.tif'])

A second approach simply produces no output:

import subprocess
merge_command = ["python", "C:\OSGeo4W64\bin\gdal_merge.py", "-o", "merged.tif", "input1.tif", "input2.tif"]
subprocess.call(merge_command,shell=True)

Any thoughts on how to resolve this issue?


texts: Extracting poly.counts through loops using R

texts:
I want to extract the count of points over Africa at the first administrative district level at monthly intervals over 26 years. My point data consist of conflict events. I am struggling to figure out how to do this without spending a lifetime copy-pasting code.
I need to capture three things: the country, administrative district name, and the count of events.
My intuition is to use poly.counts() from GISTools in some sort of a loop, such as:
for(y in 1990:2016){
    for(m in 1:12){
       poly.counts(subset(data, year == y & month == m))
    }
}

I'm not sure how to capture the output of each iteration into a single data frame. Ideally, this will produce a single data frame that shows the year, month, country, administrative district, and count. Below, I present a short "toy" version of what I want.
year     month     country     admin1     count
----     -----     --------    -------    ------
2003     1         Nigeria     Lagos      67
2003     1         Nigeria     Oyo        12
2003     1         Nigeria     Kaduna     54
2003     2         Nigeria     Lagos      12
2003     2         Nigeria     Oyo        24
2003     2         Nigeria     Kaduna     56
2003     3         Nigeria     Lagos      15
2003     3         Nigeria     Oyo        22
2003     3         Nigeria     Kaduna     14


texts: ArcGIS Raster Calculator error 000539 ImportError: No module named numpy?

texts:
I am running ArcGIS 10.2.1 on Win7 64bit.
I have finished doing an NDVI, and after defining a vegetation presence threshold at 0.05, opened the Raster Calculator (via Map Algebra in the Arc Toolbox) to generate a new raster with a 0 value for no vegetation and a 1 value for vegetation. I'm getting the following error:

Executing: RasterCalculator ""NDVI_ik_06.tif" > 0.05"

ERROR 000539: Error running expression: rcexec()
Traceback (most recent call last):
  File "<expression>", line 1, in <module>
  File "<string>", line 2, in rcexec
  File "c:\program files (x86)\arcgis\desktop10.2\arcpy\arcpy\__init__.py", line 24, in <module>
    from arcpy.toolbox import *
  File "c:\program files (x86)\arcgis\desktop10.2\arcpy\arcpy\toolbox.py", line 356, in <module>
    from management import Graph, GraphTemplate
  File "c:\program files (x86)\arcgis\desktop10.2\arcpy\arcpy\management.py", line 22, in <module>
    import _management
  File "c:\program files (x86)\arcgis\desktop10.2\arcpy\arcpy\_management.py", line 14, in <module>
    import _graph
  File "c:\program files (x86)\arcgis\desktop10.2\arcpy\arcpy\_graph.py", line 27, in <module>
    import numpy
ImportError: No module named numpy

Failed to execute (RasterCalculator).
Failed at Tue Apr 07 11:24:48 2015 (Elapsed Time: 0.09 seconds)

Any ideas?


texts: Getting field value in Print composer label

texts:
I am using QGIS print composer. I have so many images to export with a specific label. The text of the label is existing in a field of a layer.

I have tried to use this expression attribute( get_feature( 'Metro_Ring','S','1'),'Route') but the result is NULL.

where:

Metro_Ring --> name or ID of the table

'S' --> the field for selection.

1 --> the value that selects the row of interest.

'Route' --> the name of the field of interest.

I wonder if there is any error or a different method to get the value of the field 'Route' in the label.


texts: Invoking Context-sensitive Help for Python AddIn Button?

texts:
In an experimental Python AddIn that I've created using ArcPy from ArcGIS 10.3.1 for Desktop I have put some test help into my config.xml for a button called Button4:
<Button caption="Button4" category="Experimental Addin" class="ButtonClass6" id="Experiment_addin.button4" image="" message="Button4 Message goes here" tip="Button Tip goes here">
    <Help heading="Button4 Help Heading goes here" Content="Button4 Help Content goes here" />
</Button>

When I make and install my Python AddIn, the message and tip show up by mousing over the button as expected.
However, I have not figured out how to invoke the Help Heading and Content mentioned in Creating a Python add-in tool:

The final section of the tool is the help heading and content. These
properties allow you to supply information that appears when a user
invokes context-sensitive help. These are pop-up topics that remain
on-screen until the user clicks somewhere else.

Prior to 10.1 I know that this could be done using the now deprecated What's This button, but how is context-sensitive help invoked since 10.1?
I thought it might have been by mousing over the button and hitting F1, but that does nothing.  If I mouse over anything but a toolbar and hit F1 it brings up the ArcGIS 10.3.1 Help.


texts:
I have a Postgres/postgis table of polygons. When added to QGIS the total feature Count is 64671 in layer tree and attribute table. But the attribute table also says Filtered: 43573. I have not applied any filter.

What is total Count (what is 64671 and what is 43573)? What filter is applied?

Does not seem like it is the total single part Count either:


texts: Is there any specific case study to follow when developing a postcode system?

texts:
I'm developing a postcode system to a country that doesnt exist one (in fact we have one to reflect a main post office with a really large area to cover -less granularity- hence no one really uses it).
I would like to know if there is any project case I can follow something similar to this?

I'm not a developer or programmer but has a team to carry out the work. We are still in the very early stage of the process.


texts: KML and Spot and Swipe widget in flexviewer 2.2

texts:
I found the KML widget for Flex viewer but I was not able to make it work properly. Can anybody have any link for that.
I came across SPOT and SWIPE tool code, but the links to download the code redirected to something else. I thing the links are dead. So any body have that code or widget for SPOT and SWIPE tool.

Thanks in advance


texts: Aligning coordinate systems: gpx, LiDAR, shapefile

texts:
I am working on a project based on a lidar data set (for use in LidarViewer in Linux, a basic x,y,z viewer) that was downloaded from #1 below at Opentopo.org and looking to add onto this data set with a shapefile based layer and a gpx based layer.

In ArcMap 10.0 I am converting the shapefile (road perimeter around lake) into points, then Defining the Projection in Arc Toolbox's Projections and Transformation.  Since Opentopo stated UTM Zone 10, I set that to the roads shapefile, exported, completed my .shp to .las(lidar file) conversion and sent it to LidarViewer.  In LidarViewer my LiDAR data (only a section of the lake's perimeter) is floating the in the middle of the lake (with perspecitve to the road points).

So, I switched the Tahoe to the original NAD27 in ArcMap, went through the process again with no change in position. Do I have to have a shapefile based on NAD83 for this to work, or is there another way?

Lidar viewer imports .las or .txt/.csv files and only reads a x,y,z no coordinate system whatsoever.

The other related question is with the .gpx data that I am downloading from the #3 below, it is a .gpx file of a bike path around the perimeter of the lake. I am using the following link to convert .gpx into a text file (to process and import into Lidar Viewer) and the site provides the ability to output long/lat plus UTM coordinates (the output included both zone 10 and 11 South ) This .gpx file also does not project correctly. I was thinking that I should possible try another converter?

gpsvisualizer .com  convert_format=gpx

Any input regarding the above two questions are appreciated, I can also further clarify any issues that may not seem clear.

Thanks!

1) opentopography-org -> Data -> South Lake Tahoe, Nevada

Horizontal Coordinates: UTM Zone 10 N NAD83 Meters [EPSG: 26910]  -  Vertical Coordinates: NAVD88

2) tahoe usgs gov/DLG   UTM NAD27 zone 10 projection

3) bicycling. trimbleoutdoors. com/ViewTrip.aspx?tripID=153934


texts:
I've just unclicked the export map options (advanced options?) and can't get them back. I'd assumed you'd be given the choice of displaying them again, but they are nowhere to be seen. I now can't choose what file type or resolution to use, which is a huge pain.

http://forums.arcgis.com/threads/42377-Export-map-no-options-for-details-filename..

Of the options above, I've tried deleting normal.mxt, with no luck. There's  no 'ExportDlg' folder, so I don't have that option.

I'm going crazy with this. I've Arc 10.1 running on Windows 8. Any workaround / solution would be hugely appreciated. I wont have the time to uninstall.

Here's a screenshot:


texts: SUDS (SOAP) from the Python console in QGIS

texts:
I´m trying to install a SOAP API called SUDS, on a Windows 7 OS, to be able to run SOAP from the QGIS Python console and pyQGIS/Python scripts.

The only Python I have installed is the one comming with QGIS 2.01 standalone. How can I install an external Python library like SUDS, to be used within QGIS Python scripts and console?


texts: Seeking data source for worldwide urban/statistical areas as polygons?

texts:
I am looking for a shapefile (or any format really) that has worldwide coverage of major urban areas / statistical areas as polygons with names.

Similar to what the US census Bureau' Core Based Statistical Areas but with worldwide coverage.

See https://www.census.gov/geo/maps-data/data/cbf/cbf_msa.html

I need polygons, because I am trying to segment another layer of coordinates with it.

I know similar questions have been asked but it seems that most have the limitation of needing free data.  I'm happy to pay for good data if its available.


texts: SQLite db with FDO geometries: how to assign ID of point in layer A to closest point in Layer B with SpatiaLite?

texts:
How can I assign ID of point in layer A to closest point in Layer B with SpatiaLite (find nearest neighbors for each point in A)?

For each point stored in table fdo_fc1 I would like to select one, nearest point from table fdo_fc2.

Table fdo_fc1 has field id1 and contains 17 records. Table fdo_fc2 has field id2 and contains 33 record. (test data; real data has about 100,000 records)

I work with a  sqlite database created with AutoCAD Map 3D. To query my database I use spatialite_gui.exe

My result table should look like this:

I have tried a lot but I can not find a solution.

I simplified code I found (improve performance using  a maximum distance is not interesting for me now).

E.g. the accepted answer  to How to find nearest neighbors between two tables with point locations in SpatiaLite? do not work for me:

Example 1: I tried:

SELECT t1.id1 As id1, t2.id2 As id2, min(Distance(t1.GEOMETRY,t2.GEOMETRY)) AS DIST FROM fdo_fc1 As t1, fdo_fc2 As t2

But result only contains one record form table fdo_fc1:

Example 2: I tried:

SELECT t1.id1 As id1, t2.id2 As id2, Distance(t1.GEOMETRY,t2.GEOMETRY) AS DIST FROM fdo_fc1 As t1, fdo_fc2 As t2

The result contains the distance only for first record of fdo_fc1:

Here is map of my test data:

EDIT:

I tried the solution of @user30184. I have noticed the following:

I created the sqlite database with AutoCAD Map 3D. The names of the layers/tables are fc1 and fc2.
When I open the database with spatialite_gui.exe this message opens:

So now I have the new tables fdo_fc1 and fdo_fc2:

If I use fdo_fc1 and fdo_fc2 I only get one record in result. If I use fc1 und fc2 I get the correct number of records in result but the distance is not calculated:

This query returns still one record:

SELECT t1.id1 AS id1, t2.id2 AS id2, Min(Distance(t1.GEOMETRY,t2.GEOMETRY)) AS DIST FROM fdo_fc1 AS t1, fdo_fc2 AS t2 GROUP BY t1.id1 ORDER BY id1;

This query returns no distances:

SELECT t1.id1 AS id1, t2.id2 AS id2, Min(Distance(t1.GEOMETRY,t2.GEOMETRY)) AS DIST FROM fc1 AS t1, fc2 AS t2 GROUP BY t1.id1 ORDER BY id1;


texts:
I have a Python script that I have exported from a QGIS model and I want to specify the layer it creates as output.

If defaultValue=None it outputs 'Reportpmbc': 'Centroids_9c779c9c_4bcc_471b_971e_fbea3bb27c7d', in the log window and a 'Centroids' layer is created in the project.

I want it to create a temporary layer 'Report-PMBC' as below but it won't run.  (“.” files are not supported as outputs for this algorithm)

What am I missing in the statement below? I think I have to create a temporary layer in memory and write to it but I am not sure how to do that.

 self.addParameter(QgsProcessingParameterFeatureSink('Reportpmbc', 'Report-PMBC', type=QgsProcessing.TypeVectorPoint, createByDefault=True, defaultValue='Report-PMBC'))


texts: how to solve "could not access file '$libdir/librouting'" problem in ubuntu 12.04

texts:
I am trying to install pgrouting in my ubuntu 12.04 machine using the instructions of this link. However, when I try to execute the following command based on step 3, it generates some error.

command

      psql -U postgres -f /media/resource/pgrouting/pgrouting/core/sql/routing_core.sql routing

Output:

      psql:routing_core.sql:21: ERROR:  type "path_result" already exists

      psql:routing_core.sql:22: ERROR:  type "vertex_result" already exists

      psql:routing_core.sql:32: ERROR:  could not access file
      "$libdir/librouting": No such file or directory

      psql:routing_core.sql:43: ERROR:  could not access file "$libdir/librouting": No such file or directory

      psql:routing_core.sql:53: ERROR:  could not access file "$libdir/librouting": No such file or directory

      CREATE FUNCTION

How can I solve this problem?


texts:
I have imported a dataset from a url http://api.gdeltproject.org/api/v1/gkg_geojson?QUERY=TERROR%22&TIMESPAN=1440&OUTPUTFIELDS=url,name,sharingimage and I want to sync it now using API.

I followed the instructions here https://carto.com/docs/carto-engine/import-api/sync-tables and typed in

curl -v -H "Content-Type: application/json" -d '{"http://api.gdeltproject.org/api/v1/gkg_geojson?QUERY=TERROR%22&TIMESPAN=1440&OUTPUTFIELDS=url,name,sharingimage":"https://public.url.to.file/sample_file", "interval":"3600"}' "https://sengyapglobalmomentum.carto.com/api/v1/synchronizations/?api_key=API_KEY"

However, this error shows up

Did I carried out any wrong steps?


texts: Realising vector file mapping using QGIS Modeler?

texts:
I pretty often get lots of vector files with different column namings. I need to produce a defined vector data file with certain columns. Is it possible to build a model that allows to make a mapping? I hope my example is explaining my question:

I got a file with these two columns "City" and "Road". I need to get an output file with my two standard column names "Town" and "Street" these two new columns should contain the information of the original columns. Of course it would also work just to rename the input columns.

Input columns: "City", "Road"
Output columns: "Town", "Street"


texts: Finding points of intersection in polygon?

texts:
For two given polygons generated as shown below, how do I find the points that are intersected on the rectangle by the circle?
[The polygon was created using the Alphashape function]

(Similar to Returning percentage of area of polygon intersecting another polygon using shapely which is answered but I have not used Shapely to generate polygon for the distribution)

Code implemented:

import sys
from descartes import PolygonPatch
import matplotlib.pyplot as plt
import alphashape
import numpy as np
x = np.random.uniform(41.03641664650488, 59.06406419360736, 1000)
y = np.random.uniform(100.96028140425065, 118.84533708184047, 1000)
test = np.vstack((x,y)).T
alpha_shape = alphashape.alphashape(test, 0.)
fig = plt.figure(1, figsize=(12, 5))
ax = fig.add_subplot(111, aspect='equal')
ax.scatter(*zip(*test))
ax.add_patch(PolygonPatch(alpha_shape, alpha=0.2))
circle = plt.Circle((50,100), 4, color='r', fill=True)
plt.gca().add_patch(circle)
plt.autoscale()
plt.show()


texts:
I'm trying to extrapolate from a line segment to find a point on the line but a 3rd of the way 'back', i.e. trying to find point new, given points A and B below:

Given a line, I can interpolate it to get a position at any particular percentage along it:

=# select st_line_interpolate_point(
   st_makeline('0101000020E6100000300DC347C49418C03EE8D9ACFAA44A40',
               '0101000020E6100000FB743C66A03218C0CDCCCCCCCC7C4A40'),
   0.333);
0101000020E6100000ED45B41D537718C069C6A2E9EC984A40

I tried entering a negative number to find a point along the line in the opposite direction, but that fails as the interpolation argument has to be in the range [0, 1]

I thought about first scaling the line, but that doesn't use the centre of the line as the origin, so is useless for my purposes.


texts: Joining data from two shapefiles in QGIS

texts:
I think there must be a straight forward answer to this but I haven't managed yet:

I have two shapefiles for the same watershed they are both divided into different areas with different erosion factors. I am trying to get a new shapefile with, say, the sum of both factors (a column in each file).
The problem is that a polygon from one file isn't necessarily covering the same area as a polygon from the other file.

What would be the best approach to this task?


texts: gdal2tiles.py generating black tiles from .img file

texts:
I have a source HFA Erdas Imagine (.img) file which I'm trying to convert into tiles using gdal2tiles.py. I'm running:

gdal2tiles.py -w google -g GOOGLEKEY input_file.img

The HTML is generated correctly, and the tiles are correctly positioned. However, I'm getting black tiles instead of the actual content of the .img file:

I've verified (using QGis and MapTiler) that the .img contains the correct data, i.e. the source file has raster data.

Any ideas? Should I transform the .img to something else? Am I missing an option for the command?


texts: Exporting text object coordinates in DWG file to excel

texts:
How can I export the coordinates of each of the text objects in the layers
MovePLOTNUM and PLOTSNUM to excel sheet (columns: layer, text value of the object, x, y). Is this possible with AutoCAD or I need another software? here is the file. I tried to search the web to figure out how is this possible and could not find a way.


texts:
Trying to check the type of some Elements if they are (ITextElement , MarkerElement), when i try .getType() function it returns always "System._COM_object".

How i could get the type of an object if i don't know its type before ? because i could only detect the types that i pretend using (obj as Interface) != null.

but for the types that i don't pretend i get always null from (obj as Interface) .

Any help please.


texts: Extracting Colour (Brush) information from MapInfo Workspace

texts:
I have some MapInfo Tables with an associated Workspace. I do not have access to MapInfo and wish to use the TAB files in QGIS (I also have the equivalent SHP files). The TAB/SHP files have polygons with associated attribute data. Colouring of the polygons was done in the Workspace (.WOR) file with the relevant part reproduced here:

Set Map
  CoordSys Earth Projection 3, 104, "m", 134, 0, -18, -36, 0, 0
  Center (-1052836.855,-3057718.85)
  Zoom 6700.860332 Units "km"
  Preserve Zoom Display Zoom
  XY Units "m" Distance Units "km" Area Units "sq m"
shade 2 with AGE_CLASS values
  "" Brush (2,16711680,16777215) Pen (1,1,0) ,
  "Archaean" Brush (2,9800292,16777215) Pen (1,1,0) ,
  "Archaean to Proterozoic" Brush (2,13416577,16777215) Pen (1,1,0) ,
  "Cainozoic" Brush (2,16777168,16777215) Pen (1,1,0) ,
  "Mesozoic" Brush (2,5680257,16777215) Pen (1,1,0) ,
  "Mesozoic to Cainozoic" Brush (2,13696976,16777215) Pen (1,1,0) ,
  "Palaeozoic" Brush (2,10535167,16777215) Pen (1,1,0) ,
  "Palaeozoic to Cainozoic" Brush (2,13697023,16777215) Pen (1,1,0) ,
  "Palaeozoic to Mesozoic" Brush (2,13693695,16777215) Pen (1,1,0) ,
  "Proterozoic" Brush (2,16756912,16777215) Pen (1,1,0) ,
  "Proterozoic to Mesozoic" Brush (2,16752895,16777215) Pen (1,1,0) ,
  "Proterozoic to Palaeozoic" Brush (2,16765183,16777215) Pen (1,1,0)
  default Brush (1,0,16777215) Pen (1,2,0)   # color 1 #

Is it possible to convert the Brush values to RGB?


texts: IMapControlEvents2.OnKeyDown not hit on Alt + key

texts:
I want to catch the Alt key plus an addionally hit key in the MapControl's OnKeyDown. I did already set MapControl's KeyIntercept to esriKeyIntercept.esriKeyInterceptAlt but the event handler doesn't get raised on Alt or Alt + key. It does get hit on other keys.
Any suggestions?

_mapControl.KeyIntercept = (int)(esriKeyIntercept.esriKeyInterceptAlt);

Handler:

void _mapControl_OnKeyDown(object sender, IMapControlEvents2_OnKeyDownEvent e)
{
  Trace.WriteLine("hit");
}


texts:
I have a GPS log where speed, heading are obtained every 100ms. I would like to calculate yaw rate with this. Can I get the method, algorithm or formula to calculate the same?

I understand using Discretized Second Order Low-Pass Filter with Differentiator formula must be used but unsure of the formula.