texts: Fastest way to numerate intersecting polygons in ArcGIS

texts:
I need to numerate a thousands of polygons in groups, so each group will get unique number. Every intersecting polygons is unique group (see picture).

So polygons have field "groups" with values:

polygons 1,2,3 is group 1;

polygons 4,5 is group 2;

polygons 6,7,8,9 is group 3 and so on...

My current approach is very slow it takes about 5 minutes to enumerate layer with about 10 000 polygons (it is 4 100 groups).

Is there any faster way to do this?

Here is what I do:

Dissolve all polygons - that is my way to find groups
(arcpy.Dissolve_management)
Give unique numbers to every polygon in dissolved layer (arcpy.CalculateField_management)
Iterate every feature in dissolved layer (arcpy.SelectLayerByAttribute_management)
Select polygons that intersect with original (not dissolved) layer.
Give them group number (arcpy.CalculateField_management)

Here is a part of the code:

#1
arcpy.Dissolve_management(fc, dissolve_shp, "", "", "SINGLE_PART")

#2
arcpy.CalculateField_management(dissolve_shp, "ID", '!FID!+1', "PYTHON")

#3
number_of_rows = arcpy.GetCount_management(dissolve_shp)
number_of_rows = int(number_of_rows.getOutput(0))
lyr_dissolve = r"in_memory\Temporary_layer_DISSOLVE_DEL"
lyr_merge = r"in_memory\Temporary_layer_MERGE_DEL"
arcpy.MakeFeatureLayer_management(dissolve_shp, lyr_dissolve)
arcpy.MakeFeatureLayer_management(fc, lyr_merge)

for i in range(number_of_rows):
    i2=i+1
    if float(i2)//100 == float(i2)/100:
        print i2, "of", number_of_rows

    arcpy.SelectLayerByAttribute_management (lyr_dissolve, "NEW_SELECTION", '"ID" = '+str(i2))

    #4
    arcpy.SelectLayerByLocation_management (lyr_merge, "INTERSECT", lyr_dissolve)

    #5
    arcpy.CalculateField_management(lyr_merge, "TERM_P_NUM", i2, "PYTHON")


texts: Cell size in Grid (Nearest Neighbour) in QGIS

texts:
I am using Grid (Nearest Neighbour) in QGIS to obtain raster layer from point shapefile.
However, on execution of the tool and obtaining the output, the default resolution of cell is 62 m. However, I want to change the resolution to 1 m x 1 m. How do I do that. ? I am unable to find any option for that.


texts: ArcGIS Raster Symbology - Break contour at 0

texts: In ArcGIS 10, how can I break raster symbology (stretched) at a particular value (i.e. 0).  For example, I have a DEM representing both topography and bathymetry and I want a break at the 0 m contour line (shoreline) so the green colors represent positive elevation (land) and the blue colors represent values less than 0 (water depths).  Similar to what Global Mapper offers.


texts:
I am not sure if this question makes sense; I have a 3D model segment of land. Because it is a 3D model, it does not have a CRS(??). I process it in CloudCompare and eventually export a raster of some parameter. I then import it into QGIS along with a pre-made shapefile of the same area (CRS: WGS84). However, the two do not match up because the raster does not have the same CRS or extent. I have tried using R to project the raster to the shapefiles CRS, but I don't have a starting CRS, so R wouldn't know how to convert it. So how does one put a CRS into a 3D model, preferably in R.

Files that I am using (scene_mesh_decimated_textured)

I was given another file (image 2) that already had a CRS attached to it (WGS84). I did my processes and when imported the raster into QGIS, it worked. I don't understand how the 3D model could have a CRS with it? Is it because it is a .las file type?


texts: How to add an attribution on a GeoJSON Layer from Leaflet?

texts:
I need to use a GeoJSON Layer on my Leaflet map. Here is a sample of my code:

function onEachFeature(feature, layer) {
    if (feature.properties && feature.properties.popupContent) {
        layer.bindPopup(feature.properties.popupContent);
    }
}

myGeoJsonLayer = L.geoJson(data, {
    pointToLayer: function (feature, latlng) {
        return L.circleMarker(latlng, geojsonMarkerOptions);
    },
    onEachFeature: onEachFeature
});
myGeoJsonLayer.addTo(map);
TOC.addOverlay(myGeoJsonLayer, "My GeoJSON Layer");

All is working.

Now I'd like to add an attribution on my layer but how?


texts:
If I am creating data from satellite imagery, is this classed as primary or secondary data?

If I am creating slope data from a DTM, is this then classed as primary or secondary data?


texts:
A lot of our datasets have a nodata value that is, or close to, -3.4028234663852886e+38.

Can anyone think of an easy way to say this out loud when we're discussing things?

Even "minus three times ten to the thirty-eight" is a bit of a mouthful. Something like "min float" would be shorter (I think that's what this value is supposed to represent), but is not strictly correct.

The best we've come up with so far is "minus three whatever", which is short, but not exactly scientific.


texts: Kriging with a non-uniform (not space filling with concentrated data in certain regions) sampling design?

texts:
I am using Kriging for optimization where local searches are often carried out. However, because of these local searches my sample design becomes poor and my Kriging prediction as well as variance estimates are poor. I read that this is because of ill conditioning of correlation matrix.
Choosing a small subset of data is not preferred for me because then you are just losing information.
Question: Are there any other alternatives (including subset selection given that I am able to utilize most of my data)? I am using a nugget term and that doesn't seem to help much.


texts: I've recently been upgraded to MapInfo 12.5 I'm trying to trace a MapInfo .MBX tool called labeler - I've used it in the past to transfer the geographic position of current MapInfo labels to a separate text layer. I've looked on the PBinsight community downloads page, and it doesn't seem to be there to download.  It also doesn't appear in the MapInfo Tools menu when I use MapInfo and it is not in the Program files directory on my pc. Any ideas where I can locate this file to load it into my pc?


texts: Problem accessing English Environment Agency API with QGIS

texts:
I am a GIS novice attempting to use QGIS (version 3.12.3-București) to link to a UK Environment Agency flood zone dataset using the various data links provided here:

https://data.gov.uk/dataset/cf494c44-05cd-4060-a029-35937970c9c6/flood-map-for-planning-rivers-and-sea-flood-zone-2

Using the base WFS/WMS URLs provided just returns an error, but following the advice at 'https://support.environment.data.gov.uk/hc/en-gb/articles/360001340657-How-do-i-add-WFS-layers-in-ArcGIS-and-QGIS-' and adding parameters, I do manage to import some partial data using this URL:

https://environment.data.gov.uk/spatialdata/flood-map-for-planning-rivers-and-sea-flood-zone-2/wfs?SERVICE=WFS&VERSION=2.0.0&REQUEST=GetFeature&TYPENAME=Flood_Map_for_Planning_Rivers_and_Sea_Flood_Zone_2&BBOX=450660,332709,461502,347119,EPSG:27700&srsName=EPSG:27700

The layer data imported isn’t correct (it doesn’t match any watercourses and only covers a small area near Nottingham - see screenshot), but it at least shows that the API link works!

Does anyone have any experience of accessing this specific data using and has made it work correctly with QGIS, or can at least suggest what else I might try?


texts: Get all features in extent not working correctly in PageLayoutControl

texts:
I want to get all features in the current extent of my PageLayoutControl's ActiveView. I got it working for a MapControl as described here.

With my PageLayoutControl it seems that I don't get the correct extent. When I fill the found features it looks like this:

Here is how I get the features:

Dictionary<String, IFeatureCursor> featureCursors = new Dictionary<String, IFeatureCursor>();
IMap map = this._pageLayoutControl.ActiveView.FocusMap;

// Get ActiveView from MapFrame
IActiveView tmpActiveView = null;
IGraphicsContainer graphicsContainer = (IGraphicsContainer)_pageLayoutControl.PageLayout;
graphicsContainer.Reset();
IElement element = graphicsContainer.Next();
while (element != null)
{
  if (element is IMapFrame)
  {
    IMapFrame mapFrame = (IMapFrame)element;
    tmpActiveView = (IActiveView)mapFrame.Map;
    break;
  }
  element = graphicsContainer.Next();
}
if (tmpActiveView == null)
  return featureCursors;

IEnvelope envelope = this._pageLayoutControl.ActiveView.ScreenDisplay.DisplayTransformation.VisibleBounds; // FittedBounds;
// PageLayout has no SpatialReference. So use the View of the MapFrame
IPoint pointMax = tmpActiveView.ScreenDisplay.DisplayTransformation.ToMapPoint((System.Int32)envelope.XMax, (System.Int32)envelope.YMax);
IPoint pointMin = tmpActiveView.ScreenDisplay.DisplayTransformation.ToMapPoint((System.Int32)envelope.XMin, (System.Int32)envelope.YMin);
envelope.XMax = pointMax.X;
envelope.YMax = pointMax.Y;
envelope.XMin = pointMin.X;
envelope.YMin = pointMin.Y;

// Get Features in all layers
for (int iFeatLayer = 0; iFeatLayer < _pageLayoutControl.ActiveView.FocusMap.LayerCount; iFeatLayer++)
{
  ILayer thisLayer = _pageLayoutControl.ActiveView.FocusMap.Layer[iFeatLayer];
  IGeoFeatureLayer featLayer = thisLayer as IGeoFeatureLayer;
  if (featLayer != null)
  {
    IFeatureClass featureClass = featLayer.FeatureClass;
    System.String shapeFieldName = featureClass.ShapeFieldName;

    // Create a new spatial filter and use the new envelope as the geometry
    ISpatialFilter spatialFilter = new SpatialFilterClass();
    spatialFilter.Geometry = envelope;
    spatialFilter.SpatialRel = esriSpatialRelEnum.esriSpatialRelEnvelopeIntersects;
    spatialFilter.set_OutputSpatialReference(shapeFieldName, map.SpatialReference);
    spatialFilter.GeometryField = shapeFieldName;

    // Do the search
    IFeatureCursor featureCursor = featureClass.Search(spatialFilter, false);
    featureCursors.Add(thisLayer.Name, featureCursor);
  }
}
return featureCursors;

What did I miss? :)


texts:
I'm testing qgis triangulation. So far I managed to extract a single polyogon, simplify it, and convert to points.
The problem is that the triangulation produces something unexpected.

This is the points

And this is after the triangulator.

My expectation was that the process could reproduce the shape of the country.
My overall objective is to export the shapefile to a 3d program, since I don't see export options to 3d formats I must resort to points and triangles.
Can anyone suggest a different approach?


texts: Setting name for composer template using PyQGIS?

texts:
This scripts works fine for me, but every time its called, a new composer-template is created (Composer 1, Composer2 etc).
If the user saves the project, also the composers are saved and the list grows up...

So I want to remove the new item OR name it, so only one of them is in the list:

from qgis.PyQt.QtXml import QDomDocument
from qgis import gui

scale = 10000
templateFile='c:/1temp/mylayout.qpt'
activetool = 15# 'MoveItemContent'

composerview = iface.createNewComposer()
#load template
f= file(templateFile, 'r')
templateContent = f.read()
f.close()
doc=QDomDocument()
doc.setContent(templateContent)
composerview.composition().loadFromTemplate(doc)
#show composer
composerview.show()
#setup map
map = composerview.composition().getComposerMapById(0)#QgsComposerMap
#extent
map.zoomToExtent(iface.mapCanvas().extent())
#scale
if scale >0:
    map.setNewScale(int(scale))
#set activetool
if activetool is not None and activetool<>'':
    composerview.setCurrentTool(gui.QgsComposerView.Tool(activetool))
#finally
composerview.composition().refreshItems()


texts:
I have been struggling with something for a while that I figured would be simple. I have a WMS layer with many (millions) polygon features. This is a PostGIS layer (WMS). Each polygon has an ID in an attribute. I have a web page set up where a user can enter an ID and click a Search button. I would like the interface to find the polygon and zoom to it. With the research I have been doing, it appears I need to do this using a WMS layer. I have tried to import the layer as a WFS and it does not seem to be working.

My question is two part.

1) Do I need to do this with a WFS or can I stick with WMS?

2) Do you have any samples of code to do this with the data coming from this PostGIS store? I hate asking someone else to write the code, but if you can use what I put down below, that would be great just because I have been banging my head. I'm a .NET desktop programmer so this has been a learning curve.

Thanks for your help!


texts:
I have the following string in Arcmap:

"Node 1_1234.23,Node 2_567.88,Node 3_45.25"

I want to return the numbers between the "_" and the "," into three columns

e.g  "1234.23" in column A
     "567.88" in column B
     "45.25 in column C

Pls assist with an expression I can use in field calculator.


texts: I am trying to download QGIS Essen. I run the application and it says it has installed the program but it hasn't. I am using windows 7 and have previously downloaded QGIS successfully on another computer. Does anyone know why it would not be installing?


texts: Reading and storing settings in a QGIS project with QgsProject.instance()

texts:
I am trying to save full path to project settings and read it out and put it in the LineEdit function.

Save the path is working but to read it out and handle it with LineEdit doesn't?

Here is what I have done so far:

def select_output_file(self):
    filename = QFileDialog.getSaveFileName(self.dlg, "Selecteer output file ","", '*.ext')
    QgsProject.instance().writeEntry("mytools", "last_path", (str(filename)))
    self.dlg.lineEdit.setText(filename)

This works for me. When I read it out in the python console of QGIS I get this.

last_path = QgsProject.instance().readEntry("mytools", "last_path", "Output file")[0]
print last_path
D:/users/xxxxxx/downloads/filename.ext

So the path is correct, I think. But now I want to read it back like this next code:

last_path = QgsProject.instance().readEntry("mytools", "last_path", "Get last path")[0]
self.dlg.lineEdit.setText(last_path)

And here is where I don't get the path in the LineEdit field?

When I use this next code it works.

last_path = "C:/user/xxxx/"
self.dlg.lineEdit.setText(last_path)

And also this line works in python console.

last_path = QgsProject.instance().readEntry("mytools", "last_path", "Get last path")[0]


texts: Copy attributes from one feature to empty row of another feature in QGIS 3

texts:
How do you copy all of the non-geometry attribute data from one feature to another "empty" feature? For example, how would I copy all of the non-geometry attribute data from the "Seattle Mariners" polygon to the empty Alaska polygon? Using copy/paste pastes both the attributes and the geometry. I do not want to copy the geometry. Just the non-geometry attributes.

Here are instructions for how to do this in ArcMap. Is there a similar "copy attributes" option in QGIS?


texts: Extract Multi Values to Points gives ERROR 999999

texts:
I am trying to extract data from a raster based on coordinate points using the Extract Multi Values to Points tool but am repeatedly getting back ERROR 999999. I have ensured both layers are in the same coordinate system (GCS_WGS_1984)and I even clipped the raster file so it is smaller and easier to work with.
I saw that there was a bug with this tool when 64-bit rasters are used but I am using a 32-bit raster.
Any suggestions on what the issue could be?


texts: create border and labels for groups of features in arcmap

texts:
I am relatively new to GIS, so please bear with me if my question is a basic one.

I have a simple map of the parishes (counties) in Louisiana, and each of the parishes has an associated regional area which is denoted by a number. Using the symbology tab in properties, I was able to color each of the regions. However, now I need to do two things:

1) Embolden the borders between each of the regions to help make the distinctions more clear (while leaving the borders between the parishes as is)
2) Label each of the regions with their regional number. When I currently try to do this, a little number label pops up in each of the parishes rather than just one number for each region.

I am using Arcmap 10 on a windows 7 PC. Thanks for your help!


texts:
I am currently working locally on my terminal.
I have Windows 10 x64, 12Go RAM, 428Go SSD, an Intel(R) Core(TM) i5-6200U CPU @ 2.30GHz, 2400 MHz, 2 hearts.

I need to use PostGIS to make some queries and filter my data.

My database contains 39 tables. My biggest table has 4 millions rows.
I run some cross join with my tables.

The trouble I have is that my queries are taking too much time (maximum query took 6h and returned only null values).

I developed a Python script to split my queries and run them each by each instead of all of them in a row, but it also failed.

So I am looking for a way to have more power. I have understood so far that I should use a server to do this, but I do not know any server or service of this kind.

UPDATE 1:

I have created the indices for all my tables. Then I ran the query creating the most trouble, which is:

UPDATE table_1
SET field_1 = liste.my_bool
FROM(
    SELECT T1.id AS name,
        bool_or((ST_Contains(T2.geom, T1.geom) OR ST_Overlaps(T1.geom, T2.geom))) AS my_bool
    FROM table_1 AS T1
    CROSS JOIN table_2 AS T2
    GROUP BY T1.id
 ) AS liste
WHERE table_1.id = liste.name;

The table_1 has 52630 entities and the table_2 2025.
The query returned null values for my field_1 after one hour of loading. I then tried the subquery only and got this error message:

ERROR:  ERROR:  could not write the block 8568766 from temporary file : No space left on device
SQL state: 53100

I want field_1 to be True if the geometry of table_1 shares total or partial surface with any entity of table_2. This is why I need a CROSS JOIN.


texts: How to add legends in map canvas when exported to PDF ol5

texts:
Here is my code to export map in PDF and map is successfully exported to PDF
Now i want to add legends in map canvas but i do not know how to do that.

    var exportElement =  document.getElementById("export-pdf");
      exportElement.addEventListener('click', function(e) {

            exportElement.disabled = true;
            document.body.style.cursor = 'progress';

      var format = document.getElementById('format').value;
      var resolution = document.getElementById('resolution').value;
      var buttonLabelElement = document.getElementById('button-label');
      var label = buttonLabelElement.innerText;
      var scale = 8000;

    var mapView = map.getView();
    var mapProjection = mapView.getProjection();
    var mapResolutionAtEquator = mapView.getResolution();
    var viewCenter = mapView.getCenter();
    var mapPointResolution = ol.proj.getPointResolution(mapProjection, mapResolutionAtEquator, viewCenter);

    var dim = dims[format];
    var width = Math.round(dim[0] * resolution / 25.4);
    var height = Math.round(dim[1] * resolution / 25.4);
    var size = map.getSize();
    var extent = map.getView().calculateExtent(size);

    var mapView = map.getView();
    var mapProjection = mapView.getProjection();
    var mapResolutionAtEquator = mapView.getResolution();
    var viewCenter = mapView.getCenter();
    var mapResolutionAtEquator = mapView.getResolution();
    var mapPointResolution = ol.proj.getPointResolution(mapProjection, mapResolutionAtEquator, viewCenter);
    var mapResolutionFactor = mapResolutionAtEquator / mapPointResolution;
    var pdf_name = dist + "_" + tahsil + "_"+village+ ".pdf";
    map.once('rendercomplete', function(event) {

         buttonLabelElement.innerText = label;
          var canvas = event.context.canvas;
          var data = canvas.toDataURL('image/jpeg');
          var pdf = new jsPDF('landscape', undefined, format);
          pdf.addImage(data, 'JPEG', 0, 0, dim[0], dim[1]);
          pdf.save(pdf_name);
          map.setSize(size);
          map.getView().fit(extent, size);
          map.renderSync();
          exportElement.disabled = false;
              document.body.style.cursor = 'auto';
    });

    var printPointResolution = (scale * 25.4) /(resolution * 1000);
    var printResolutionAtEquator = mapResolutionFactor * printPointResolution;
    var printZoom = mapView.getZoomForResolution(printResolutionAtEquator);

    map.setSize([width, height]);
    mapView.setZoom(printZoom);
    map.renderSync();
   }, false);

see outout pdf looks too blur at 150 dpi or 300 dpi and even label is not clear.


texts: Querytask on multiple inputs - ignore when value is empty

texts:
I am trying to select some features by attributes. This function works correctly:

  function executeQueryTask() {
    query.where = "(xy.zzz.jprl_dop.dv = " + dom.byId("mySecond").value + ") AND (xy.zzz.jprl_dop.DC = " + dom.byId("myOne").value + ")";
    queryTask.execute(query, showResults);
}

But how to fix situation when one of partial query in function is not set ?

My HTML:

One: <input style="width:50px;" type="text" id="myOne" value="" />

                    Second: <select id="mySecond">
                        <option value="">-</option>
                        <option value="1">one</option>
                        <option value="2">two</option>
                        <option value="3">three</option>
                    </select>

Can someone advice me how to run query in my function when value=" " in option element ?


texts: Automatically assigning attributes when defining one attribute in table using QGIS?

texts: If it is possible to assign attributes when you define one attribute in a table, like automatic selection additional Attributes, then how is it done?


texts: Include MetadataURL (Metadata Link) in layer coverage programmatically

texts:
I'm using GeoServer's Java API library, and am trying to add add the MetadataURL to a coverage layer programmatically.  I can add it from within GeoServer's layer setup admin page.  I just don't see anywhere in the API where I can set it.

For example, this is what I see in the layer when I have the Metadata Link set.  How would I add this in using the Java API?

How the actual URL is made and it's content is trivial.

    <MetadataURL type="ISO19115:2003">
        <Format>text/plain</Format>
        <OnlineResource
         xlink:type="simple" xlink:href="http://localhost:8080/gogc/wcs?
         request=DescribeCoverage&version=2.0.1&service=WCS&CoverageId=HYCOM-
         GLOBAL:Region06.salinityWater.DBSL"/>
    </MetadataURL>


texts: Exporting multiple WCS layers to GeoTIFFs to use them in rastercalculator

texts: I would like to know if anyone could tell me in what way I can export multiple WCS files to GeoTIFFS in QGIS. I use WCS files, but these cannot be used in the raster calculator. I first need to export them and save them as GeoTIFFS for example, but since I have over 30 layers I would like to know if there is a quick way of doing this.


texts: Skip locked rows for editing using ArcObjects?

texts:
Can I lock a row so that the cursor does not bring it because it is in use by another cursor ?

ITable table=GetTable();
ICursor crntCrsr=null;
ICursor secondCrsr=null;
crntCrsr=table.Search(filter,false);

 // This cursor will contain one row. I want to lock it so that if
var firstRow= crsntCrsr.NextRow();

// now If I use the same statement as above, I will get the same row again.
// I want to skip this row because it is already being processed by crntCrsr
secondCrsr=table.Search(filter,false);
var secondRow= secondCrsr.NextRow();

// basically I do not want secondRow to be the same as firstRow


texts:
As my client software uses lat/lon coordinates when communicating with my (spherical mercator) postgis database i decided to ST_Transform every geometry to WGS-84.

However i noticed that ST_Distance for WGS-84 returns units as degrees (i need meters). So i decided to use ST_DistanceSpheroid(geom1, geom2, 'SPHEROID["WGS 84",6378137,298.257223563]') but this method seems to be extremely slow.

Therefor i switched back to using spherical mercator again and transforming my input and output to and from wgs-84 as that performs a lot better.

Am i using the correct method or is this a known issue?


texts: gdal_merge tries to create image with negative dimensions

texts:
I am trying to use gdal_merge.py to create a mosaic from a set of geotiffs, however the command fails with the error message: "ERROR 1: Attempt to create -3510x1817 dataset is illegal,sizes must be larger than zero.".  I have tried with gdalbuildvrt, but the .vrt also has negative dimensions.

The source geotiffs are not north-south aligned, but instead were generated from .jpgs using

gdal_translate -a_srs WGS84 -a_ullr tlx tly brx bry src dst

with the bounds calculated from the (known) image footprint extents, and so have a range of different orientations.

The geotiffs I'm trying to mosaic are here: https://www.dropbox.com/s/hvhhpsrb913l3xb/geotiffs.zip?dl=0

gdalinfo on a sample geotiff:

Driver: GTiff/GeoTIFF
Files: geo_pict20140910_131133_0.tif
Size is 1616, 1080
Coordinate System is:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]
Origin = (-121.459108654999994,38.582958816199998)
Pixel Size = (-0.000000523938738,-0.000000765216574)
Metadata:
  AREA_OR_POINT=Area
  EXIF_BrightnessValue=(4.22031)
  EXIF_ColorSpace=1
  EXIF_ComponentsConfiguration=0x1 0x2 0x3 00
  EXIF_CompressedBitsPerPixel=(3)
  EXIF_Contrast=0
  EXIF_CustomRendered=0
  EXIF_DateTime=2014:09:10 13:11:32
  EXIF_DateTimeDigitized=2014:09:10 13:11:32
  EXIF_DateTimeOriginal=2014:09:10 13:11:32
  EXIF_DigitalZoomRatio=(1)
  EXIF_ExifVersion=0230
  EXIF_ExposureBiasValue=(0)
  EXIF_ExposureMode=0
  EXIF_ExposureProgram=4
  EXIF_ExposureTime=(0.001)
  EXIF_FileSource=0x3
  EXIF_Flash=16
  EXIF_FlashpixVersion=0100
  EXIF_FNumber=(4)
  EXIF_FocalLength=(16)
  EXIF_FocalLengthIn35mmFilm=24
  EXIF_GPSAltitude=(67)
  EXIF_GPSAltitudeRef=00
  EXIF_GPSImgDirection=(95)
  EXIF_GPSImgDirectionRef=T
  EXIF_GPSLatitude=(38) (34) (57.16)
  EXIF_GPSLatitudeRef=N
  EXIF_GPSLongitude=(121) (27) (34.32)
  EXIF_GPSLongitudeRef=W
  EXIF_GPSVersionID=0x2 0x3 00 00
  EXIF_Interoperability_Index=R98
  EXIF_Interoperability_Version=0x30 0x31 0x30 0x30
  EXIF_ISOSpeedRatings=2500
  EXIF_LightSource=0
  EXIF_Make=SONY
  EXIF_MaxApertureValue=(2.96875)
  EXIF_MeteringMode=5
  EXIF_Model=NEX-5T
  EXIF_Orientation=1
  EXIF_PixelXDimension=1616
  EXIF_PixelYDimension=1080
  EXIF_ResolutionUnit=2
  EXIF_Saturation=0
  EXIF_SceneCaptureType=0
  EXIF_SceneType=0x1
  EXIF_Sharpness=0
  EXIF_Software=NEX-5T v1.01
  EXIF_WhiteBalance=0
  EXIF_XResolution=(350)
  EXIF_YCbCrPositioning=2
  EXIF_YResolution=(350)
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (-121.4591087,  38.5829588) (121d27'32.79"W, 38d34'58.65"N)
Lower Left  (-121.4591087,  38.5821324) (121d27'32.79"W, 38d34'55.68"N)
Upper Right (-121.4599553,  38.5829588) (121d27'35.84"W, 38d34'58.65"N)
Lower Right (-121.4599553,  38.5821324) (121d27'35.84"W, 38d34'55.68"N)
Center      (-121.4595320,  38.5825456) (121d27'34.32"W, 38d34'57.16"N)
Band 1 Block=1616x1 Type=Byte, ColorInterp=Red
Band 2 Block=1616x1 Type=Byte, ColorInterp=Green
Band 3 Block=1616x1 Type=Byte, ColorInterp=Blue


texts: Histogram of ratio between different pixels in Google Earth Engine

texts:
I would like to generate a histogram of the ratio between pixels within a polygon called "MODISofstation". Thereby all the pixels within my polygon should be compared to the pixel where my Point "pt" is located in and I can see the variability within the polygon. I am only looking at B8 for this example.
The first part of my code works fine and I receive all the values from each month but from the part "//RATIO" on, I fail to make it work to calculate the monthly mean, the ratio and bring it into a histogram.
Below you can find the code snippet

var start = '2020-06-01';
var finish = '2020-06-30';
var pt = ee.Geometry.Point([-49.31582, 69.56833]);
// Map the function over one year of Landsat 8 TOA data and take the median.
var l8 = ee.ImageCollection('LANDSAT/LC08/C01/T1_TOA')
.filterDate (start,finish)
.filterBounds(pt)
.filter("WRS_ROW <= 122")
.map(maskL8);

//Map. addLayer(B8)
//Map. addLayer(pt)
Map.centerObject(pt, 16);

print(l8);

// B8 VALUE OF PT

var myB8 = l8.select("B8")
print("myB8",myB8);

var getB8 = function(image) {

  // Reducing region and getting value
  var value_B8 = ee.Image(image)
    .reduceRegion(ee.Reducer.first(), pt)
    .get('B8');

  return value_B8;
};

var count = myB8.size();

var B8_list = l8.toList(count).map(getB8);

print("B8 list", B8_list);

var img = l8.first();

Map.addLayer(img);
Map.addLayer(pt);

var allDates = l8.aggregate_array('system:time_start');

var allDatesSimple = allDates.map(function procDates (ele) {

  return ee.Date(ele).format().slice(0,10);

});

var paired = allDatesSimple.zip(B8_list);

print (paired);

// RATIO

//Monthly mean of pt pixel
var MonthlyMeanPt = myB8.mean;
print (MonthlyMean)

var ratio= ee.Image.divide(Monthlymean)
Map.addLayer (ratio);

//Histogram within my polygon MODISofstation

var histogram = ui.Chart.image.histogram({
  image: ratio,
  region: MODISofstation,
  scale: 15,
});
histogram.setOptions({
  title: 'Histogram of Ratio'
});

print(histogram);


texts: I have a shapefile with polygon type geometry. I want WKB format of the geometry values. Is there any such tool in QGIS so that without using gdal/ogr I can do this?


texts: I have raster image that projected by WGS1984 UTM Zone 37N. But my shapefile was Adindan UTM Zone 37N. Due to this, it doesn't overlay the same. I tried to define projection for the raster image to Adindan UTM zone 37N, But not match. So, how can I display the raster image in Adindan UTM Zone 37N without change WGS1984 UTM Zone 37N?


texts: How to remove  text from legend in layout view?

texts: I'm producing a map in the layout view which I will eventually export as a jpeg. I've added a legend which includes 15 classes but when I insert a legend in layout view its including the word VALUE at the top of the legend column - I had a look at the properties and can't seem to figure out how to remove this - I'm using ArcMap 10.


texts:
I am trying to automate a process where I am converting point shapefiles to polygons in QGIS using Points2One. I am fairly new to Python coding and have not been able to find any example code online of this being done.

Is it possible to do this using Python code?

I am using QGIS version 2.14.


texts: Finding crossing lines in ArcGIS Desktop?

texts:
I have a network of lines (pipes) and I want to be able to show the difference between an intersecting line and a crossing line.

I was wondering if there is a tool for generating points on all crossings in ArcGIS?

If I use "Intersect" I can find all intersecting and crossing points. I only need the crossing points.


texts: Automating point feature class placement inside polygons using specific parameters?

texts:
I am trying to generate a point feature class of possible plot centers for a research project. The plot centers will be inside of defined forest stand polygons which are aggregated in various ways across the landscape. What I am trying to achieve is this: approximately (+/- 100 m) equidistant points within and between forest polygons, so each polygon has two points inside and then one or both points are the same distance from a point in a separate polygon so that each pair looks something like this (the green/orange shapes are supposed to be two forest polygons, the blue points are the plot centers, and the red lines are just to show equal distance between points-I don't  need to actually create a line class connecting the points, the euclidean distance between them just needs to be similar.

I would also like to be able to define a minimum buffer distance between the points and the polygon edges. I have lots and lots of polygons, but for most of them this configuration will not be possible, so I'm looking for a way to identify potentially usable polygon pairs. I have access to ArcGIS Pro and I use ArcMap over Citrix if that's relevant. I have limited experience with Python but am not afraid to learn if that's the best solution.


texts: Is there any way to convert/read/import Geosoft Montaj grid in Mac OS?

texts:
I'd like to read data from a .grd file, made by Oasis Geosoft Montaj.

There is a free viewer that can convert to an ascii file, that can be read in e.g. QGIS, but the viewer program only works in Windows.

Are there any options to import Geosoft grids on a Mac? Preferable free, open source.

(It is not the same format as e.g. GMT .grd, which is a kind of netCDF.)


texts: Saving MultiLineString to Shapefile and KML/KMZ with Python

texts:
What is the fastest way to save MultiLineString to Shapefile (without attributes) and KML/KMZ? Preferred open-source libraries like fiona, shapely etc.

from osgeo import ogr

multiline = ogr.Geometry(ogr.wkbMultiLineString)

for (X1, Y1, X2, Y2) in coordinates:
    line = ogr.Geometry(ogr.wkbLineString)
    line.AddPoint(X1, Y1)
    line.AddPoint(X2, Y2)
    multiline.AddGeometry(line1)

When I print the multiline on the console:

print(multiline.ExportToWkt())

I get something like this:

MULTILINESTRING ((10.0 15.0 0,10.5 15.5 0),(20.0 25.0 0,20.5 25.5 0))

I would like to save these two lines to Shapefile and KML/KMZ as two separate lines.


texts:
I'm new to QGIS and Maperitive. I wondered if there was a way to export all road /street names for a particular area as editable text in Illustrator (SVG)?

I have tried using Maperitive but for some reason when I export to SVG most of the words that I can see when zoomed in in Maperitive are either lost or missing letters.

Is there perhaps a mrules I could use that doesnt lose all the words?

Maybe there is a way to do this is QGIS with a particular plugin or layer?


texts: Any experience with OSX and Mac OSX Maverics using OSGeo?

texts:
There is a new Apple Release of OSX called Mavericks and my questions are:

Will GDAL run?
Will QGIS run?
Will Tilemill run?
Will Mapnik run?
Also are there any known issues?

A similar question about Mountain Lion.


texts: How to calculate mean value of a raster for each polygon in a shapefile

texts:
I wrote a Python script who calculate, for rasters, mean value for each polygon of a shapefile. But I think that it's a "war machine" (it takes many treatment time) and I would like to simplify it.

Is somebody know a simple way to do it?

I'm using geopandas for shapefile manipulations.


texts:
I am trying to do conditional statement in ENVI to calculate multiband image:

result = (b1)^2 - (b2)^2
if (b1 eq 0) or (b2 eq 0)
result * 0
else
result * (result gt 0.0)

Is there a way to combine the statements in Band Math? Or I have to write in IDL?


texts: Leaflet create marker from ajax call data

texts:
Get all the basic works with Leaflet and want to add markers using ajax call

    $("button").click(function(){
    $.get("GetLocationsServlet", function(data){
        markers = data.split(";");
        var length = markers.length;

        for (var i = 0; i < length; i++) {

            alert(markers[i]); // data print out ok
            L.marker(markers[i]).addTo(mymap); // Error
        }
    });
});

Got error when adding the marker to the map

Uncaught TypeError: Cannot read property 'lat' of null
at Object.project (leaflet.js:5)

The data return from the ajax call is like this

[49.205036,-123.002846, "Riverway Sport Complex"];[49.22193,-123.01249, "Imperial@Patterson"];[49.214756,-122.9882, "Royal Oak @Rumble"]

It worked when I tried to manually add the marker by hard coding the data to the L.marker() function. i.e.

L.marker([49.205036,-123.002846, "Riverway Sport Complex"]).addTo(mymap);


texts: I have two polygon layers:  one indicating countries, and other on top indicating lakes. I want to create a new layer that keeps the countries but have the lakes cut out of them. What tool do I use?


texts: How to get rid of double labels on ArcGIS 10

texts:
I'm displaying several point shapefiles.  I want them all labeled with the field "ident" which is a 3 digit number starting at 001.  Two of the three shapefiles show the labels correctly but the 3rd one shows a double label for each point.  Instead of just saying "027" it says "027027".  In the Layer properties window, I tried unchecking the box that says "Label features in this layer" but this removes the labels completely, then when I recheck "Label features" in the menu that comes when I right click the layer in the table of contents, the double layers re-appear.

Any suggestions would be much appreciated!

Molly


texts: Calculate slope from ALOS Global Digital Surface Model

texts:
I have downloaded the relevant tiles from this model for the area I am studying, but cannot create meaningful slope data using the "Slope" tool of ArcGIS 10.5

I understand that the data requires projecting - I think the correct procedure is to use the UTM or Lambert projections due to the units being in arcseconds, however my study area covers multiple UTM grids leaving only Lambert an option (I think).

I tried projecting the raster from WGS 84 to "Europe_Lambert_Conformal_Conic" using the Coordinate Transformation: ED50 to WGS 84 (1), but this resulted in massively blown out of proportion data, perhaps due to an incorrect cell size? I used "Bilinear" over the "Nearest" resampling technique but didn't specify a cell size as I'm not sure how to work this out.


texts: assign_vertex_id：problems when creating the node ids

texts:
Before using pgrouting, I have to create the node ids.
I loaded my polyline layer into PostGIS Manager：

I type this in SQL query：

select assign_vertex_id('split_4', 800, 'geom', 'gid')

But as I click the "Execute" button, I got this error：

Error: Function StartPoint (geometry) does not exist
LINE 1: SELECT gid AS id, StartPoint (geom) AS source, EndPoint (geom) ...
                          ^
HINT: did not comply with the specified name and parameter types of function. You may need to add explicit type conversion.
QUERY: SELECT gid AS id, StartPoint (geom) AS source, EndPoint (geom) as target FROM split_4
CONTEXT: PL / pgSQL function "assign_vertex_id" line 23 at target EXECUTE statement FOR

I searched for ST_StartPoint, and found that it needs working with linestring, so I cut my multilinestring to a linestring(named it "ways"), and then run the SQL：

ALTER TABLE ways ADD COLUMN source integer;
ALTER TABLE ways ADD COLUMN target integer;
SELECT assign_vertex_id('ways', 800, 'the_geom', 'gid');

but got the same error above.

How to solve this and to create the node ids, then I can proceed to the next step in A Beginner’s Guide to pgRouting？I'm working with pgadmin3, postgresql8.4, PostGIS2.0 under WIN7-64bits.


texts: What is a good metric for mapping accuracy

texts:
I have recently received a dataset that supposedly mapped the altitude a given region (a DEM).  I would like to determine the accuracy of the new map by comparing it with a previously generated DEM.  I am not sure what kind of metric I can use to ensure that the performance is good enough.  (Good enough means that the horizontal deviation of a feature can not be off too far and that the vertical deviation can not be off too far)

One way to attack this would be to descritize the space into small chunks and then look at the distribution of the vertical residuals.  If that block has vertical residuals that are less than the vertical deflection then I can assign that block a TRUE value and if it does not I could assign the block a FALSE value.  I can then simply calculate the mean of the TRUE/FALSE values across all the blocks to get an overall value such as 67% or something which would tell me how good my mapping performance was.

There are huge problems with this method:

It does not seem to factor in the horizontal deflection requirement
It seems as though it would violate the iid (identical independently distributed) assumption normally required of the mean calculation.  This is the same as stating that the individual boxes would be correlated.
There does not seem to be a defensible method to choose how big the descritization box size should be.

An alternative could be that I simply place special markers with known 3D locations and then simply measure their vertical and horizontal displacement.  This seems like it could share some of the same problems as the earlier method.  (how many markers should I place?  are the measurements correlated if they are too close?)

Before running off and doing something silly, I wanted to see if the community has a golden metric or method that has good statistical properties.   Something that is continuously measured across the DEM, takes into account both the vertical deflection and the horizontal deflection, and I can look at the distribution and make a defensible argument?  Please let me know if this question is unclear and thanks for any suggestions.

P.S. Could somebody could comment on the need of the two individual mapping accuracy levels (one for vertical deflection and one for horizontal deflection)?  It seems to me that a horizontal mis-mapping would be reflected in a vertical deflection.  Is it standard to have two separate requirements for mapping?


texts: PostGIS 2.0 function returns pseudo-type record rather than table

texts:
I'm writing a PL/pgSQL function to take a geography column and output a fishnet that bounds the extent of the shapes.

I've tried to emulate the examples in docs and tutorials, but the function returns a pseudo-type result, not a usable table.

CREATE OR REPLACE FUNCTION ST_Fishnet(
        nrow integer,
        ncol integer,
        geo_table varchar,
        geom_column varchar)
RETURNS Table( ids text, geom geometry)
AS $$
DECLARE
    collected_geom geometry;
    xmin numeric;
    xmax numeric;
    ymin numeric;
    ymax numeric;
    xsize numeric;
    ysize numeric;
BEGIN
    EXECUTE 'SELECT st_collect("'||$4||'") FROM "'||$3||'"' INTO collected_geom;
    xmax := st_xmin(collected_geom);
    xmin := st_xmax(collected_geom);
    ymin := st_ymin(collected_geom);
    ymax := st_ymax(collected_geom);
    xsize := abs(xmax - xmin) / $2;
    ysize := abs(ymax - ymin) / $1;
    return query
        SELECT
            i+1||'-'||j+1 as ids,
            ST_Translate(cell, j * xsize + xmax, i * ysize + ymin) AS geom
        FROM
                generate_series(0, $1 - 1) AS i,
                generate_series(0, $2 - 1) AS j,
                (SELECT ('POLYGON((0 0, 0 '||ysize||', '||xsize||' '||ysize||', '||xsize||' 0,0 0))')::geometry AS cell) AS temp;
END
$$ LANGUAGE plpgsql;

But while my output should have two distinct columns, ids and the geom, I'm getting pseudo-type record that looks like:
(id, geom)

Very much appreciate any pointers as to what I'm missing.

Thanks.


texts: Generating and hosting vector tiles from GeoJSON - no errors, but no data

texts:
I wish to construct a pipleline that takes GeoJSON through to self-hosted MapBox vector tiles specification.

My plan is:

Tippecanoe converts geojson to .mbt
mbutil converts mbt to tiles
tileserver-php exposes endpoint for tiles.

This seems to be possible, given the docs. I have put it all together. No errors are thrown by any process and tileserver is showing that my tileset exists. However, there is an empty preview of the data and HTTP GET requests to the Z/X/Y endpoint return an HTML page (Welcome to Maps hosted with TileServer-php v2.0). (Perhaps I need to curl with specific headers in order to test if it's working?)

My GeoJSON is an array of linestrings & is structured like this:

{
  "type": "FeatureCollection",
  "bbox": [-119.437,33.298,-116.724,34.583],
  "features": [
    {
      "type": "Feature",
      "geometry": {
        "type": "LineString",
        "coordinates": [
          [-117.9756412,33.803147],
          [-117.9775983,33.8031329]
        ]
      },
      "properties": {
        ...
      }
    }, ...

Using tippecanoe - appears to be OK:

$ tippecanoe -o out.mbt in.json
For layer 0, using name "my_layer"
99641 features, 2896256 bytes of geometry, 2482729 bytes of metadata, 1532578 bytes of string pool
tile 5/5/12 size is 674507 with detail 12, >500000
tile 6/10/25 size is 537511 with detail 12, >500000
tile 6/11/25 size is 608121 with detail 12, >500000
tile 7/21/51 size is 665802 with detail 12, >500000
tile 7/22/51 size is 685644 with detail 12, >500000
tile 7/22/51 size is 538985 with detail 11, >500000
tile 8/44/102 size is 754537 with detail 12, >500000
tile 8/43/102 size is 750288 with detail 12, >500000
tile 8/44/102 size is 652762 with detail 11, >500000
tile 8/43/102 size is 633573 with detail 11, >500000
tile 8/44/102 size is 513572 with detail 10, >500000
tile 9/88/204 size is 518065 with detail 12, >500000
tile 9/87/204 size is 765483 with detail 12, >500000
tile 9/87/204 size is 697630 with detail 11, >500000
tile 9/87/204 size is 583396 with detail 10, >500000
  100.2%  14/2773/6531

Tileserver docs say it works better if I run mbutil over the mbt file first. I do this with the command mb-util out.mbt tiles. There's a lot of output in the format similar to:

flipping
INFO:mbutil.util:3774 / 3774 tiles exported

Afterwards, the directory tiles has a 3-levels deep nested structure with lots of files ending with png. Despite the extension, the OS cannot render these as raster images. But given what they're intended to do I assume they're not PNGs as I'm used to and that this is expected.

Finally, I push the tiles along with the Tileserver code onto a PHP-enabled host in ElasticBeanstalk. The app runs and presents me with my tileset, rendered as a completely empty (black) MapboxGL view.

Here, let me show you.

So it seems nothing threw any errors, but I don't have any data. How can I find where this all went wrong?


texts: Count Blanks, Spaces, and Hardcoded Nulls

texts:
I want to get counts of  Blanks, Spaces, and Hardcoded Nulls in fields using arcpy.SearchCursor. Is it possible? When I try to filter am getting None values only.
I tried the below code.
 cursorfields = arcpy.ListFields(featureclass_path)
    for field in cursorfields:
        cursor = arcpy.SearchCursor(featureclass_path)
        for row in cursor:
            if (str(row.getValue(field.name)) == "Null"):

            if(str(row.getValue(field.name)) == "NULL"):

            if(str(row.getValue(field.name)) == "<Null>"):

            if(len(str(row.getValue(field.name))) == 0):

            if(not(str(row.getValue(field.name)) and str(row.getValue(field.name)).strip())):


texts: How can I merge two lines to one via OpenLayers?

texts:
I have a tool on the map that lets the user draw a line.
Now I want to give the user the option to split the line to multiple lines, to join them back together to one line later.

How can I do that ?


texts:
how can i identify orientation of a building in reference to another rectangle (eg building located  at North side or south side in reference with campus).

in FME there is a way to find out?


texts: How to change spatial resolution using gdalwarp?

texts:
I just realized that I have Geo-TIFF's in the same coordinate system but different spatial resolution. DOP20 (20cm=1px) and DOP10 (10cm=1px) mix.
Now wanted to use gdalwarp to transform the DOP20 TIFF into DOP 10.
The commandline should be simple

gdalwarp.exe -tr 20000 20000 -r bilinear in.tiff out.tiff

but gdalwarp complains:

Creating outputfile that is 0Px0L Error 1 Attempt to create 0x0 dataset...

I thought it is a simple graphics operation, "blowing" a 10000x10000 pixel TIFF into a 20000x20000 pixel TIFF.

Can anybody help me understand what's wrong here?


texts: How to query aggregate max, min, median latitude of polygons in PostGIS?

texts:
I am just starting to learn how to use PostGIS (1.5.4) and my SQL skills are essentially those of a competent beginner. Having perused the PostGIS documentation and searched these forums, I'm not finding the answers to what are surely some very simple questions.

Given a single table loaded with polygons from shapefiles, how can I query the latitudinal maximum and minimum and the centroids of single and aggregate rows? I'll note that I have found the ST_Centroid() function, but haven't figured out how to interpret its output (CRS=WGS84).


texts: Rule for labeling ways not being part of OSM relation in QGIS

texts:
A brief overview of what I want to achieve: Create a map showing all cross-country ski trails in a specific region with some customizations/styling on them.

The data come from OpenStreetMap exclusively. Ski trails can be either relations or ways.

I imported the ways and relations from OpenStreetMap using QuickOSM plugin in QGIS 3 with the following Overpass query:

<osm-script output="xml" timeout="25">
    <union>
        <query type="way">
            <has-kv k="piste:type" v="nordic"/>
        </query>
        <query type="relation">
            <has-kv k="piste:type" v="nordic"/>
        </query>
    </union>
    <union>
        <item/>
        <recurse type="down"/>
    </union>
    <print mode="body"/>
</osm-script>

Now I have 2 layers in QGIS, one for lines (OSM ways) and one for multiline strings (OSM relations).

What I want to do is to add a label on ways (showing the value of the attribute piste:name or, if it doesn't exist, name) only if those ways are not part of a relation from the second layer. This is to avoid duplicating the label since these relations already have their name displayed with another labeling rule.

I tried to write rules using geometry functions like contains() but it seems that I can't reference an object from another layer.

Is there a way to achieve what I want to do in QGIS?


texts: Table not included as layer with result map service

texts:
I published a geoprocessing service that outputs (1) a feature layer and (2) a stand alone table. The goal is to join these two on the client side... given the size of the table, it takes too long to join them on the server and use "Copy feature" to create a permanent layer.

According to this documentation, GPRecordSet data types should be included in the resulting map service (I am assuming a table is this type). However, it is not included after a successful run of the task. Here's subset of the code:

result_layer_var = arcpy.MakeFeatureLayer_management( geography, "Result_layer" )

###... Do some spreadsheet work using Pandas ...###

result_spreadsheet_path = os.path.join( scratch_folder, "Result_spreadsheet.csv" )
pandas_table.to_csv( result_spreadsheet_path )
result_table_var = arcpy.TableToTable_conversion( result_spreadsheet_path, scratch_gdb, "Results_table")

arcpy.SetParameter(3, result_layer_var)
arcpy.SetParameter(4, result_table_var)

Any thoughts on why the table isn't being included in the resulting map service?

More infomation:

GP settings

Using "Table View" As parameter type

Using "Table" As parameter type

Results from both -- notice that E_Output (table) is not included in map service

Why is the GPRecordSet (table) subject to the "exceededTransferLimit: true" but the GPFeatureRecordSetLayer (feature layer) is not... they both have the same number of records.


texts:
I try to make maps of the different administration units in Denmark but have problems with finding the correct CRS.

I downloaded GeoJSON files from this repository.

When I import one of those to QGIS, it looks like this

But I expected it to look more like this. (I mean the general form, not the colours or borders.)

I did some research which CRS would be best for Denmark, and in this list, I found out, that "ETRS_1989_DKTM1" should be fine. But unfortunately, I still don't get the result I want.

What CRS should I use? Or is it something else that I did wrong?


texts: Updating QGIS version through OSGEO4W - Unable to get setup.ini

texts:
I want to update the QGIS versions that I am using through OSGEO4W. I am currently running QGIS 2.0.1-Dufour (and also have an QGIS 191a229 version) on Windows 7. However, when I go through the advanced installation in OSGEO4W Setup and either select "Direct Connection" or "Use IE5 settings" I get the following message: Unable to get setup.ini from http://download.osgeo.org/osgeo4w.
This is not the first time I have updated QGIS versions and I have always gone through "Direct Connection" without problems - this message has never popped up before. As a result, no progress is made and I have no other option than to shut OSGEO4W. Does anyone know why this is happening and what is the best way around this?
I am thinking of either installing a standalone version of QGIS or deleting the whole OSGEO4W folder and starting over again.


texts: WMTS from GeoServer to OpenLayers is Unknown TILEMATRIX X

texts:
I tryed use default WMTS layer with GeoServer source.
Code looks like that:

      var projection = ol.proj.get('EPSG:3857');
      var projectionExtent = projection.getExtent();
      var size = ol.extent.getWidth(projectionExtent) / 256;
      var resolutions = new Array(10);
      var matrixIds = new Array(10);
      for (var z = 1; z < 10; ++z) {
          // generate resolutions and matrixIds arrays for this WMTS
          resolutions[z] = size / Math.pow(2, z);
          matrixIds[z] = z;
      }
      const wmts_osm = new ol.source.WMTS({
        url: 'http://local:8080/geoserver/gwc/service/wmts',
              layer: 'workspace:test_CFD_layerset',
              matrixSet: 'geogrid4CFD',
              format: 'image/png',

              tileGrid: new ol.tilegrid.WMTS({
                        origin: ol.extent.getTopLeft(projectionExtent),
                        resolutions: resolutions,
                        matrixIds: matrixIds
              }),
              style: 'default',
              wrapX: true
      })

      var wmtslayer = new ol.layer.Tile({
        source: wmts_osm
      })

but WMTS requests in Chrome debuger are 400:

<ExceptionText>Unknown TILEMATRIX 6</ExceptionText>

What is the problem?


texts: Load PostgreSQL table with composite primary key to QGIS

texts:
I would like to load a PostgreSQL table without geometry and with a composite primary key in QGIS (there is no unique column, see below).
This yields “ERROR” in all cells. It is not resolved by adding a SERIAL column.

The same occurs when loading using DB Manager. But the issue is resolved when I “Edit table” in DB Manager, using “SELECT * from table” and “Load as new layer”. If I select all records, DB Manager adds a unique _uid_ column. If I select one subgroup from (id_name, sub_id), I can select the time_ column as unique values which also works.

What am I missing for QGIS to load the table directly from PostgreSQL? Is a unique column as primary key required, or should I look elsewhere?

I am using a PostgreSQL/PostGIS database v9.4.6 and QGIS 2.14.9 ltr.

Table info:
Fields:
Name  - Type  - Length  - Null
id_name  - varchar - (30)  -  N
sub_id  - varchar - (30)  -   N
time_  - timestamp  - 8  - N
mtaw  - float8  - 8  - Y
mmv  - float8  - 8  - Y

Constraints:
Name  - Type  - Column(s)
timeseries_values_pkey  - Primary key  - id_name, sub_id, time_
timeseries_values_id_name_fkey  - Foreign key  - id_name, sub_id


texts: Compute zonal statistics for portion of polygons

texts:
I am working with a shapefile containing the district polygons of the political map of India and a raster file containing the elevation data of the same country. My goal is to compute the average and variance altitude within 10 km from the borders of each district. Using Zonal Statistics I was able to compute the statistics I am looking for at the district level, the further step I need is disregard the pixels lying in the central part.

I have tried to get rid of the internal parts of the districts with add rings, but the user should draw the internal polygons by himself. This is not what I am looking for: I want to have an internal polygon with exactly the same angles as the original one but with smaller sides, computed automatically for all India.

How can I create, alternatively, a mask in the shapefile for the area that I don't want (or that I want, the one that most fits the case) to take into account and then use it to compute mean and variance with the Zonal Statistics plugin?

The ideal answer would provide me a PyQGIS code but it's fine also using icons in QGIS.


texts:
When scripting with arcpy.GeocodeAddress_geocoding (http://resources.arcgis.com/en/help/main/10.2/index.html#//002600000006000000), there's always a major lag to start up the geocoder.  The issue here is that even geocoding 1 address typically takes 2-3 minutes.

Since the function just takes the path to the locator, I'm not seeing any way to pre-load the locator into memory (this is command line or gp service, not within arcmap).

Environment:

ArcGIS Server 10.2.2 64-bit python called from command line
Navteq Street Map Premium 2013 - USA_Zip4.loc (composite locator)
Same results on several very fast machines, so computer's not really the issue here

Looking for workarounds or creative solutions here.... ArcObjects solution available to python?  Changing settings on locator?  Would prefer non-service based solutions (publish as geocode service and call).


texts: Add NZGD2000 Coordinate Reference System

texts:
I'm not sure if this is the right place to ask, but...
I've been playing with GeoTools 8.0M-1, in particular trying to read some GeoTiffs. Each time, it crashes with the error No code "EPSG:2193" from authority "EPSG" found for object of type "EngineeringCRS". The GeoTiff I'm reading is an NZGD2000 map. A Google for NZGD2000 and 2193 finds http://resources.esri.com/help/900/arcgisexplorer/sdk/html/612c4dad-84c2-7ec5-2e23-7db2bdfecc7d.htm ?

How can I get that information into the GeoTools, so that it will carry on and read the image via the GeoTiffReader class?

Thanks

Frank


texts: Selecting features whose text is equal to the text after certain word in QGIS

texts:
I have origin-destination lines and I'm looking for intra-city lines or lines whose origin and destination are in the same city.
A sample intra-city feature would be "Quezon_City to Quezon_City" and inter-city would be "Quezon_City to CALOOCAN_CITY".
I'd like some help in writing an expression that selects all features whose string before "to" is the same with the string after it using Python or Expression Builder.


texts: Loop through Features in a Layer with ArcPy

texts:
I am trying to loop through a list of features in a layer to:

select a feature
select all features adjacent to that feature
export the selection as a new table, then convert it to an excel file.

I have the code to do this for an individual selection:

# Select a zip 3 zone
arcpy.SelectLayerByAttribute_management("test", "NEW_SELECTION", """ "zip3_str" = '432' """)
# select all the adjacent features to that zip 3 zone
arcpy.SelectLayerByLocation_management("test", "BOUNDARY_TOUCHES")
# MakeTableView_management (in_table, out_view, {where_clause}, {workspace}, {field_info})
arcpy.MakeTableView_management("selection", "zip3")
# convert it to an excel file
arcpy.TableToExcel_conversion("zip3", r"J:\Strategic and Org. Dev. Department\PantryTrak\ARCGIS\Projects\20160921 - Zip 3 Relationship File\zip3.xls")

This is how I am attempting to write the loop but it is not working:

cursor = arcpy.SearchCursor("test")
count = 432

for row in cursor:
    arcpy.SelectLayerByAttribute_management("test", "NEW_SELECTION",""" "zip3" = count """)
    arcpy.SelectLayerByLocation_management("test", "BOUNDARY_TOUCHES")
    arcpy.TableToExcel_conversion("zip3" , r"J:\Strategic and Org. Dev. Department\PantryTrak\ARCGIS\Projects\20160921 - Zip 3 Relationship File\\")
    count = count + 1

The program gives me an error for the line:

arcpy.SelectLayerByAttribute_management("test", "NEW_SELECTION",""" "zip3" = count """).

It says """ "zip3" = count """ is an invalid expression, although it works without the variable. I am not sure why it doesn't work when I try to use a variable.


texts: What is the overall architecture of ArcSDE, and how does it fit into the enterprise?

texts:
I'm finding many gaps in my understanding of ArcSDE, particularly as it relates to Oracle.

I just can't quite wrap my brain around the overall architecture of it, and how it fits into an enterprise. Is ArcSDE installed on the Oracle server or on every client? Both?


texts: Mosaic UAV image step by step using ArcMap?

texts:
I have 20 pieces of UAV image, and using ArcMap i want to mosaic it all to become one raster but step by step,

what i mean is mosaic image 1 & 2, image 3 & 4 ...... then result of 1&2 mosaic with result of 3&4......, then mosaic them all (i do this to prevent crashing so my time didnt wasted).

Should i build statistic and pyramid after every mosaic step, or i can just skip build pyramid and stat for every mosaic i build, and just build mosaic on the last step (when mosaic all of it)?


texts:
I am having some trouble creating a model that will iterate through folders and subfolders. I want to be able to have the top level folder as a parameter and then have the code go through each of the sub folders and complete an action.

I have working code for the action I'm looking to complete and I am able to create an iteration to go through one folder, but I want to be able to go through 18 folders and run my script on every .mxd within each folder.

Here is the model that runs within one folder.

I want this to run through the 18 folders that I have.

Any ideas? I can run this either as a model or a script.


texts: ESRI Proportion tool/command and ArcFM split configuration

texts:
In my project we are using ArcGIS 10.1 and ArcFM 10.1.
I have requirement that

"If customer use Proportion Tool to split abandon liner feature then action need to assign new GISID (attribute) to shorter segment(s) only. (none of other attributes should need to be updated)"

To satisfy this requirement I should need to build new AutoUpdater.

Currently I know if we use ESRI split Tool then it's fires split events of ArcFM (BeforeFeaturesplit, OnFeaturesplit, OnFeaturecreate(for new segement), AfterFeatureSplit).

but While using the Proportion Toolit's just fires only
OnFeatureCreate event.


texts:
I'm wanting to georeference an old chart in the same manner described in a published article:

"All charts were transformed to the British National Grid. The 1904 chart was georeferenced using nine fixed landmarks such as church spires, of which the National Grid coordinates were known. Errors due to georeferencing were <0.01 cm, corresponding to a maximum positioning error of 7.5 m on 1:75,000 maps. The accuracy of the digitizer itself was claimed to be +/- 0.127 mm."

My own image has a cell size of 22, has a scale of approx. 1:127,000 and was georeferenced onto a map of scale 1:25,000 using 10 evenly-distributed GCPs. My Total RMS Error is 36.

How do I identify the maximum positioning error as described in the article? I.e. can I convert RMS Error into meters?

Interpretation of RMS Error is well described in this thread:
Generally accepted root mean square (RMS) error for rectifying topographic maps. However, I'm not sure how the author obtains an "RMS of 5 m". Is the RMS Error expressed in meters, or is there an additional conversion needed?


texts: What is the best way to visualize subsurface drillhole and borehole geology data within ArcGIS?

texts:
I've heard some folks in mining are using Target for ArcGIS extension by Geosoft to
visualize and analyze subsurface drillhole and borehole geology data within ArcGIS.

Does anyone have experience with this tool? What are your impressions? What is the learning curve like? Is there any other comparable extension for ArcGIS with similar functionality?


texts: Turning off specific attribute values in shapefile using QGIS

texts:
I realize from Splitting shapefile into separate files for each feature using QGIS or SAGA that I can split a shapefile by attribute, but I'd prefer to be able to just exclude certain data when I'm running various analyses on that data.
For example, I have a column "COLORS" in the attribute table, and I want to exclude only the cell value 'RED' within that column.
How do I do this?


texts: Projected layer misaligned with Mapserver

texts:
I have a couple or raster files in epsg:23031 reference that i want to reproject to other reference systems (epsg:4326 and google projection) using the 'on the fly reprojection' of Mapserver.

No problems using the default projection. But when i try to overlap my layer (in epsg:4326) with google earth there is a 100m misalignment between the two of them.

Any ideas?

Here is the map file:

MAP
 NAME topo
 IMAGECOLOR 125 125 125

 SIZE 600 400
 RESOLUTION 96
 RESOLUTION 300

 MAXSIZE 4096
 STATUS ON

 PROJECTION
  "init=epsg:23031"
 END

 OUTPUTFORMAT
  NAME png
  DRIVER "GD/PNG"
  MIMETYPE "image/png"
  IMAGEMODE RGB
  EXTENSION "png"
 END

 OUTPUTFORMAT
  NAME GEOTIFF_RGB
  DRIVER "GDAL/GTiff"
  MIMETYPE "image/tiff"
  IMAGEMODE RGB
  EXTENSION "tif"
 END

 EXTENT 443275 4647717 488074 4691799

 WEB
  IMAGEPATH "C:\OSGeo4W/tmp/ms_tmp/"
  IMAGEURL "/ms_tmp/"
  METADATA
  "wms_srs" "EPSG:4326 EPSG:900913"
  END
 END
 LAYER
  NAME orto
  PROJECTION
   "init=epsg:23031"
  END
  TYPE RASTER
  DUMP TRUE
  STATUS ON
  TILEINDEX "orto.shp"
  TILEITEM "LOCATION"
  EXTENT 443275 4647717 488074 4691799
 END
END

I'm using mapserver 5.6.4 that comes with osgeo4w tools


texts:
Working with QGIS3 (QGIS 3.2.3)
Is it possible to change what the default point symbol is when adding a new point layer to a new project?
Context: I have to work with massive point layers very often. The first action in QGIS I always do is change remove the stroke from marker symbol (Simple marker -> stroke style -> no pen). Stroke width 0 or stroke color = fully transparent would work too. I also always change the marker size. It would save me a lot of time if the default could be changed.


texts: Account for cell size differences in WGS 84 rasters

texts:
My data

I downloaded climate data from the WorldClim website, provided as rasters at a global scale using WGS 84 and various resolutions. I'm using the 15min resolution and I focus on Mediterranean countries (North Africa, Southern Europe and Western Middle-East).

My aim

I want to analyze climate differences between polygons but for several reasons I don't extract climate data for each polygon (using the extract function from the raster R package or QGIS Zonal Statistics). Instead, I rasterized those polygons and I compare raster cells from one polygon with raster cells from another polygon.

My problem

Raster cells don't have the same area, which means that small raster cells have more weight than they should, which could distort my analysis.

How could I take account of this difference in raster cells area?

Can I use another projection? I've read here that reprojecting rasters might not be a good solution. Also  I'm not sure that I can give different weights to those cells in the statistical analysis I use to compare cells.

Sorry if this topic had already been discussed, I did not find any solution on this website.

EDIT: I am seriously considering reprojecting my data using an equal-area projection like Lambert Azimuthal Equal-Area projection centered on my study area. Is this solution such a bad idea?


texts: Placing labels inside specified polygons/shapefiles in QGIS

texts:
Is it possible to place labels inside specified polygons in QGIS?

Case: I have a material which consist of water areas and various small islands. They are separated shapefiles. I would like to label the islands so that the label is placed always on the water area (-shapefile) next to the island. The islands are not equal in shape and size so the normal x and y based placement method is out of the question.


texts: Good algorithm to detect "wrong" features in a vectorial coverage

texts:
I am working with a bunch of shapefiles. Each shapefile consisting in rectangles with an integer attribute representing the strength of a signal. I serve these vector files as overlays on a base map, applying different colors based on the integer attribute (like a heat map, or a coverage map).

I have a problem with data: it contains "holes". The term "hole" here can be easily explained with an image:

That is a screenshot of a portion of my data loaded in QGIS. I've applied a simple style: just use orange for a range of values between X and Y. Otherwise, no color applied.

The white squares surrounded by blue circles are probably measurement errors. Each square represent 250m in the real world. It doesn't really make sense getting a measure of X in a point (i.e.: -55 dbM) and 250 meters beyond a value way lower (or higher), out of the represented range.

So, I'm basically looking for an algorithm to "detect" this pattern: One vector containing a integer value of -10000 (the "anomaly") but surrounded by vectors containing "normal" and coherent values. I want to avoid a simple dumb loop, checking all the vectors based on their position. I manage several GBs of data, and this can take forever.

Is there any known algorithm to detect this kind of problem in vectorial files?

UPDATE: There are some suggestions in comments. But I missed one important part: the space that surrounds the shapefile (the big white space around the shapefile) is considered "water", and a value of "no coverage" is allowed. So, I just need to get rid of the small holes (in blue in the image), but keep the surrounding space unchanged. That makes things more difficult...

Many thanks in advance


texts: Use Arcgis.com map package as a web map?

texts:
In Arcmap 10 there is a menu item to create a map package and upload to Arcgis Online (File > Create Map Package > upload to...). That works fairly smoothly, if a tad slow sometimes. However once placed in system it seems the only way to access the map is to download the map package. What we really want to do is create the map in Arcmap Desktop, upload to Arcgis Online, and then embed that map in a website.

If we take exactly the same data and create the map on the website, upload the feature classes individually, apply styling and symbology using the website tools, share, and so on. The resulting map can be opened in a) Arcgis.com map viewer, b) ArcExplorer, and c) Arcmap Desktop.

In the "My Content" page a map created from Arcmap Desktop is listed as a "Map Package", and a map created entirely using the website is a "Web Map" (what we want).

Am I missing something or is this really the way it's supposed to work? If so it's a terrific duplication of effort.


texts: gdal translate -co, which options are available to PNG?

texts:
I can create a PNG with:

gdal_translate -of PNG -outsize 10% 10% image.bsq image.png

I see there is a -co (creation options) parameter I can pass to gdal translate. Where can I find which options are available for PNG?


texts: Iterative merging of .shp per folder through multiple folders using ArcGIS Pro ModelBuilder

texts:
In ModelBuilder for ArcGIS Pro, I'm trying to create an iterative model for merging pointfiles according to their names. This is a pre-processing step to handle echo-soundings, spread over one or multiple files per lake, for the creation of lake-bathymetric maps.
The idea is to keep the point.shp files named according to "LakeName_n"(where LakeName = actual lake name and n=numerical order of files) in a folder named "LakeName" which is then kept in a main working directory. To use the model you'd only have to give the path to the main working directory for it to then go through each of the "LakeName" folders and merge (all) the LakeName_n.shp for each lake individually.
I have made a nested model that collects the values/files that are then to be used in the merge (img. 1 & 2).

I have gotten so far as to be able to either merge all the .shp files from all of the folders together or, by using wildcard, managing to do it for just one of the "LakeName" folders. But I can't seem to make this last atempt iterative. I.e. making a separate merged output for each lake.
I have tried using only the Name variable as input to both the Collect and Merge functions.
I have tried simply using %Name%* and also (by using parse path) %Workspace_name%* as wildcards, but the wildcard only seems to work if the "LakeName" is excplicitly spelled out for each individual lake.
I have allso tried adding "Iterate multivalues" within the parent modell without success (img. 3).

How can I do this?


texts: Dealing with thousands of features of vector layer in OpenLayers?

texts:
I am building a webgis portal too show a big number of polylines (road network) and points (villages, cities etc.).
All my data are stored in a POSTGIS and throuh geoserver and Openlayers I display them on the map.
For displaying I use two WFS vector layer (one for each geometry). I decided to use WFS instead of a WMS (WMSGetFeatureInfo) as this is what geoserver manual suggests here (regarding the GetFeatureInfo and GetFeature).

The problem comes when the number of features becomes relative big (around 10000 thousand features).
My question is how can I deal with such a big number of features.
I know one solution is BBOX but I am afraid this can slow things down by sending multiple requests.

Is there another strategy out there?


texts: JoinDataSource ArcGIS Javascript API, how to do it?

texts:
I'm looking to join a table to a polygon layer and I can't get the fields to transfer over. I'm sure there's something in my syntax that's messed up, but I can't figure out what it is. The goal of this process is to use a value from a field in the table to set the definition expression on the polygon layer (which I think I can do if the tables are joined?). It could be that I didn't create the workspaceId on the correct layer too... Or maybe I need to create a dynamic layer first? I'm utterly confused when it comes to dynamic layers.

This is what I currently have:

//Join habitat management attributes to the habitat management polygon layer
var joinDataSource = new JoinDataSource();
joinDataSource.joinType = 'left-inner-join';
joinDataSource.leftTableKey = "HabitatManagement.dbo.MgmtTracts.MgmtTractID";
joinDataSource.rightTableKey = "HabitatManagement.dbo.MgmtAttrb.MgmtTractID";

leftTableSource = new LayerMapSource({
      mapLayerId: 3
});

var rightTableSource = new LayerDataSource();
rightTableSource.dataSource = new TableDataSource({
      workspaceId: "hmgmtAttrb",
      dataSourceName: "HabitatManagement.dbo.MgmtAttrb"
});

joinDataSource.leftTableSource = leftTableSource;
joinDataSource.rightTableSource = rightTableSource;

var lyrDataSource = new LayerDataSource();
lyrDataSource.dataSource = joinDataSource;

//Add the habitat management tract feature layer
var hbMgmtTractFL = new FeatureLayer("https://stuff_here/rest/services/HabitatMonitoring/HabitatData/MapServer/3", {
       refreshInterval: 10,
       visible: false,
       source: lyrDataSource,
       outFields: ["*"]
});
hbMgmtTractFL.setMinScale(500000);
hbMgmtTractFL.setSelectionSymbol(selectionSym);

The JS runs without errors in the console, and I can see the source comes through:

In my server manager, I've set up the workspaceId on a connection with the HabitatManagement geodatabase instead of the HabitatMonitoring geodatabase. The map service I've set up the workspaceId on resides in the HabitatMonitoring folder. Would that cause an issue? The "MgmtTracts" and "MgmtAttrb" layer (the layers to be joined) reside in my HabitatManagement geodatabase which is why I created the workspaceId using that connection.


texts: Spatial 1 to 1 join by proximity in sf and R?

texts:
I'm trying to join a Simple Feature of POINT geometry type to another Simple Feature of POLYGON geometry type using st_join and nngeo::st_nn as the person asking this other question also did. The difference is that I want a 1:1 match from both sides (i.e. one point per polygon).

Above is a picture of how my data looks like when mapped. My code looks like:

matched<-st_join(buildings, addresses, join=nngeo::st_nn, maxdist=50,k=1)

Instead of a 1 to 1 match, I keep getting more than one point merged to the same polygon or more than one polygon merged to the same point depending on how order the arguments to st_join, with repeated points or polygons accordingly.

What I want this to do is to match each polygon to the closest point which also has that polygon as the closest to itself (i.e., to match polygons to points when they are mutually the closest to each other), and to randomise matching when more than one point or polygon meet those criteria.

Is there any way of doing that with sf or sp in r? I have thought about extracting the distance matrices and doing a two-sided matching between the two sets of features (like a Gale-
Shapley algorithm), but I want to give it a shot with what is already available before putting many hours into coding that.


texts: WMS of mixed cache from ArcGIS 10.2 not showing correctly in QGIS 2.14.9(Essen)

texts:
I'm trying to publish a cache in ArcGIS 10.2 to WMS and to show it on QGIS, but something there doesn't work right, the WMS is displayed very partially (only specific layers and not others, sometimes only labels).

I created a mixed cache in ArcGIS 10.2. The cache consists of several vector layers showing different elements (streets, cities, park, rivers, buildings, viewshed etc.) in different scales. When I Publish it to a service I enable the WMS option, copy the soap URL and after publishing, I go to my QGIS 2.14.9(Essen) > Add WMS Layer > New > Paste the URL> connect > choose upmost 'Layers' and add to the map. The cache that is added is very partial, showing very few features of only part of the layers, sometimes only labels.

Am I doing something wrong? Am I skipping a step? Or are mixed-caches not supported in QGIS? I have a suspicion that the difference between the scales is causing the problem (that QGIS doesn't know how to handle the different scales - that some layers are shown only in specific scales and not in others) but I can't prove it.

EDIT: I tried to use the fiddler, but he didn't show me any sensible URL (when adding the WMS to QGIS).
This is the soap url I thought I should use and it actually works when I give the url of a cache of an orthophoto. Withe vector map cache it still doesn't...


texts:
I have a set of coordinates that creates a polygon with holes in it.  I'm using qgis with mmqgis to extract.  How does MYSQL accept the WKT representation of a polygon with holes?  By holes I mean areas within the perimeter of the polygon that belong to that of another polygon, they're voids in the shape.

Is it with something like GeomFromText('MULTIPOLYGON((0,0 3,3 0,3 0,0), (1,1 1,2 2,2 2,1))')?

The first set of coordinates is the main polygon, the second set is the inner?


texts:
I created TIN to generate DTM using Mass points from Lidar data and river data that I made from stereoplotting. TIN results show there are several locations where the river elevation floats. So I have to improve the elevation data of the floating river so that it is at the correct elevation.

Is there an interactive tool when I edit the elevation of the wrong river then the TIN data changes automatically without me having to re-create the TIN data?

I have tried the TIN editor in ArcGIS but the tools only edit the TIN data.


texts: Nesting two iterators for cost distance calculation in ModelBuilder?

texts:
I am attempting to calculate cost distances iteratively for each point element in different datasets.

I created a model that does the job for a single dataset but when I try to batch it I'm getting the following error  "ERROR 000865: Input raster or feature source data: I_atlanta_OBJECTID does not exist".

This is, the output layer generated in the "feature select" submodel is not a valid input in the cost distance or in any other process I put after the iterator.

This error occurs only when I nested the model to run it in different datasets and not when I run the model in only one dataset.

Is there any consideration I have to take into account when nesting models?

I have read the documentation available with no luck so far.


texts: Trying to install SpatiaLite 4.3 with SQLite 3.8 on Windows 7

texts:
I'm new to the open-source GIS database world, and am trying to have SpatiaLite installed on my Windows 7 machine.

I downloaded SQLite 3.8.2 and unpacked the 3 files in a newly created folder c:\Program Files\SQLite. I added the folder to the PATH env variable, and managed to get that going.

To install SpatiaLite, I went to webpage https://www.gaia-gis.it/fossil/libspatialite/index, downloaded libspatialite-4.1.1.zip , which I believe is the latest SpatiaLite version for Windows. I unpacked this zip-file into the aforementioned SQLite folder, and then tried to run  .load 'SpatiaLite.dll' from the SQLite prompt. But this gives me an error.

Investigation shows that there is no such thing as 'SpatiaLite.dll', even though I thought I was just following the download and install instructions.

Where am I getting wrong?


texts: Intelligent hachure detection in ArcMap / QGIS

texts:
I have a DXF file containing residential and commercial buildings as well as the property lines. After importing it to ArcMap (10.2.1) I get something like this:

Now I want to:

Seperate the buildings from the property lines (which sometimes run through buildings)
Differ between residential (45 degree hachure) and commercial / other buildings (all other hachures)

To fulfill task 1 I already converted the polylines (each building originally consists of at least 4 lines + lines for hachure) to polygons and tried to aggregate them. But the result is always only one big polygon instead of many smaller ones.

Can anybody help me with this issue?

EDIT:

The result should look something like this (without the red patterns):

I just found another version of the DXF file which has two layers: one for the property lines and one for the buildings. So the property line problem is solved. But: All buildings are included in one layer and can't be differed by an attribute.


texts: What causes rainbow image artifact in Google Maps satellite imagery?

texts:
I discovered this "rainbow" over the Yellowstone River in Google Maps:

View the site in Google Maps: https://goo.gl/maps/y6xETmY859M2
You can see a distinct red, yellow, and blue streak hovering over this part of the river, just east of Lower Falls. Viewing the site in Google Earth Pro, it appears the image was taken in September 2015 (from a SPOT satellite I think? But I'm not certain).
I know these colors are not differences in geology because
(a) you don't see discoloration in images from earlier years and
(b) the discoloration is visible across the river itself.
First thing that came to mind was that the satellite captured an actual rainbow in the mist of the falls, but I don't think that's physically possible.
The closest I've come to describing it is resembling some type of "lens flare" effect.
What attribute of the satellite or camera might cause this to appear?


texts: Prefixing multiple values within QGIS using the Field Calculator?

texts:
I have a column in my my attribute table that contains a series of individually numbered locations. I wish to prefix all of these numbers with a letter, this letter will be the same for each of the values. E.g. 100, 200 and 300 would become A100, A200 and A300. I have a lot of data so simply going through and renaming would take a lot of time.

I have a feeling this would be achieved through the Field Calculator but unsure of which expression would be best for this?


texts: Making a clickable button in QGIS Attributes Form that is a URI from an Expression?

texts:
I would like to make a clickable button in my Attributes Form that is a URI that is dynamically generated using an Expression Calculator from fields in the data.

I have already created an action for the Attributes Table, though for my implementation, this is insufficient, hence asking about a button for the Attributes Form.

Expression Calculation:

trim(concat('bingmaps:?q=',"PARCELADDRESS",', ',"PARCELADDRESSCITY",', ', "PARCELADDRESSSTATE",', ', "PARCELADDRESSZIPCODE"))

Generated URI:
bingmaps:?q=123 Street Road, City City, State State, Zip Zip

This URI button would be clickable in my attributes form so that I can easily link the user to an appropriate map view in a third party tool.

Here's my current form, sans URI button to launch Microsoft Maps on the user's PC:

QGIS Version Information:
3.12.1-București - 121cc00ff0


texts: Qgis - need to draw exact lines (mtrs) with fixed angle

texts:
I use GIS 1.8.

I have the need to be able to draw exact lines and shapes on a bing map.
(This is to project a yet to be build greenhous project on a map to show what i would look like).

I have found the tool of Olivier Dalang (improved polygon capture) to be very usefull, since this tool enables me to draw lines while specifying the exact length and angle of the lines, and by doing so i build all walls of the greenhouse on the map).

Here is my problem.
The plugin enables me to enter the length in 'units', and drawing a line of 1000 units makes me come up with a line of about 622 meters.

I cannot find the exact way to change or find out how to draw in meters.

Any help highly appreciated!

(If there are other tools for drawing 90 degrees angled lines to a map I welcome these too!)


texts:
I'm trying to convert an osm.pbf file to a spatialite db but I can't get ogr2ogr to run.

I have qgis 2.14.2 and ubuntu 15.10.  I've reinstalled gdal-bin and I'm getting the following message from the terminal:

  ogr2ogr: error while loading shared libraries: libcharset.so.1: cannot
  open shared object file: No such file or directory

I've searched numerous places and I keep ending up back at instructions to install gdal.

Any suggestions?


texts: Statistics for geo-referenced floristic inventory

texts:
I have a dataset consisting of several thousand geo-tagged photographs of flowering plants (dozens of families, roughly 240 species).  The inventory was conducted monthly over three years along a trail in a forest.  I have identified everything, cleaned the database, and it is ready for analysis.  Since it is a point file, naturally I would like to use some of the latest statistical tools in ArcGIS/ArcGIS Pro to analyze the data and make some really awesome maps that show statistically significant results.  I would like to group the data based on the similarity (from taxonomy and traits - family, genus, growth habit, etc), frequency of occurrence (counts), and time (survey months) - but not necessarily in that order.  My general hypothesis is that there are distinct blooming seasons with distinctly different flowering "communities" (different families and species); though there would of course be some temporal and species-compositional overlap in these seasons.  The real question I want to get at is do these communities also differentiate themselves spatially?  Example:  are communities "A" and "B" (that are temporally clustered) also clustered in different areas on the map?

The problem I have run into is all of my data is text-based.  Grouping based on similarity and then a time series analysis would not be a problem in non-GIS platforms.  I know ArcGIS and Pro do grouping, but their help page says it does not do well with (or even use) nominal, categorical, or ordinal data.

I hope I have provided enough information to get good answers.  ANY thoughts on how to incorporate this analysis into ArcGIS or ArcGIS Pro?


texts: I'm a GIS newbie. I'm porting an app that uses PostGIS to SpatiaLite and wonder what the SpatiaLite equivalent of ST_MakeBox2D might be. Can someone point me to an example?


texts: Openlayers vector layer don't draw well on chrome (single tile)

texts:
I have singleTiled layer in OL (using chrome).
The layer don't draw completely unless I pan the map slowly to the south (here is example how the layer keep drawing itself each time I pan) :

This is how I create the layer :

var vectors = new OpenLayers.Layer.Vector('vectors', {
    isBaseLayer: false,
    singleTile: true,
    layerName: 'vectors',
    displayInLayerSwitcher: true,
    styleMap: roadsStyleMap
});

Any idea why it happens in chrome only ? (in FF I don't have this problem)


texts: Is geocoding with Nominatim available for all the countries

texts: I am trying to geocode an address using Nominatim for Malaysia but it cannot recognize it using his page: https://nominatim.openstreetmap.org/ . Is Nominatim available only for certain countries?


texts: Detecting and deleting alpha mark in QGIS

texts:
I got a couple of 4 band images, band 1-3 is RGB and band 4 is NIR.

In the properties it says nothing about an alpha mark, does that mean that none of them is marked as alpha channel or how do I detect an alpha channel (and if there's one, how can I remove it)?


texts: I got task, to plot polygon on map. After finish draw only marker inside a polygon appear on map. I'm done with leaflet.draw and Draw.Event.Created. What method I can use to check which marker inside area of polygon. A function same like a gif below.


texts: How to turn on the 'on-the-fly' functionality with PyQGIS?

texts:
Before QgsMapRenderer was deprecated if we wanted to activate OTF we used code like this:

canvas = QgsMapCanvas()
canvas.mapRenderer().setProjectionsEnabled(True)

Currently, QgsMapRenderer is divided into two classes QgsMapSettings and QgsMapRendererJob. I did some research, but with no success. So how we can achieve the same now?


texts: How to colour text background by category?

texts:
How can I change the colour of text labels by category?
EG, I have a text layer with text that is categorised into 6 different categories.  I need either the text itself, or its background to match the colour of polygons that are also in that category. I've hunted through here but can't find any answers that I can understand. It looks like the way forward involves some kind of rule based style, but I can't make head nor tail of the info I've found If I was using our old method, I would have text boxes with black text, a background that was a lighter shade of the colour used for the category, and an outline to the text box in a darker shade of that colour.
Is there a way to do something similar here? I've added some screenies. With the text all the same colour and the background all the same colour, it's difficult to see which label belongs to which object.


texts: Converting a CSV to DXF using ogr2ogr from the command line

texts:
I am trying to convert a .csv file to .dxf using ogr2org from the command line using the following command:

ogr2ogr -f DXF control.dxf opua.vrt

This is the .vrt file I am using:

<OGRVRTDataSource>
  <OGRVRTLayer name="opua">
    <SrcDataSource>opua.csv</SrcDataSource>
    <GeometryType>wkbPoint</GeometryType>
    <LayerSRS>WGS84</LayerSRS>
    <GeometryField encoding="PointFromColumns" x="X" y="Y" z="Z"/>
    <Field name="id" type="String"/>
    <Field name="description"/>
  </OGRVRTLayer>
</OGRVRTDataSource>

The .csv file contains these columns: point#,x,y,z,description. Using the suggestion found here (with the modified command above) produces a .dxf file with properly coordinated points points but no point number or description and two error messages:

ERROR 1: DXF layer does not support arbitrary field creation, field 'id' not created.
ERROR 1: DXF layer does not support arbitrary field creation, field 'description' not created.

I need to know what to do with these two errors to get my point numbers and descriptions to show up next to the points in the .dxf output.


texts: Area-weighted calculation on an intersection

texts:
I'm trying to learn postgis sql by performing some familiar operations from ArcGIS.  Here, I want to perform an area-weighted calculation of certain populations on the intersection of a set of parcels with census blocks.

I have parcels parcels, parcel attributes parcelinfo, census blocks census, and census attributes tgr39035sf1blk.

This calculation would entail

Calculating the area of each census block
Intersecting census blocks and parcels
Calculating the area of each polygon in the intersection of census blocks and parcels
Calculating the area weighted value of each attribute (attribute value * (area of  polygon / area of census block)

In order to test ST_Intersection, I tried the following, but phppgadmin craps out (no error message, a proxy error):

SELECT ST_Intersection(p.the_geom, c.the_geom) as inter from parcels p, census c;

I could sum the attributes from each block that intersects with a parcel, but that is, obviously, a seriously flawed approach.  E.g.,

SELECT   sum(DISTINCT t.pop2000), sum(DISTINCT t.males)
FROM     parcels p, census c, tgr39035sf1blk t
WHERE ST_Intersects(p.the_geom, c.the_geom) AND c.stfid = t.stfid;

I'm sorry if I'm missing any essential, additional information.  I'm a novice at this.


texts: How do add attributes from shapefile to list of lat/long-coordinates using QGIS?

texts:
In QGIS I have a csv-file which contains lat/long-coordinates and population data for each coordinate. I also have a shapefile containing certain regional districts, in my case US counties.

Now I want to add a column to my csv-file which in addition to the coordinates specifies the county of the coordinate point.

Since there are no common attributes to join on, how should I proceed using QGIS?

Thanks a bunch,

T


texts: So, I have an excel spreadsheet with information regarding people within a certain zip code. I basically need to make a map that shows the amount of people within each zip code either using dots or representative colors. How should I do this? The spreadsheet is set up with columns of name, email, and zip codes.


texts: Formatting SQL queries for use in Arcpy?

texts:
I have to make an SQL query on a column having not adequat name "type_borne_id", I couldn't make SQL queries using arcpy.Select_analysis neither selectbyattribute, but when I use select_analysis toolbox in arcmap it works with the same column.

I wonder if I should add some characters, or what to do, to query this column normally? I remarked that this pattern of names "name_column_id" is vulnerable to this error.

Example of the queries that didn't work:

arcpy.Select_analysis("table" , "table_copy"  , " type_borne_id = 23" )
arcpy.Select_analysis("table" , "table_copy"  , ' "type_borne_id" = 23' )

I have tested with other fields like "id" and it works normally.

I'm using Arcgis 10.1 & ArcSDE 10.1 on PostgreSQL.


texts:
I was working in ArcHydro 9 and ArcGIS 9.3. I have a street network and I was populating From_Node and To_Node columns in Attribute Table using ArcHydro 9 Attribute Tools tab. However only 1000 arc FIDs got populated in my ArcMap. Rest From_Node and To-node was populated null.

I have seen in Stack Exchange group someone referred to editing MaxRecordCount, but I can't locate it. Neither could I locate the path in Program Files\ArcGIS\server\user\cfg.
Can any one help?


texts: How to code to skip a shapefile, if it doesn't exist?

texts:
So if shapefile 'AAA010' does not exist, how can I skip the Merge and Clip and move on to the next part of the script? If possible I would like to keep the current script as it is, but by just adding else or if, Sorry still trying to learn the ins and outs.

import os, fnmatch, arcpy

arcpy.env.workspace = r'C:\Users\sysadmin\Desktop\DataMagement\Comp_3\SA_MGCP3.gdb'
fcMerged = 'MergedFeatureClass'
fcClipped = 'fcClipped'
fcToClip = r'C:\Users\sysadmin\Desktop\DataMagement\Country.shp'

filesMatched = list()

for root, subfolders, files in os.walk (r'C:\Users\sysadmin\Desktop\DataMagement\MGCP_Raw'):

    for f in fnmatch.filer(files, 'AAA010.shp'):
        filesMatched.append(os.path.join(root,f))

    if filesMatched:
        print ' files found - merging'
        arcpy.Merge_management (filesMatched, fcMerged)
        arcpy.Clip_analysis (fcMerged, fcToClip,fcClipped)

    else:
         print 'nothing found - skipping'


texts: SAGA tool 'Relative Heights and Slope Positions' - What do results tell me?

texts:
I am currently trying to get in morphometry, using the SRTM 30m DEM. I was recommended to test the SAGA mudule 'Relative Heights and Slope Positions' (http://www.saga-gis.org/saga_module_doc/2.1.3/ta_morphometry_14.html) in QGIS.
It runs smoothly and I have first results, but I am a bit lost in understanding what results tell me. Does anybody of you have expert knowledge to explain what the following mean?

Slope height (values ~0-700)
Valley depth (values ~0-600)
Normalized height (values ~0-1)
Standardized height (values ~0-3000)
Mid slope position (values ~-0,00001-1)

Input options w, t and e are also not clear to me.

If anybody of you has an idea or an idea where to find information, I would be very happy.


texts:
I have a problem in Arcmap 9.3 When we are using 3D data in ArcMap to dissolve it the elevations are lost.

Is there someway to keep the elevation during the dissolve method?