texts:
I have implemented a GeoTIFF image layer usinga  WCS service in GeoServer 2.2.4. I have searched on the net for operations supported by WCS... it tells me about three operations:

GetCapabilities, DescribeCoverage and GetCoverage

I have used these operations using HTTP request and all are working perfectly.

My questions are:

is there any other way to use these operations other than HTTP request (if yes, please explain with some example)
does WCS support any other feature/facility/operation other than the above mentioned operations.

Is it possible to apply stretch operations (standard deviation/histogram etc) and other operations on GeoTIFF image in GeoServer using WCS?


texts: GeoServer option "followline" results in no labels on a river layer

texts:
I'm using GeoServer 2.10.1 with a river layer and I'm trying to add labels which follow the lines. (The river layer is sourced from a PostGreSQL table if that makes a difference.)

Using the line with default label SLD, the rivers are labelled:

But when I use the label following line example (which includes the <VendorOption name="followLine">true</VendorOption> parameter), the labels disappear:

   <FeatureTypeStyle>
     <Rule>
       <LineSymbolizer>
         <Stroke>
           <CssParameter name="stroke">#FF0000</CssParameter>
         </Stroke>
       </LineSymbolizer>
       <TextSymbolizer>
         <Label>
           <ogc:PropertyName>name</ogc:PropertyName>
         </Label>
         <LabelPlacement>
           <LinePlacement />
         </LabelPlacement>
         <Fill>
           <CssParameter name="fill">#000000</CssParameter>
         </Fill>
         <VendorOption name="followLine">true</VendorOption>
       </TextSymbolizer>
     </Rule>
   </FeatureTypeStyle>

Is this due to the complexity of the rivers, or is there another issue occurring here?


texts:
I am using a Feature Class Iterator in a Submodel ending with Collect Features:

The main model looks like this:

When I want to run it I only get this error:

I have no clue why he says the value is empty.
Does anyone know it?


texts: Making Raster of Northings (for use in Maxent) in ArcMap?

texts:
I'm working on building Maxent species distribution models for rockfish along the Northern California Coast. To do this I need rasters of predictor variables, and I want to capture latitudinal variation in habitat.

Basically I want to duplicate the Northing coordinate of the raster into the cell value. Seems simple, but it's been driving me nuts.

I'm running ArcMap 10.2 with the spatial analyst package.


texts: Google Earth Engine chart error: No features contain non-null values of "system:time_start"

texts:
I'm trying to plot a chart from a collection I have created. While the data from this collection maps fine, I can't seem to produce a chart for it.
var sstMonthly = temporalCollection(sst, ee.Date('2013-01-01'), 12, 1, 'month');
print('sstMonthly', sstMonthly)

print(ui.Chart.image.series(sstMonthly, BIOT, ee.Reducer.mean(), 500));

Which gives this error;
Error generating chart: No features contain non-null values of "system:time_start".

From what I've read it seems that the error is probably because it can't find the property. But nowhere states how to find that property in the image collection or what the correct syntax to add it into the code for charts is. Any help?


texts: Leaflet set the ImageOverlay switched off at the initial stage

texts:
I would like to set my overlayimage off when reload the map.

According to the link:

https://stackoverflow.com/questions/46043995/set-leaflet-overlay-off-in-the-layer-control

I tried sth like this:

   var overlayMaps = (null,{
    "Geographic names": positronLabels,
    "City fibre": cityfibre
  });

but there is no reaction on my map

In turn this example:

https://stackoverflow.com/questions/45333282/remove-leaflet-overlaymaps-layer-programmatically

Esri leaflet layers turned off by default in map layer control box

shows how to remove some layer instead of switch it off at the initial stage.

When I put null into my code, my overlay layers are gone:

    L.control.layers(baseMaps, null, overlayMaps).addTo(map);

I couldn't find a good solution on the Leaflet.js reference:

https://leafletjs.com/examples/layers-control/

Is there any option to make them switched off at the initial stages?


texts: Extract raster data along shapefile in QGIS

texts: I'm trying to extract raster values along multiple shapefiles in QGIS for further analysis. I've successfully extracted values for points, so I either need to convert my vector shapefile to points (with the same spacing as my raster pixel resolution), or I need a tool to extract values along the vector line. I've tried Zonal Statistics, but I need to extract each pixel value, not the sum/count/mean etc.


texts: Get values from select by location to append to table

texts:
I'm trying to calculate the road frontage of polygons by converting polygons to lines and splitting in their vertices, then using select by location to find lines that are within 7.5 meters of the road. I then want to retrieve the road label and append it to the lines table. This is what I have, which find the lines correctly and writes the road information to the lines table, but it doesn't get the correct road info. It doesn't seem to find the road that was used in the select by location tool. Can anyone help?

arcpy.PolygonToLine_management("BLPUs", "BLPUs_PolyLine")
arcpy.SplitLine_management("BLPUs_PolyLine", "BLPUsPolyLineSplit")
arcpy.MakeFeatureLayer_management("BLPUsPolyLineSplit", "BLPU_PolyLineSplit")
arcpy.SelectLayerByAttribute_management("BLPU_PolyLineSplit", "NEW_SELECTION", '"LEFT_FID" <> -1 ')
arcpy.DeleteRows_management("BLPU_PolyLineSplit")
in_layer = "BLPU_PolyLineSplit"
out_layer = "BLPUs_PolyLineSplit"
arcpy.CopyFeatures_management(in_layer, out_layer)
arcpy.AddField_management("BLPUs_PolyLineSplit", "RoadLabel", "TEXT")
arcpy.AddField_management("BLPUs_PolyLineSplit", "DescTerm", "TEXT")
arcpy.AddField_management("BLPUs_PolyLineSplit", "Nature", "TEXT")
arcpy.MakeFeatureLayer_management("BLPUs_PolyLineSplit", "BLPU_PolyLineSplitt")
with arcpy.da.UpdateCursor("BLPU_PolyLineSplitt", ["RoadLabel", "DescTerm", "Nature"]) as cursor:
    field1 = "RoadLabel"
    field2 = "DescTerm"
    field3 = "Nature"
    for row in cursor:
        arcpy.SelectLayerByLocation_management("BLPU_PolyLineSplitt", "WITHIN_A_DISTANCE", "C:\Users\Ant\Documents\ArcGIS\OSMaster.gdb\ITN_RoadLinks", "7.5 Meters")
        TABLE = "ITN_RoadLinks_Copy"
        searchRows = arcpy.da.SearchCursor(TABLE,["RoadLabel", "DescTerm", "Nature"])
        for searchRow in searchRows:
            searchRow1 = str(searchRow[0])
            searchRow2 = str(searchRow[1])
            searchRow3 = str(searchRow[2])
        row[0] = searchRow1
        row[1] = searchRow2
        row[2] = searchRow3
        cursor.updateRow(row)
arcpy.CopyFeatures_management("BLPU_PolyLineSplitt", "BLPUs_FrontageLines")


texts: How to stop writeOGR from abbreviating Field Names when using "ESRI Shapefile" driver

texts:
I'm currently using the following script to add some attribute data from a table to lots of individual shapefiles:

library(rgdal)
specieslist <- read.csv("SpeciesList1.txt", header=F)
attdata <- read.table("TestAtt.csv", sep = ",", header=T)
for (n in 1:dim(specieslist)[1])
{
speciesname <- specieslist[n,1]
shp <- readOGR("Mesoamerica_modified_polygons", speciesname)
shp$ENGL_NAME<-attdata[n,2]
writeOGR(shp, "PolygonsV2", speciesname, driver="ESRI Shapefile")
}

I get the following warning(s) at the end:

1: In writeOGR(shp, "PolygonsV2", speciesname, driver = "ESRI Shapefile") :
Field names abbreviated for ESRI Shapefile driver

When viewing the attribute table of the shapefiles after this process, the Field name has been shortened to 'ENGL_', but I want it to stay as 'ENGL_NAME'. Is there a way to turn this abbreviating off?

Any help much appreciated.


texts: Extracting desired values from KML file using QGIS?

texts:
I have a kml file of roads that I would like to convert into a shapefile using QGIS, I understand that this should be a straightforward task however I am running into the issue of extracting necessary data from each feature in the kml file. Essentially I want all of the information from each kml feature (see google earth screenshot).

I want each of those attributes (FID,  AMA_NAME, Shape_Length, etc...) to be displayed in the attribute table in QGIS as separate values. However, when I bring this kml file into QGIS and view the attribute table all of this information is grouped within the description field for each feature (see below QGIS screenshot).

Is there a way to get the information for each feature from the description field into an attribute table or other form that can be queried in QGIS?


texts: How do you represent parallel/meridian bounding boxes in a geography column?

texts:
This PostGIS one-off query returns true:

psql=> select ST_Covers(
    ST_SetSRID(ST_MakeBox2D(ST_Point(-180, -30), ST_Point(20, 90)), 4326),
    ST_SetSRID(ST_Point(-77.4874420166016, 39.043758392334), 4326));
st_covers
-----------
t

But if you cast both geometries to geographies, it returns false:

psql=> select ST_Covers(
    ST_SetSRID(ST_MakeBox2D(ST_Point(-180, -30), ST_Point(20, 90)), 4326)::geography,
    ST_SetSRID(ST_Point(-77.4874420166016, 39.043758392334), 4326)::geography);
st_covers
-----------
f

I understand from this old question and this other old question that this is because geography uses great circles to connect points in a polygon, not parallels and meridians.

The question is, how do you convert a parallels-and-meridians bounding rectangle to a geography object that covers exactly the same region of the earth?  In the real thing, these rectangles go into a table which is then joined against another table with a long list of points, and that table has to use geography for other reasons, so I need these parallel-and-meridian rectangles as geographies so the indices work properly.


texts: I have downloaded MODIS images for a specified area.. Its my first time that I am working on UHI (Urban Heat Island).. I need to know the methods to work on it.. I am using ENVI for this.


texts: Mapping with time frame symbology in QGIS?

texts:
I am using QGIS 3.8, and I am making a map of archaeological sites. My table looks something like this:
NAME - Coordinates X - Coordinates Y - Site_type - Data_type - Century_begin - Century_end

where century_begin is the first century since my data is relevant and century end is the last (e.g. Archaeobotanical analyses only relevant from century 5th to 7th).

I want also to produce now a temporal map, in order to show the sites relevant to each century.
E.g.
5th c. (Display sites)
6th c. (Display sites).

Do you know how I can obtain this sort of map?

I tried to use the categorized symbols, but the problem here is that each site has a coverage larger than 1 century (it's always something like 5th-6th).

Solved with help from @csk, but with just one problem. There is a site always showing. This is how my table looks like (It is not editable because I imported it as delimited text - I will Google how to transform it into an editable object).


texts: Set Data Source Option in ArcGIS Pro - Similar to set data source for mxd

texts:
Where can I find following option to change data source of all layers in map in ArcGIS Pro?
This is how I access data source in ArcMap


texts:
Given that AGS does not allow enabling security on individual services, what is the best approach for hosting a mix of both secure and open services?

Update: we have abandoned AGS-based security because of its many flaws and capricious inconsistencies. Beware!


texts: Why my tooltip cannot shown when I hover over layer in leaflet.js?

texts:
I want to use tooltip on my polyline to show its features and I use this code:

var boroughs = L.geoJson(null, {
  style: function (feature) {
    return {
      color: "black",
      fill: false,
      opacity: 1,
    };
  },
  onEachFeature: function(feature, layer) {
    layer.bindTooltip(feature.properties.KECAMATAN, {
      sticky: true
    });
    boroughSearch.push({
        name: layer.feature.properties.KECAMATAN,
        source: "Boroughs",
        id: layer.feature.properties.KODEKEC,
        bounds: layer.getBounds()
      });
    }
});
$.getJSON("data/kec.geojson", function (data) {
  boroughs.addData(data);
});

and I got a message in the console that said:

  TypeError: layer.bindTooltip is not a function

I don't know what's wrong with it, and I followed step by leaflet docs.


texts: openlayers wms loaderror event get response text of

texts:
1.How do I get the http response text & error code from OpenLayers 2.12 WMS loaderror event ?

2.The WMS request returns returns http error code 500 with some response text, however, as per image below, I can't get those details.

3.Workaround: I send a HEAD request in the 'loaderror' handler

                    intramapsConfig = {
                        name : "MapImage WMS",
                        options : { "ratio": 1, "singleTile": true,
                                    "removeBackBufferDelay": 500 ,
                                    "transitionEffect" : "resize",
                                    "tileOptions": {
                                    "eventListeners": {
                                            'loaderror': function(evt) {
                                                console.log('loaderror');
                                                debugger;
                                                $.ajax({
                                                    type: "HEAD",
                                                    async: false,
                                                    url: evt.object.url,
                                                    success: function(message,text,response){
                                                        debugger;
                                                    },
                                                    error: function(error) {
                                                        debugger;
                                                    }
                                                });

                                            }
                                        }
                                }
                        }

4.The image below is the event object sent by onImageError function.


texts:
I want to make dynamic filter on Leaflet, but I got a little problem.
so first this my dropdown:

<select id="selc">
<option value="baik">baik</option>
<option value="sedang">sedang</option>
<option value="berat">berat</option>
</select>

and this my filter:

filter: function(feature, layer) {
    <?php
    $asa = $this->db->query("SELECT du_kode,du_kondisi FROM data_umum inner join data_umum_head on data_umum.duh_id=data_umum_head.duh_id where data_umum_head.maplink_id = '$ml->maplink_id'")->result();
    foreach($asa as $dew){ ?>
        woe = <?= $dew->du_kode?>;
        if(feature.properties.ID==woe){
        wawa = '<?= $dew->du_kondisi?>';

        $('#selc').on('change', function() {
        console.log($('#selc').val(),feature.properties.ID,wawa);
        if(wawa == $('#selc').val()){
            return (feature.properties.ID)
            }
        });
      }
<?php } ?>
},

My problem is why my on.change is not working to return (feature.properties.ID)

but if I do manually and change the code like this

filter: function(feature, layer) {
<?php
    $asa = $this->db->query("
    SELECT du_kode,du_kondisi FROM data_umum
    inner join data_umum_head on
    data_umum.duh_id=data_umum_head.duh_id where
    data_umum_head.maplink_id = '$ml->maplink_id'")->result();
                        foreach($asa as $dew){ ?>
                            woe = '<?= $dew->du_kode?>';
                            if(feature.properties.ID==woe){
                                wawa = '<?= $dew->du_kondisi?>';
                                /*$('#selc').on('change', function() {
                console.log($('#selc').val(),feature.properties.ID,wawa);
                                    if(wawa == $('#selc').val()){
                                        return (feature.properties.ID)
                                    }
                                });
                                */
                                console.log(wawa);
                                if(wawa == 'berat'){
                                    return (feature.properties.ID)
                                }

                            }

                        <?php } ?>
                    },

and my results be like this

This is the result I expected, but how can I use it with onchange?


texts:
I'm working with qgis 2.8 (windows).

I am able to generate maps using qgis atlas. But i have a problem with the scale bar.

When i generate atlas with "fix scale" setting (1:35.000 for instance), the scale length changes for each "slide", each exported map.

Could you tell me why?

I hoped to solve the problem using "predefined method". In vain.

Below are two screenshots showing my problem :


texts: 3D topographic progress compensation by the least squares method

texts:
I'm looking for an explanation of the least squares method used in the case of a correction of 3D point network.

We have reference points with known coordinates XYZ. We calculate intermediate points based on chaining from the starting point (reference point) until we reach the end point (other reference point). We then calculate the difference between the measured coordinates and known coordinates of the end point, so the difference is the cumulative error. We now have to split the cumulative error between all measured points by a least square method.

So please does someone have a document explaining the mathematical formulas, what to do, and sample application of this method?


texts: Auto-incrementing field in feature class using ArcGIS Desktop?

texts:
So I'm relatively new to ArcGIS and ArcCatalog.  I have created a field name of PROJECTID in a feature class in a Personal Geodatabase.  The data type is TEXT.  This will be a ten digit ID.  The first five characters will be a acronym.  The sixth character will be a digit that represents the type of feature.  It will be a number.  The last four characters will be numbers that should be incremented starting at 0001.

Example:
    WATER-1-0001

Screen shot of the field:

So my question is how can I do this in ArcCatalog.  I've seen some Python code that will autoincrement like the following:

rec=[OBJECTID]
def autoIncrement(a):
 global rec
 return "water" + str(a) + "-" +  format(rec, '04d')

Again I'm a newbie when it comes to ArcGIS .
Now I am getting an error NameError name 'OBJECTID' is not defined.


texts: how to use ArcGIS cloud hosted map with custom data from mysql

texts:
I have map hosted on ArcGIS.com which I can use and embed in my webpage. I have few tables with latitude-longitude in mysql. is there anyway I can add it as marker overlay layer on map taken from ArcGIS.com ?

Note: My server does not support spatial tools like PostgreSQL and Geoserver. SO I am trying to find an alternative. any help and guidance will be appreciated.


texts: Creating linestring geometry from jsonb column

texts:
I am trying to create a line from an array of x,y coordinates inside of a jsonb field. Here is the format of the field:

I created a linestring geometry column, but having difficulty with the syntax. Do I need to make points from the x,y first?  The j.line part is where I think i need to make points?

UPDATE waze.wazejams SET geom = ST_SetSRID(ST_MakeLine(**j.line**),4326)

I'm very green with PostGIS so go easy :)

UPDATE: I'm thinking I need to create a multipoint geometry to house the multiple vertexes along the line.


texts: How do you preform multi spectral image classification on non Landsat images in QGIS

texts:
My aim is to classify vegetation types, using areal/ satellite imagery, in a small (4km^2) area of the North York Moors (UK). This area contains two small rivers and their catchment areas; one of the rivers is in an area of planted forestry and the other open heather moorland. I plan is to classify the image to better discriminate between the vegetation types to teach students about the effect vegetation has on rivers.

I have an area in mind and have done this kind of classification exercise before, however only when following tutorials and using Landsat data in ArcGIS. I now need help locating multi spectral imagery of high enough quality for my 4km area (so not Landsat) and help to classify it using QGIS.

My question is therefore two fold:

a) where can I purchase/ acquire such high quality data?

b) how do I preform classification on a multi spectral image which is not Landsat in QGIS (or is it even possible)?


texts: Why am I getting 'Script error' - Google Earth Engine?

texts:
I want monthly average radiance composite images using nighttime data from the Visible Infrared Imaging Radiometer Suite (VIIRS) Day/Night Band (DNB) for Mopti region. I'm able to export the result to drive, however, after some time has elapsed I get this message 'Script error' in the print console. How can I rectify this?
var dataset = ee.ImageCollection('NOAA/VIIRS/DNB/MONTHLY_V1/VCMCFG')
                  .filter(ee.Filter.date('2019-11-01', '2019-12-31'));
var nighttime = dataset.select('avg_rad');
var nighttimeVis = {min: 0.0, max: 60.0};
Map.centerObject(mopti, 8);
Map.addLayer(nighttime, nighttimeVis, 'Nighttime');

print(nighttime.first());

Export.image.toDrive({
  image: nighttime.first().clip(mopti),
  description:'Mopti_NightLight_November',
  folder: 'VIIRS_Night_Light_2019_Monthly_Data',
  scale: 400,
  region: mopti
});


texts:
I have been struggling to install PostGIS on my Redhat Linux based CentOS 7 (x64) machine and so far I am still running into problems. I used this step by step tutorial developed by Regina at PostgreSQL Online to install PostgreSQL 9.5 and PostGIS and everything seems to be fine until I started installing PostGIS. I am running into dependencies problems! After installing PostgreSQL 9.5, when I typed this command:

# sudo yum install postgis2_95 postgis2_95-client

I am getting these dependency errors:

Loaded plugins: fastestmirror, langpacks
Loading mirror speeds from cached hostfile
* base: ftp.crc.dk
* epel: ftp.crc.dk
* extras: ftp.crc.dk
* updates: ftp.crc.dk
Resolving Dependencies
--> Running transaction check
---> Package geos-devel.x86_64 0:3.5.0-1.rhel7 will be installed
--> Processing Dependency: geos = 3.5.0-1.rhel7 for package: geos-devel-3.5.0-1.rhel7.x86_64
--> Running transaction check
---> Package geos.x86_64 0:3.4.2-2.el7 will be updated
--> Processing Dependency: libgeos-3.4.2.so()(64bit) for package: grass-6.4.4-7.el7.x86_64
--> Processing Dependency: libgeos-3.4.2.so()(64bit) for package: grass-libs-6.4.4-7.el7.x86_64
---> Package geos.x86_64 0:3.5.0-1.rhel7 will be an update
--> Finished Dependency Resolution

Error:
Package: grass-libs-6.4.4-7.el7.x86_64 (@epel)
    Requires: libgeos-3.4.2.so()(64bit)
Removing: geos-3.4.2-2.el7.x86_64 (@epel)
libgeos-3.4.2.so()(64bit)
Updated By: geos-3.5.0-1.rhel7.x86_64 (pgdg95)
    ~libgeos-3.5.0.so()(64bit)
Available: geos-3.4.2-1.rhel7.x86_64 (pgdg94)
libgeos-3.4.2.so()(64bit)
Available: geos-3.2.2-2.el6.x86_64 (elgis)
~libgeos-3.2.2.so()(64bit)
Available: geos-3.3.0-1.el6.x86_64 (elgis)
~libgeos-3.3.0.so()(64bit)
Available: geos-3.3.1-2.el6.x86_64 (elgis)
~libgeos-3.3.1.so()(64bit)
Available: geos-3.3.6-1_0.el6.x86_64 (elgis)
~libgeos-3.3.6.so()(64bit)
Available: geos-3.3.8-2.el6.x86_64 (elgis)
~libgeos-3.3.8.so()(64bit)

Error: Package: grass-6.4.4-7.el7.x86_64 (@epel)
Requires: libgeos-3.4.2.so()(64bit)
Removing: geos-3.4.2-2.el7.x86_64 (@epel)
libgeos-3.4.2.so()(64bit)
Updated By: geos-3.5.0-1.rhel7.x86_64 (pgdg95)
~libgeos-3.5.0.so()(64bit)
Available: geos-3.4.2-1.rhel7.x86_64 (pgdg94)
libgeos-3.4.2.so()(64bit)
Available: geos-3.2.2-2.el6.x86_64 (elgis)
~libgeos-3.2.2.so()(64bit)
Available: geos-3.3.0-1.el6.x86_64 (elgis)
~libgeos-3.3.0.so()(64bit)
Available: geos-3.3.1-2.el6.x86_64 (elgis)
~libgeos-3.3.1.so()(64bit)
Available: geos-3.3.6-1_0.el6.x86_64 (elgis)
~libgeos-3.3.6.so()(64bit)
Available: geos-3.3.8-2.el6.x86_64 (elgis)
~libgeos-3.3.8.so()(64bit)
You could try --skip-broken to work around the problem.

Could anyone help me out to resolve these dependency issues or point out whether I am doing something wrong?


texts: Combine multiple records into one based on ID

texts: I have a table that has the building width for each side and I want to combine them together in GIS. I have attached an example in which I have 4 rows for one building, and I want only one row per building as shown on the picture. Is there a method in GIS?


texts:
I stumbled upon your question regarding the mosaic feature in qgis, the thing is I'm also using Quantum and this is my first encounter with GIS this semester, so if you can please help me with a technical question.
For the watershed analysis I know the steps as a concept, like how I should first get the flow direction, then flow accumulation then the stream network etc.. but technically I don't know how to locate these functions in Quantum, if you have any knowledge about that matter can you please share it with me?


texts: I need to store raster data on a cluster with CDH 6.2.0 either on HBase 2.1.0 or on Hadoop HDFS 3.0.0. I also need to connect rasters to GeoServer (2.16) to publish map layers. I tried to integrate the elements with GeoMesa (2.4.0) but the installation wasn't successful. According to the documentation only HBase 1.3.x and 1.4.x are supported and I didn't even find any reference to Hadoop 3.0.0. Is this solution compatible with my cluster?


texts: How to use "Merge vector layers" in QGIS using Geopackages as output with FIDs duplicated

texts:
I'm using QGIS 3.2.3 in Windows and I created a tessellation with polygons (Voronoi tiling / Thiessen polygons) with the built-in tool "Voronoi Diagrams", (Call this layer  LAYER1). Afterwards, I used LAYER1 to subdivide another layer stored as a Geopackage with many polygons (call this layer LAYER2) using the built-in tool "Clip" and the option "Iterate over this layer" to create a separate layer for every feature in LAYER2 with the subdivision of LAYER1 and set as output for each one a new Geopackage (call this layers OUTPUTS).

Now I want to merge all this OUTPUTS in a new layer, so I get a tesellation of polygons in LAYER2, for that I tried to use the tool "Merge Vector layers" that I suppose replace and extends the old tool "Merge Shapefiles", but as I'm using now geopackage this process fails (If I set the output format as Shapefile it runs successfully, because that format doesn't use a unique ID).

I suppose that, this problem is caused because all features in all OUTPUTS must to have different FID to be merged, but that's no my case because the FID in the outputs are based on FIDs of LAYER1 so the same polygon can divide many others in LAYER2.

Is there a solution to this problem using QGIS built-in tools? The easy way could be  just use as output format Shapefile. But I like to maintain the consistency of my project and avoid to rename some fields in my layers, or edit manually the fids in the outputs.


texts: Open-source software for managing terabytes of different formats of geospatial data?

texts:
I have a problem that I am sure many GIS Teams or Departments have which is how to find, keep track of and manage multiple different datasets from multiple sources. I have read the responses from a similar question asked two years ago

Managing large amounts of geospatial data?

I was wondering if there are any other software solutions, I know about VoyagerGIS, which looks good, but it is linked with ESRI and wondered if there was a similar Open Source version.

To give you an overview of the current situation:

Terabytes of Imagery which changes every 6 weeks, I would like to keep historic
1 Terabyte of MapInfo TABs
ESRI Shapefiles again changes every 6 weeks
100s Gbs of CSV, GeoJSON, KML, SQLite, Excel spreadsheets, JPEGs, TIF files (some georeferenced, some with world files), these come in from multiple different sources.

I guess my perfect solution would be a GUI interface to allow me to update metadata, search based on filename or even attribute, a visual image (like seen in VoyagerGIS), it show metadata, a list of attributes. However, I would like to keep the original format of the data rather than loading it all into PostGIS.


texts: Issue with using TopoJSON server API to return TopoJSON object with properties

texts:
I am creating a node application and using TopoJSON server API reference (see the link below) to convert a GeoJSON object to a TopoJSON object.

TopoJSON server API reference

I used these three lines of code below which returns the TopoJSON topology but without returning the properties from the GeoJSON object.

var collection = {type: "FeatureCollection", features: […]}; // GeoJSON
var topology = topojson.topology({collection: collection}); // convert to TopoJSON
console.log(topology.objects.collection); // inspect TopoJSON

However, I also want to include the properties of each topology to the final TopoJSON object, and according to the API, I need to use this function below copy all properties from the GeoJSON object:

function propertyTransform(properties, key, value) {
  properties[key] = value;
  return true;
}

But my issue is I don't know how to use this propertyTransform function as an option inside of the topojson.topology() function in order to return the TopoJSON object with properties.

Does anyone know how to do it or provide me some example code?

Thanks a lot!

Vanni


texts:
I am working on a forked OSM/leaflet/overpass API application which has its own max bounds declared

var mapBounds = { south: 50.8025, west: 0.3724, north: 50.8785, east: 0.5290 };
var LBounds = L.latLngBounds([mapBounds.south, mapBounds.west], [mapBounds.north, mapBounds.east]);

followed by

var map = new L.map('map', {
    maxBounds: LBounds.pad(0.5),
    maxBoundsViscosity: 1.0,
    minZoom: 10,
    maxZoom: 20,
...

The query to Overpass API is done via Nominatim

geocoder: L.Control.Geocoder.nominatim({
    geocodingQueryParams: {
        bounded: 1,
        viewbox: [mapBounds.west, mapBounds.south, mapBounds.east, mapBounds.north].join(',')
    }
})

I have set the mapBounds to allow world view

var mapBounds = { south: -90, west: -180, north: 90, east: 180 };

and then do something like

var mapBounds='bbox=' + map.getBounds().toBBoxString();

in order to restrict Overpass API requests to the current viewport only
I couldn't make it work yet.


texts:
I would like to retrieve the list of the different WRB soil classes by country (https://data.isric.org/geonetwork/srv/fre/catalog.search#/metadata/c4dc161c-d62d-11ea-a1a3-292680b15169).
I would like to get this data in order to know the diversity of soils by country.
Can you give me some line of thought ?
I looked for if a file (CSV format) exist but I did not find anything. Perhaps an R package like SoilGrids should be used in order to retrieve the list of the different classes by country.


texts: I have a data with an unknown coordinate system. The coordinates has 10 digits left of the decimal. Its unit is equivalent to about a cm. The data is for Tehran(Iran), boundaries of Tehran are about 35.54 to 35.84 N and 51.19 to 51.84 E in decimal degrees. The boundaries in this coordinate system are 1499559806.00 to 1496487278.00 N and 604413102.00 to 609510913.00 E.


texts:
I have a problem passing variables through infoTemplate. Its like there's something funny happening with the scope. Certain variables work but when I get too crazy with them my map doesn't load.

Example1 works:

template = new InfoTemplate({
    title: welltitle,
    content: "ID: "+coords[i][0]+"
"+coords[i][0]+"
"+datatext
});

Example2 does not:

template = new InfoTemplate({
    title: welltitle,
    content: "ID: "+coords[i+1][0]+"
"+coords[i+1][0]+"
"+datatext
});

EDIT coords is a multidimensional array of values

EDIT NEW SYMPTOM! map crashes in firefox but not Chrome...this may be just a cross-browser issue


texts: Efficient WMS or TMS Loop for Radar data?

texts:
I'm looking for an efficient way to loop radar data from a Tilecache or WMS Server.  Right now, I make an array of layers and increment through them, but this seems inefficient since I end up with a lot of layers for what should be just one layer. I saw this post: https://gis.stackexchange.com/a/20836/7387, however, I'm not sure how to go about implementing what is suggested...

          {Destroy Layers} |10|11...19|20| {Start Loading Layers}
|---------------------------------------------------------------------------|
                           ↑
        Timespan           Current Possition

Can you do this on a single layer? A simple example to follow would be very helpful. Thanks!


texts:
Looking to create a mobile iOS app to collect consumer feedback on 10 separate locations using a 5-star rating system.  There would be a total of 4 questions, each of which would require a separate rating.

I currently have my app set up with a map displaying 10 point features and geofences around each, so that when a user enters the geofenced area they receive a push-notification asking them to complete a survey.  From this point on is where I need some assistance.

I am not sure how to go about creating a 5-star survey that links back to my feature layer so that when results are collected they can be displayed back to the user via callout box when they select a feature on the map.  Also looking for recommendation on where to host my data so that it is best suited for this (i.e. Parse, CartoDB, etc).


texts:
Currently, the pgrouting 1. 5 version has a shortest path function which actually traces a path from nearest vertex of source to nearest vertex of destination but. As, Dijkstra's algorithm works on the principle of vertex to vertex routing. But, if the location (point geometry) lies at somewhere in the mid of the road geometry and I want the path to be traced until that location exactly, is there a special function or way to do the same in pgrouting?
In short, currently node to node routing happens. Whereas, I am looking for point to point routing.

Please refer: http://www.mapgears.com/downloads/foss4g2013/foss4g2013-pgrouting-for-dummies.pdf As per this document as well as issues raised in the mailing list, this script helps in achieving that. But, I figured out that there are multiple issues in the routing_core_smart.sql file. Some tweaking needs to be done to make it work. And that is what is unknown.


texts: How to add Image to KML file using JAVA?

texts:
I have implemented an Android app which adds user's current location with Lat , Long, Title (Location Title) , Remark & Image for the location (captured with Camera).

Now , What I want to achieve is add Image to KML file & to generate KML file without Image I'm using below code & it is working completely fine .

Here is code to generate KML file which includes lat , long , title & remark :

public void generateKMLFile()
{
try {
    DocumentBuilderFactory factory = DocumentBuilderFactory.newInstance();
    DocumentBuilder builder = factory.newDocumentBuilder();
    TransformerFactory tranFactory = TransformerFactory.newInstance();
    Transformer aTransformer = tranFactory.newTransformer();
    Document doc = builder.newDocument();
    Element root = doc.createElement("kml");
    root.setAttribute("xmlns", "http://earth.google.com/kml/2.1");
    doc.appendChild(root);
    Element dnode = doc.createElement("Document");
    root.appendChild(dnode);
    Element rstyle = doc.createElement("Style");
    rstyle.setAttribute("id", "restaurantStyle");
    Element ristyle = doc.createElement("IconStyle");
    ristyle.setAttribute("id", "restaurantIcon");
    Element ricon = doc.createElement("Icon");
    Element riconhref = doc.createElement("href");
    riconhref.appendChild(doc.createTextNode("http://maps.google.com/mapfiles/kml/pal2/icon63.png"));
    rstyle.appendChild(ristyle);
    ricon.appendChild(riconhref);
    ristyle.appendChild(ricon);
    dnode.appendChild(rstyle);
    Element bstyle = doc.createElement("Style");
    bstyle.setAttribute("id", "barStyle");
    Element bistyle = doc.createElement("IconStyle");
    bistyle.setAttribute("id", "barIcon");
    Element bicon = doc.createElement("Icon");
    Element biconhref = doc.createElement("href");
    biconhref.appendChild(doc.createTextNode("http://maps.google.com/mapfiles/kml/pal2/icon27.png"));
    bstyle.appendChild(bistyle);
    bicon.appendChild(biconhref);
    bistyle.appendChild(bicon);
    dnode.appendChild(bstyle);

    Cursor c = mDatabaseManager.getAllLocations();
    while (c.moveToNext())
    {
        LocationBean lb = new LocationBean();
        lb.locationTitle = c.getString(c.getColumnIndex(DatabaseManager.LOCATION_TITLE));
        lb.lat = c.getDouble(c.getColumnIndex(DatabaseManager.LOCATION_LATITUDE));
        lb.log = c.getDouble(c.getColumnIndex(DatabaseManager.LOCATION_LONGITUDE));
        lb.remark = c.getString(c.getColumnIndex(DatabaseManager.LOCATION_REMARK));

        Element placemark = doc.createElement("Placemark");
        dnode.appendChild(placemark);
        Element name = doc.createElement("name");
        name.appendChild(doc.createTextNode(lb.locationTitle));
        placemark.appendChild(name);
        Element descrip = doc.createElement("description");
        descrip.appendChild(doc.createTextNode(lb.remark));
        placemark.appendChild(descrip);
        Element styleUrl = doc.createElement("styleUrl");
        styleUrl.appendChild(doc.createTextNode( "#" +lb.locationTitle+ "Style"));
        placemark.appendChild(styleUrl);
        Element point = doc.createElement("Point");
        Element coordinates = doc.createElement("coordinates");
        coordinates.appendChild(doc.createTextNode(lb.log+ "," + lb.lat));
        point.appendChild(coordinates);
        placemark.appendChild(point);

    }

    Source src = new DOMSource(doc);
    Result dest = new StreamResult(new File("/sdcard/PlaceMarkers.kml")); // temporarily directly saved to sdcard
    aTransformer.transform(src, dest);
    System.out.println("Completed.....");

}
catch (ParserConfigurationException e) {
    e.printStackTrace();
} catch (TransformerConfigurationException e) {
    e.printStackTrace();
} catch (TransformerException e) {
    e.printStackTrace();
}

}

It generates KML file & I have imported it in Google Maps to check if it works & it is imported perfectly with location title & remark with marker on map.

How do I add an image to this KML file?

Note : location image is stored into SQLite database as BLOB data.


texts: Mysterious Invalid Syntax in Field Calculation

texts:
I'm trying to calculate a field in Model Builder using the Calculate Field tool. I'm using the following code block and expression in the calculation, but am getting an "invalid syntax (, line 1)" error.

# Code Block:
def get_block(block_id):
    block_only = str(block_id)[1:]
    if block_only[-1].isalpha():
        return block_only.zfill(3)
    else:
        return block_only.zfill(2)

# Expression:
get_block(%block_id%)

%block_id% is an in-line substitution that points to a string variable I create as the output of a custom Python script tool. Any ideas what the error is?


texts: Setting Default Label options in ArcMap?

texts: In ArcMap 10.3 I am looking to set the default style for my labels as white with a black halo for any newly added or created layer. Is this possible? I am getting tired of altering each one as they are brought in or created. I have worked with the "labeling" tool but still have to change any incoming layer to the halo style I need.


texts: Getting coordinates of Polyline in ArcGIS Desktop?

texts: I would like to know how to get the coordinates of one of the vertexes of a polyline in ArcMap 10.3.1.


texts: Curve Number equation in ArcGIS Raster Calculator

texts:
I have a floating-point raster of Curve Numbers (70-91) and I'm attempting to use the curve number equation to create a weight raster for flow accumulation, to determine runoff for a single rainfall event. I'm familiar with this curve number equation:

where
Q = total runoff
P = preciptation (1 inch in this case)
0.2S = initial abstractions
0.8 = potential maximum retention after runoff begins

The 2 outside the first bracket means to the power of 2. Firstly, I am unsure how to calculate the power of 2 in the Raster Calculator. Also, would precipitation be an input raster or just a value (1 for example)?


texts: What causes ERROR 000732 from MosaicToNewRaster?

texts:
I'm trying to run the MosaicToNewRaster tool as part of a python script. This is my code:

input_rasters = arcpy.GetParameterAsText(5), arcpy.GetParameterAsText(6), arcpy.GetParameterAsText(7), arcpy.GetParameterAsText(8)
output_location = "C:/Birkbeck_Python/Prog_Assign_1.gdb"
arcpy.MosaicToNewRaster_management(input_rasters, output_location, "R_Dataset", "", "",    "", "1", "", "")

I get the following error message, anybody got any ideas? It works fine when I hardcode it.

Thanks

ExecuteError: Failed to execute. Parameters are not valid.
ERROR 000732: Input Rasters: Dataset '';'';'';'' does not exist or is not supported
Failed to execute (MosaicToNewRaster).


texts: Adding shapefile names in message to progress window/dialog of ArcGIS for Desktop from ArcPy?

texts:
How do you use python in a ArcGIS 10.3 script to list/print the shapefiles with a message from the result of a geoprocessing tool in the progress window/dialog box. For example;

abc.shp
1234.shp
zyx.shp
987.shp

The above 4 datasets have been projected to: WGS_1984_World_Mercator and saved in C:\GIS_Result

I have been looking at this example; http://resources.arcgis.com/en/help/main/10.2/index.html#//03q30000000q000000


texts: Distance matrix (metres instead of degrees) in QGIS?

texts:
I have a problem related to the distance matrix in QGIS (2.12.0). Instead of meters it calculates distances in degrees. I found that some other users faced also this problem. I tried to solve it according to the instructions given by others. Unfortunately it didn’t work in my case. I don’t know where the problem is.

Regarding project properties, I set enable on the fly CRS transformation and chose ETRS89 / Poland CS92 as my project CRS. For Settings->Options->CRS I chose Automatically enable on the fly reprojection by default; and use project CRS for new layers.

I uploaded CSV file with x,y coordinates (the sample: 52.29685 20.93881    DORYCKA 5, Warszawa). And then saved this CSV as ESRI shp (with project CRS). The objects were visualized on the map.

Unfortunately, as I said, I was able only to measure distance matrix in degrees - is it possible to do so in metres?


texts: How to travers fields of records without introducing them all in the function?

texts:
In my previous question I wanted to reach a field in the field calculation; however, it was pretty easy to introduce a few more inputs for that function.

Now, I am facing a problem which needs around a hundred inputs.
It might be a very naive question but I am wondering if there is a way to introduce an input to the function which is an object of the record type and then by a for loop I can traverse through the fields and sum them up.

My data includes the number of people for each age from 0 to 100 year-old and I want to calculate the number of people less than 30, above 30, less than 50 and above 50. So, having 50 inputs for the function and changing them for each new field is time-consuming.  I think there is supposed to be a more efficient way.


texts: Performing Random-Forest Classification of 10cm Imagery for species-distribution in R (no point-shapes)?

texts:
I'm trying to classify 10cm-resolution images with the RF-package.

I do this to separate some species in the forest, seen on my images.

While I was searching for an answer, I found the following code, bit I don't understand why the images have to be transformed into point-shapes, so I'm looking for a way to perform a RF-classification just with Raster (*tif, *png).

The code I found (made by Mr. Evans):

    require(sp)
    require(rgdal)
    require(raster)
    require(randomForest)

# CREATE LIST OF RASTERS
    rlist=list.files(getwd(), pattern="img$", full.names=TRUE)

# CREATE RASTER STACK
    xvars <- stack(rlist)

# READ POINT SHAPEFILE TRAINING DATA
    sdata <- shapefile("inshape.shp")

# ASSIGN RASTER VALUES TO TRAINING DATA
    v <- data.frame(extract(xvars, sdata))

# RUN RF MODEL
## assumeing that sdata has a variable called 'train' with e.g. 0 and 1 values
    rf.mdl <- randomForest(x=v, y=as.factor(sdata$train))
## or for regression:
    rf2 <- randomForest(x=v, y=sdata$train)

# CHECK ERROR CONVERGENCE
    plot(rf.mdl)

# PLOT mean decrease in accuracy VARIABLE IMPORTANCE
    varImpPlot(rf.mdl, type=1)

# PREDICT MODEL
    predict(xvars, rf.mdl, filename="RfClassPred.img", type="response",
        index=1, na.rm=TRUE, progress="window", overwrite=TRUE)

My Questions:

At "Read Point Shapefile Training Data", why can't I use just raster?
I don't get the "jump" to point-shape
In my understanding, I need to create a raster-stack and some training-areas, which are imported as raster's as well. Then I perform rf-classification and plot the results.

Am I right?


texts:
I'm trying to convert a CQL filter to OCG. I followed this example: https://blog.ianturton.com/geotools/code/ogc/2015/03/18/Generating-OGC-Filters.html

In line with org.geotools.xml.Configuration configuration =
               new org.geotools.filter.v1_0.OGCConfiguration(); I get

  "Incompatible types error: Required: org.geotools.xml.Configuration
  Found: org.geotools.filter.v1_0.OGCConfiguration"

What could be the problem, since I'd say that the example provided in the link is working.

My dependency for that part is:

<dependency>
  <groupId>org.geotools.xsd</groupId>
  <artifactId>gt-xsd-filter</artifactId>
  <version>21.0</version>
</dependency>


texts: Viewing multiple rasters from single geopackage in QField

texts:
I'm running into an issue viewing rasters in QField with multiple rasters in a single geopackage.
I export rasters (tif, ecw) to geopackage and they open in QGIS. They're listed as "copy" in QField Sync properties. When I open the QField project, I get an "Unable to load some layers" error saying "the following layers could not be loaded, please review those and reconfigure the QGIS project".
I tried importing them into the QGIS project in different ways (add saved file to map from export dialogue, added from QGIS Browser, dragged from Windows File Explorer), none of which load in QField.
The same rasters and export settings will open in QField if I save each raster to its own geopackage.
Does this have to do with the raster export settings or how the rasters are imported in to the QGIS project?


texts: problems of display modis ndvi with TIMESAT

texts:
I use MODIS MRT tool mosicia tow tiles MODIS NDVI data, then use TIMESAT image view to display the mosaic data, encounting file size waring :

and after drawing data, the sequences of two tiles is wrong, original data display sequence is right to left, but display result is left to right.

but when i use ARCGIS or envi to display the mosaic data, the display sequence is ok.

so i want to know what is wrong with my mosaic data?


texts:
Using the QGIS application, I am able to do this.
However, I am unsure of how to do it using Python.

I have this so far:

from qgis.gui import *
from PyQt4.QtCore import *
from PyQt4.QtGui import *
from qgis.core import *
import time
import sys

def main():
    qgis_prefix="/usr"
    QgsApplication.setPrefixPath(qgis_prefix, True)
    app = QgsApplication([], True)
    app.initQgis()

    canvas = QgsMapCanvas()
    canvas.setCanvasColor(Qt.white)

    fileName = "/..."
    fileInfo = QFileInfo(fileName)
    baseName = fileInfo.baseName()
    path = fileInfo.filePath()
    rlayer = QgsRasterLayer(path, baseName)
    if not rlayer.isValid():
            print "Layer is not valid!"
    QgsMapLayerRegistry.instance().addMapLayer(rlayer)
    canvas.setExtent(rlayer.extent())
    canvas.setLayerSet([QgsMapCanvasLayer(rlayer)])

    canvas.refresh()
    app.exec_()
    app.exitQgis()
main()

Also, I can't get this program to terminate without manually doing it.


texts: Is JAI really that usefull for raster images on Geoserver?

texts:
I've been going back and forth with JAI for quite a while by now. I've never been really impressed with it to be honest. Is it really all that jazz that I hear people talk about?

The problem with JAI is that I'm stuck with a 32bit JVM, thus not fully utilizing the abundance of RAM that I have.

Has anyone done any head-to-head comparisons between 32bit JVM + JAI vs. 64bit JVM?

I have a Geoserver instance with 10k+ geotiff images, and only these geotiff images.
Which one of these setups would give me the best performace in your opinion?

JVM 32 bits
Native JAI installed
-Xmn1000m
-Xms1400m
-Xmx1400m
(among others)

JVM 64 bits
No native JAI
-Xmn2g
-Xms10g
-Xmx10g
-d64
(among others)

Thanks for your responses!


texts: QGIS plugin multithreading and loading memory layers

texts:
I've been working on a small plugin recently and I wanted to add a progress bar to it, since there are tasks in this plugin that take some time to complete. I looked at this and this, but so far I've encountered a problem which I can't seem to solve.

Whenever I run the plugin without moving it to another thread, it runs perfectly fine, everything works (but the GUI is unresponsive for a short while). But when I move the worker (the python class that does the heavy lifting) to another thread, suddenly, I can't seem to add the QgsVectorLayer (which is a memory layer) to the iface.LegendInterface anymore. Everything else still works fine, but the layer just does not show up in the legend in QGIS. So far, this is my code, I've left out the specific methods and imported code, since that all works if I don't move it to a thread (so I figure that can't be the problem):

# import some modules used in the example
from PyQt4 import QtCore
from qgis.core import *
from qgis.gui import *
import utility_functions as uf

import traceback

class Worker(QtCore.QObject):

def __init__(self, iface, data, layer, plugindir, outfile):
    QtCore.QObject.__init__(self)
    self.iface = iface
    self.data = data
    self.layer = layer
    self.plugindir = plugindir
    self.outfile = outfile
    self.killed = False
    self.PathLayer = None

finished = QtCore.pyqtSignal(object)
error = QtCore.pyqtSignal(Exception, basestring)
progress = QtCore.pyqtSignal(float)
layerCreated = QtCore.pyqtSignal(object) # tried for debugging, does not work

def run(self):
    ret = None
    try:
        # Create edgelist
        edgelist = self.selectEdges()

        # Select features on layer and copy to new layer
        features = self.selectFeatures(edgelist, self.layer)

Everything works perfectly fine up until here. In the method 'createNewLayer', the memory layer should be (1) created using QgsVectorLayer, and (2) Loaded using QgsMapLayerRegistry.instance().addMapLayer(layer)

        self.createNewLayer(features)

        # Apply layer style
        if not self.PathLayer:
            lyr = uf.getLegendLayerByName(self.iface, 'Path')
        else:
            lyr = self.PathLayer
        self.applyStyle(lyr)

        if self.killed is False:
            self.progress.emit(len(self.data))
            ret = 1
    except Exception, e:
        # forward the exception upstream
        self.error.emit(e, traceback.format_exc())
    self.finished.emit(ret)

But somehow, the moving to another thread messes with the layer and it can't be added. What's even weirder: I can see the layer when I type into the QGIS python console: QgsMapLayerRegistry.instance().mapLayers(). So I'm guessing that the way QThreads work, messes up the QgsMapLayerRegistry OR the iface.LegendInterface.

I'm working on macOS Sierra, QGIS 2.14, but I've also tried the plugin on a W10 machine, QGIS 2.18, same problems.


texts: Using QGIS Raster Calculator with two different rasters

texts:
I have several maps like this one:

And with the proximity tool and a coastline shapefile I made this output, representing the distance to shoreline:

Now, I want to combine both using the Raster Calculator. However, everytime that I make an operation, the output is just NaN.

Do you have any hint on why is that?

Solved! Both maps had different projections and the "on the fly reprojection" was activated.


texts: Inconsistent export of transparent layers

texts:
I have a map and an inset map. The content of the inset map's data frame is a direct copy and paste of the content of the main map. Both layers have (in order from bottom to top) a base map layer, a polygon feature class with 60% transparency, and a point feature class. When I export the map to any file format (pdf, png, bmp, jpeg), the transparency of the polygon feature class is preserved on the main map but on the inset map the polygon feature class is rendered opaque (0% transparency).

Why aren't both being rendered transparent?

Arcmap:

Export:


texts: How do I change this exif data into long and lat?

texts:
I've been using ImageMagick in a Linux terminal but I'm not sure how to read the coordinates it's giving me.

If you can't see the picture, the data is:

Properties:
    date:create: 2020-05-21T12:33:53+01:00
    date:modify: 2020-05-21T12:33:53+01:00
    exif:ExifOffset: 38
    exif:GPSDateStamp: 2015:12:09
    exif:GPSInfo: 68
    exif:GPSLatitude: 23/1, 21/1, 35928/625
    exif:GPSLatitudeRef: N
    exif:GPSLongitude: 98/1, 0/1, 14886/625
    exif:GPSLongitudeRef: W
    exif:GPSMapDatum: WGS-84
    exif:GPSTimeStamp: 11/1, 21/1, 42/1
    exif:GPSVersionID: 2, 2, 0, 0
    exif:PixelXDimension: 1000
    exif:PixelYDimension: 660
    jpeg:colorspace: 2
    jpeg:sampling-factor: 2x2,1x1,1x1
    signature: 1c3ce39629e01d4cf738f153e810a468264c1dc110a9af7ebe37b83e81969ee3
    xmp:GPSLatitude: 23,21.958080N
    xmp:GPSLongitude: 98,0.396960E
    xmp:GPSMapDatum: WGS-84
    xmp:GPSTimeStamp: 2015-12-09T11:21:42.37Z
    xmp:GPSVersionID: 2.2.0.0
    xmp:PixelXDimension: 1000
    xmp:PixelYDimension: 660


texts: Using Editor tab Clip function in ArcPy?

texts:
I am using ArcGIS Desktop 10.6, basic license, Spatial Analyst, python 2.6, 2019 AD.

I have a polygon feature class with representing a lake, with multiple features overlapping (all in the same feature class). The shoreline feature covers the entire area of the lake, the 1 foot depth feature covers all of the area at 1 foot, etc, etc. Opening an editing session and starting at the deepest depth, I can use the Editor toolbar Clip function to clip each area from the others. The end result is a flat representation of the lake, no overlapping depths.

To accomplish this with arcpy, I think I need to sort my list descending from the largest depth and clip the below features. I've attached a sample feature class that this script would intake, and another ("flat") that should be output.

import arcpy
from arcpy import env

try:
    # Set workspace environment
    env.workspace = "C:/data/default.gdb"

    # intake lake
    in_dataset = arcpy.GetParameterAsText(0)

    # sort the lake attribute table by the depth
    sort_fields = ["depth", "DESCENDING"]

    # Use Peano algorithm
    sort_method = "PEANO"

    # sort the fc, run Editor Clip on each attribute, starting with the largest
    arcpy.Sort_management(in_dataset, sort_fields, sort_method)

    with arcpy.da.SearchCursor(in_dataset, 'SHAPE@') as cursor:
        for row in cursor:
            #don't know how to call this part
            EditorClipTool(row, "DISCARD")

    print arcpy.GetMessages()

except arcpy.ExecuteError:
    # Print error messages
    print arcpy.GetMessages(2)

except Exception as ex:
    print ex.args[0]

input:
https://drive.google.com/open?id=1ShZ5OQensnskLbgIzMWroVM2GIKLukou

desired output:
https://drive.google.com/open?id=1QLbLr30yY136wvAUaEXEHSrM33ZKR3cG


texts: Create a rainfall raster with gstat or idw in R

texts:
I am a beginner in R, I need to use gstat or idw function to get a rainfall raster, but I have an error in the formula. How can I change it?
This is my code:

setwd("C:/Users/A/Desktop/DatiPrecipitazioni/dati")
mydata<-read.csv(file="prova.csv",header=TRUE,sep=";",dec=".",na="")

dataset<-read.csv(file="Dataset.csv",header=TRUE,sep=";",dec=".",na="")

DTM=raster("DTM.tif")

vals<-mydata$HPrec

simplerasterPREC = function(dataset,DTM,vals){

 coordinates(dataset) = ~x+y

 proj4string(dataset)<- proj4string(DTM)

 idw_formula = gstat(id="vals",formula = as.formula(paste(vals,"~1")), data = dataset,nmax=10, maxdist=50000, set = list(idp=2))

 rast_int <- interpolate(DTM, idw_formula)

 rast_int <- mask(rast_int,DTM)

 return(rast_int)
}

If I see what happens line by line I have this error:

coordinates(dataset) = ~x+y

proj4string(dataset)<- proj4string(DTM)

idw_formula = gstat(id="vals",formula = as.formula(paste(vals,"~1")), data = dataset,nmax=10, maxdist=50000, set = list(idp=2))

Error in terms.formula(formula) : invalid term in model formula


texts: QGIS 2.18 Raster Calculator not working properly

texts:
I've been trying to make rasters of the area within a certain elevation band of some elevation rasters I have. Simple enough. But when I do the raster calculation, the output is not what it should be (like in a raster with values 1000-2000 if I do "raster">1500 there should be some some areas >1500 that are one color and other areas that are another color, but the resultant raster is one solid block of a single color)

I attached some pictures of what I did and the result to hopefully show the problem better. As you can see, the resultant vector should be a mix of black and white, because some of the original raster is <1525, but it is not

I'm just trying to figure out how to get these kinds of calculations to work

This was my result

nan is from divide by zero, right? Why would this pop up here? I'm not even dividing anything


texts: How to solve loading performance issues when handling multiple KMLs?

texts:
I am using over 10 KML in Google Maps. When I increase the number of KMLs they load too late.

I use checkboxes to hide/show KML layers. With up to 5 KML layers, it worked fine. When I increase the size of the KMLs it's not working normally.

When I check "show", it shows but "hide" does not work ...

How can I solve this?


texts: GeoServer WMS using NetCDF plugin and GRIB2 time dimension

texts:
I am attempting to use GeoServer with the NetCDF plugin to create a WMS-T (WMS with time dimension) service. My input is a GRIB2 file containing a single parameter (temperature) and multiple forecast hours.

I have no issues when using a NetCDF file, but I can't get the time dimension to work when using a GRIB2 file. With the GRIB2 file, the only time I can get to display is the initialization time, not the forecast times.

Is this a known issue with the GRIB2 implementation of the NetCDF plugin? If not, does anyone have an idea of what I could be doing wrong? I have also tried the GRIB2 plugin, but that did not work either.

Here is the GRIB2 file that I am using:
https://gofile.io/?c=9yHAEA


texts: How do I cluster by features in CartoDB?

texts:
I have a set of points for a number of buildings in NYC. I would like to group those points into clusters per their respective neighborhood (ex: Chelsea, Upper East Side etc.). Each of my points have a designated neighborhood.

Is there a way to tweak the CSS to make it happen? I don't quite understand how the Wizard is currently clustering my points:

    /** cluster visualization */

#data{
  marker-width: 12;
  marker-fill: #FD8D3C;
  marker-line-width: 1.5;
  marker-fill-opacity: 1;
  marker-line-opacity: 1;
  marker-line-color: #fff;
  marker-allow-overlap: true;

  [src = 'bucketE'] {
    marker-line-width: 5;
    marker-width: 24;
  }

  [src = 'bucketD'] {
    marker-line-width: 5;
    marker-width: 34;
  }

  [src = 'bucketC'] {
    marker-line-width: 5;
    marker-width: 44;
  }

  [src = 'bucketB'] {
    marker-line-width: 5;
    marker-width: 54;
  }

  [src = 'bucketA'] {
    marker-line-width: 5;
    marker-width: 64;
  }
}

#data::labels {
  text-size: 0;
  text-fill: #fff;
  text-opacity: 0.8;
  text-name: [points_count];
  text-face-name: 'DejaVu Sans Book';
  text-halo-fill: #FFF;
  text-halo-radius: 0;

  [src = 'bucketE'] {
    text-size: 12;
    text-halo-radius: 0.5;
  }

  [src = 'bucketD'] {
    text-size: 17;
    text-halo-radius: 0.5;
  }

  [src = 'bucketC'] {
    text-size: 22;
    text-halo-radius: 0.5;
  }

  [src = 'bucketB'] {
    text-size: 27;
    text-halo-radius: 0.5;
  }

  [src = 'bucketA'] {
    text-size: 32;
    text-halo-radius: 0.5;
  }

  text-allow-overlap: true;

  [zoom>4]{ text-size: 16; }
  [points_count = 1]{ text-size: 0; }
}


texts: Adding Points to line at certain distance using QGIS?

texts:
USING QGIS 2.18
This question seems to have been asked a ton with lots of options.i just need to add some points/nodes/vertex to a line at a recurring distance along that line.

There are some options for use with PostGIS (Getting points along line at defined distance using PostGIS?) as well as an idea to add buffers (Creating points along line at varying distances using QGIS?) at each distance needed, but this would not be useful on a long line.


texts: Access ArcGIS map created on localhost using JavaScript

texts:
I've created an ArcGIS map on my local ArcGIS map server using ArcMap desktop and I am able to open up the map from "View In: ArcGIS Online map viewer". When viewing the map, the url is something like:

www.arcgis.com/home/webmap/viewer.html?url=http%3A%2F%2Flocalhost%3A(...etc) --> simply the local url pointing to the map on my local rest service.

Can I modify this map using JS from here? How should I access this map and what is the map id in this case or do we even need it?

Also when I do "View In: ArcGIS JavaScript", it only shows the frame and the map doesn't come up


texts: Is there a limit on number of fields imported in QGIS? If yes, how to disable it?

texts: I'm working on a dataset with around 200 fields. However, when I import the data (in the form of .csv file) as a layer, the attribute table does not display all the fields – around 100. Is there a limit on the number of fields that can be imported in QGIS? More essentially, how do I remove this limit (or set it to a comfortably high number) using python?


texts: Precise st_area() of small polygon need geography


texts:
I am using leaflet 1.0.0 to display the GeoJSON. The GeoJSON contains data for the whole bounding box but I want to display the feature in some location only.

Is it possible to clip the GeoJSON with another GeoJSON or something like that to display on the leaflet map?


texts:
I try to open modis aqua data (HDF5 format: A2015208112000.L2_LAC_OC.nc) on QGIS 2.8 ver. And an error appears: No extra dimensions found.

Could someone help me to solve this problem?


texts: v.generalize creates SHP instead of GPKG

texts:
I'm trying to use generalize.smooth, but the time I try to save in GPKG (GeoPackages) it always saves in SHP.
The reason for saving in GPKG is that the file has more than 5GB and is impossible to save in shape.
Generalize has the save option in GPKG, but the output is always in SHP.


texts: I'm using the intersection tool (processing toolbox) to intersect polygons. I have a huge number of small polygons to intersect and it needs so much time to do do that. Maybe it is caused by the great extent of the polygon shape. Is it possible to set a maximum distance which defines the max distance between polygons they will be intersected. I would like to put this option into a python script if it is possible.


texts: How to extract a shapefile with a SQL expression from a database of shapefiles using PyQGIS or OGR

texts:
I have a shapefile Water_Districts.shp that consists of several water district polygons. I am trying to use the script below to 'extract' a specific polygon with the name Cache La Poudre River listed under the NAME column in the Water_Districts.shp attribute table. I have been using the python script below, but am not sure how to incorporate the SQL expression in the -where exception of the ogr2ogr command. I am getting the error: FAILURE: SetAttributeFilter(NAME = Cache La Poudre River) on layer'Water_Districts' failed.

import gdal
import subprocess

directory = r"F:\IrrigatedLands\FC_test\Water_Districts\Water_Districts.shp"
output_shp = r"F:\IrrigatedLands\FC_qgis\boundary.shp"

subprocess.call(["ogr2ogr", "-where", "NAME = Cache La Poudre River", output_shp, directory])


texts: Routing with osm2po - starting from middle of the street

texts:
I'm currently looking at using osm2po's java library to build a simple application that builds a shortest path between my favorite bars :)

However, I'd like to start at any location, not at an intersection.

I saw a response that's it's theoretically possible (http://www.rqna.net/qna/nwkrrs-can-i-start-and-end-an-osm2po-route-in-the-middle-of-a-road-segment.html), but not really sure how to do this with the latest java API.

Is this functionality exposed? I've seen findClosestVerexId(lat, lon), but what I really want is something like findEdge(lat, lon), and then a way to get source and target from the edge.

Thanks!


texts:
arcpy.TableSelect_analysis(in_table, out_table, {where_clause})

"Selects table records matching a Structured Query Language (SQL) expression and writes them to an output table."

Using the arcpy.TableSelect_analysis() function I have this simple script:

import arcpy
from arcpy import env

env.workspace = r'\\gisfile\GISstaff\Jared\ModelBuilder\JaredTest.gdb'

#variables
intable = r"\\gisfile\GISstaff\Jared\ModelBuilder\JaredTest.gdb\ElectionResults_Nov2016"
outtable = r"\\gisfile\GISstaff\Jared\ModelBuilder\JaredTest.gdb\Results3"
where_clause = "ContestTitle = 'BALLOTS CAST - TOTAL'"

arcpy.TableSelect_analysis(intable, outtable, where_clause)

It selects all records named 'BALLOTS CAST - TOTAL' from the 'ContestTitle' column of my 'ElectionResults_Nov2016' table and writes them in the newly created 'Results3' table. 'BALLOTS CAST - TOTAL' is only one of 58 separate table record names.

How can I write all 58 to their own tables?


texts: MMQGIS Google Geocoder limiting me to less than 10 addresses

texts:
I've been using MMQGIS geocoding with Google Maps as the geocoder for a little while now, never any issue. As of yesterday, have gotten a few strange issues. I was unable to geocode at all yesterday, as it was giving me an HTTP 503 error when trying to connect to the google api. I figured this was just a result of maintenance or whatnot, so I used the openstreetmap geocoder, which is unfortunately far less accurate.

I tried using the google geocoder today, and it's even stranger, it successfully geocodes some of my addresses, but usually quits before I get to 10 of them, saying that I've reached my limit for the day. Not only is this impossible, it's also said this on different amounts of addresses geolocated (I've had 0, 2, 6, 7, and 10 so far). Has anyone else been having this issue, or does anyone have any possible solutions?


texts:
We have a map with streets as basemap.

We have been trying unsuccessfully, for the last two days to change the basemap to tiled layer with no success at all.

I am posting a snippet because script is long. Please let me know if you need additional script.

On the surface, it looks simple but it has been extremely challenging.

Your generous help is greatly appreciated in advance.


texts:
I have a shapefile from the Census Bureau on Zip Code Tabulation Areas (http://www.census.gov/cgi-bin/geo/shapefiles/index.php). I wanted to get, for each zip code, the list of neighboring zip codes.

More specifically, I would like to export a file that gives a dataset of the form,

ZIPCODE,NEIGHBOR1,NEIGHBOR2,NEIGHBOR3,NEIGHBOR4
90001,90002,90010,.,.
60640,60660,60625,60613,60618

So for each zip code, there would be columns which would be populated with neighboring zip codes. The column would be set to missing for zip codes that do not have as many neighbors.

Basically, I want to do this (Finding all neighbors for all polygons in Spatialite) in qgis. I guess I'm asking whether the command in the link would generate what I want described above, and if it does how I can run it. I have PostGIS installed at this point, but I do not know how I can load the vector on it using QGIS and PSQL...

Basically, I am very unfamiliar with QGIS, except being able to load a vector layer using .shp file. I would like to know how to get to the part where I can use the command like

SELECT p1.ROWID, p1.ZCTA5CE10, p2.ZCTA5CE10
FROM "US_zcta5_2010" p1, "US_zcta5_2010" p2
WHERE Touches (p1.Geometry, p2.Geometry)=1;


texts:
I have a point layer and three different network layers and am trying to create OD-Matrices for all the points using the QNEAT3 plugin.

It did work out fine on two of the three networks, but the third always returns just NULL for all the relations.
The networks that work are buslinenetworks, one taken from a local bus companies opendata service, the other one I created myself. The third represents the streets of the city. The data is taken from OpenStreetMap.

When I simply try routing between two points, it works out fine and returns the entry, network and exit costs for the relation. Only when using the algorithms to create a martrix I get empty results for all the relations.

Did anyone ever have a similar problem or might know where to find the problem?


texts: How to make a Standalone QGIS application run as an exe file

texts: I had made a simple standalone application and that works fine. But the user have to execute a batch file(inside in which I had set the QGISHOME path, Python path and call my main python file) to run the application. I am assuming there are ways to execute the application as an EXE file. If so what are the steps to convert to exe?


texts: Add new field in PostgreSQL row while using ogr2ogr for shapefile to PostgreSQL

texts:
I am using the following ogr2ogr command for putting my shapefiles into PostgreSQL table. This command runs fine.

ogr2ogr.exe -append -lco GEOMETRY_NAME=geom -lco SCHEMA=public -f "PostgreSQL" PG:"host=localhost port=5433 user=postgres dbname=abc123 password=xxxxxx" myshapefile.shp -nln shapetable

Now I want is to add extra field/fields like (countryName=xyz / RoadName= street123/etc) while using the above command when exporting my shapefile to PostgreSQL table.

I was trying to find the answer in another link but it was for adding extra field while exporting GeoJSON into PostgreSQL table.


texts: Google Earth Engine: reprojecting over an image collection

texts:
Mapping over an image collection (NO2 in this case), for .clip, is simple:
var NO2NL=NO2.map(function(im){
   return im.clip(NL);
});

Replacing .clip with .resample(bilinear) works too, however using .reproject(scale=1000) always produces an error: Variable must be declared with "var": scale. Any idea how I can reproject over an image collection?
Entire script for reference:
//Set variables

var countries = ee.FeatureCollection("FAO/GAUL/2015/level0");
var NL = countries.filterMetadata('ADM0_NAME', 'equals', 'Netherlands');
var start_period = ee.Date('2018-07-01');
var end_period = ee.Date(new Date().getTime());

var NO2 = ee.ImageCollection('COPERNICUS/S5P/OFFL/L3_NO2')
    .select('tropospheric_NO2_column_number_density')
    .filterDate(start_period, end_period);

var NO2NL=NO2.map(function(im){
   return im.clip(NL);
});

//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

// Date slider parameters
ee.Dictionary({start: start_period, end: end_period})
  .evaluate(renderSlider);

function renderSlider(dates) {
  var slider = ui.DateSlider({
    start: dates.start.value,
    end: dates.end.value,
    period: 14, // Every 14 days
    onChange: renderDateRange,
    style: {
    position: 'top-left',
  }
  });
  Map.add(slider);
}

function renderDateRange(dateRange) {
  var image = NO2NL
    .filterDate(dateRange.start(), dateRange.end())
    .median();

//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

 // Visual parameters and mapping
  var vis = {min: 0, max: 0.0001, palette: [
    'black', 'blue', 'purple', 'cyan', 'green', 'yellow', 'red']};
  var layer = ui.Map.Layer(image, vis, 'Tropospheric NO2');
  Map.layers().reset([layer]);
}
Map.setCenter(5.5, 52.2, 7.5);

Eventually I hope that the time slider shows only dates on which the data is available. I also want to use ECMWF data for wind speed and join the datasets so that only days with wind speed <10kph are shown!


texts: Unexpected results with Raster Terrain Analysis Plugin

texts:
In QGIS 2.18.18 on Windows 10, I have downloaded the NASA file srtm_05_01.tif that is available at http://dwtkns.com/srtm. I use the Terrain Analysis Slope method on that file with the following dialog box values

I use a z factor of 111,120 because the CRS is WGS84 EPSG:4326.  The resulting raster file has the following histogram

where all degrees have 250,000 as a count.  When I use the resulting raster in the raster calculator and select for an evenness of < 5 (i.e. "slope@1" < 5) I get the same output as selecting for an evenness < 87.

If this is not the expected behaviour, or if the srtm_05_01 file is unusual in some way, can someone suggest a file that would work well with the slope method?


texts: Extracting RGB values from aerial photo and assigning to LAS file with FME?

texts:
I need to extract the RGB Values from imagery and assign to a .las file in order to display 3d landscape with "real" colours. I have an unlicenced LAStools, FME 2017 and ArcGIS 10.3.1 at hand.

I have tried with las2las, but the tools don't take the option --color-source as was documented in liblas.org. With FME I can visualize this effect but not export to .las.

Has anyone done this?


texts: How to send multiple value for one param to mapserver through openlayers (encode_uri)

texts:
Here's my filter on Mapserver's Side :

PROCESSING "NATIVE_FILTER=(ct in ('%ct%'))"

It allows me to select several values for my column. My validation looks like :

VALIDATION
    'ct' '^[a-zA-Z\-]+$'
    'default_ct' 'a'
END

On Openlayers side, I get stuck to send several value though. If I do this :

var wmsSource = new ol.source.TileWMS({
   url: 'http://localhost:5000/cgi-bin/mapserv.exe?map=/ms4w/apps/tutorial/htdocs/essai.map&SERVICE=WMS&VERSION=1.1.1%20&REQUEST=GetCapabilities',
        serverType: 'mapserver',
        params: {'LAYERS': 'layer','ct':"value1,value2",'TILED': true}
  });

My comma is encoded while asking mapserver as you can see on that url :

http://localhost:5000/cgi-bin/mapserv.exe?map=/ms4w/apps/tutorial/htdocs/essai.map&SERVICE=WMS&VERSION=1.1.1%20&REQUEST=GetCapabilities&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&FORMAT=image%2Fpng&TRANSPARENT=true&LAYERS=later&ct=value1%2Cvalue2&TILED=true&WIDTH=256&HEIGHT=256&CRS=EPSG%3A3857&STYLES=&BBOX=-5009377.085697312%2C5009377.085697312%2C-9.313225746154785e-10%2C10018754.171394624

My comma has been transcoded in %2C. Is there a way to prevent that, or to solve it?

EDIT : That was actually a validation error. I made the following changes :

var wmsSource = new ol.source.TileWMS({
   url: 'http://localhost:5000/cgi-bin/mapserv.exe?map=/ms4w/apps/tutorial/htdocs/essai.map&SERVICE=WMS&VERSION=1.1.1%20&REQUEST=GetCapabilities',
        serverType: 'mapserver',
        params: {'LAYERS': 'layer','ct':"('value1','value2')",'TILED': true}
  });

and in my mapfile :

VALIDATION
    'ct' '^[a-zA-Z\-\,\\\(\)\']+$'
    'default_ct' 'a'
END

As I need to always return something (see that question) my FILTER look like that :

PROCESSING "NATIVE_FILTER=((ct in %ct%) or (%ct% in 'a'))"

I've got that error :

  msPostGISLayerWhichShapes(): Error (ERREUR:  erreur de syntaxe sur ou près de « 'a' »
  LINE 1: ...ainer','fishing')) or (('container','fishing') in 'a')) and ...

(== Error syntax.)


texts:
I'm trying to add a new attribute to a vector layer in my plugin. But it does not work.. Why?

selectedLayerIndex = self.dlg.comboBox.currentIndex()
selectedLayer = layers[selectedLayerIndex]
name = "newAttribute"
selectedLayer.dataProvider().addAttributes([QgsField(name,QVariant.Double)])
self.selectedLayer.updateFields()
#or
self.selectedLayer.commitChanges()


texts: QGIS crashes when using dxf to shp converter

texts: this is my first question on Stackexchange.  I use a GPS to capture data in .dxf, manipulate it in a CAD program that I know how to use, then the dxf2shp converter (I think Python is mentioned) to let me generate .shp files for the end user. Every time I've done this recently I've got a meesage from QGIS that it has stopped working and will close.  I've tried re-installing QGIS (1.8.0), but still have the same message.  Am I doing something wrong or is this a known problem.  Please note I'm not really computer literate, so if you could avoid assuming I know what a technical term means, i'd be grateful.


texts: Is it possible to import OSM data based on coordinates into QGIS3?

texts:
I'm trying to import Openstreetmap building polygons to QGIS 3.4 based on coordinates: I have coordinates of sampling spots and would like to import all buildings within the vicinity of say 1 km of these spots.

I have tried to just import all building polygons for an entire region (in Germany) with the idea to make a buffer around sampling spots and save only buildings that are subsequently in the vicinity but my (half-year old) computer is to slow to compute that.

I have also tried to use the QuickOSM plugin but it apparently only accepts the name of cities or villages to load data and not coordinates.

How do I (batch) download building polygons from OSM around certain coordinates?


texts: Why doesn't my vector layers show up when I put them under my raster layer in QGIS 3.12?

texts: My vector layers do not appear in the raster year, do you know how could I fix this?


texts: Clipper tool problem in QGIS when it was used to clip JPEG

texts:
If I draw a vector，that vector can also be used to clip an image (JPEG)，right?

In QGIS I got this error:

Traceback (most recent call last):
File "C:/PROGRA~1/QGIS2~1.14/apps/qgis-ltr/./python/plugins\GdalTools\tools\widgetPluginBase.py",
line 88, in onFinished self.base.onFinished(exitCode, status)
File "C:/PROGRA~1/QGIS2~1.14/apps/qgis-ltr/./python/plugins\GdalTools\tools\dialogBase.py",
line 205, in onFinished msg = unicode(self.process.readAllStandardError())
UnicodeDecodeError: 'ascii' codec can't decode byte 0xd6 in position 42: ordinal not in range(128)

Even changing the saving style into jpeg, it still doesn't work.


texts: Build lookup table from multiple input attributes

texts:
In FME, how could I build a lookup table with multiple input attributes, such as the one presented here: https://community.safe.com/s/article/advanced-schema-mapping-schemamapper which looks like:

and which comes from an input table defined by these fields:
MaterialTypeAttribute: The name of an attribute which defines the pipe material.
MaterialTypeValue: Potential values of pipe material.
DiameterAttribute: The name of an attribute which defines the pipe diameter.
DiameterValue: Potential values of pipe diameter.
StatusAttribute: The name of an attribute which defines the pipe status.
StatusValue: Potential values of pipe status.
?


texts: Automatically georectify raw Landsat scenes

texts:
I am working with 20 raw Landsat images downloaded from GloVis, in WGS 84 (EPSG 32634) coordinate system. I need to geometrically correct them to S_JTSK (EPSG 5514).
Usually, georectification is done manually using set of control points (same location in georeferenced map and image I want to georectify). However Landsat resolution is so poor (30x30m) that I can´t find easily these points.

In this reason I would like to know if there is some automated approach to georectify Landsat scenes?

maybe using scene corners as a Control points?

but how to find them in my S_JTSK map?

I can´t use ArcGIS Tool ´Georeferencing a raster automatically´ because my Landsat image has already a coordinate system.

I am working with ArcGIS, ERDAS and R.


texts: No data found error when reconciling version in sde

texts:
I got an error when reconciling a version in the sde database and it seems sde cannot find the data in one table with a specific STATE_ID (for example STATE_ID = 12345 for table abc).

In the delta table of table abc, all records have a same STATE_ID different from 12345(say 67890). Before the issue occurs, I:

reconciled and posted all versions and compressed database
unregistered table abc and made some schema changes
registered table as versioned.

We use Oracle 11g behind sde. Anyone have an idea how I can resolve this and reconcile the version?


texts:
I have coordinates for an area in Japan, where the survey has used the coordinates -8073.896, 63828.792. I'm not sure what coordinate system this is however I'm trying to convert this using online software to correct for JGD2011 however this plots points with a ~30km offset from where they should be.

Does anyone know how to fix this, or convert the coordinates so that they can be used in Google Earth?


texts:
I am new to working with spatial data in R and I have run into a few issues, that I was hoping the community could help me with.
Here is the situation:
I generated a raster file of a country with a range of probabilities. Within that country, I have state-level data -  # of farms.
My goal is to randomly seed a certain number of points (farms) (# of farms  is stored in CSV file, “Table2”) in each state of the country, based on the probability, indicated in the raster file. As such, I clipped my raster probability file up into 16 different raster files (the number of states in this country).

I have my code written, for what I describe above. But this code needs to be performed over each of the 16 layers - the ‘ncells’ and “TotalFarm’ value needs to be updated for each file.

install.packages("raster")
install.packages("sp")
install.packages("rgdal")

setwd("~XXXXXX”)
r <- raster("Suitability_Prob_v30.tif") #Reading in 16 rasters

r_31 <- raster("Suitability_Prob_v31.tif")
r_32 <- raster("Suitability_Prob_v32.tif")
r_33 <- raster("Suitability_Prob_v33.tif")
r_34 <- raster("Suitability_Prob_v34.tif")
r_35 <- raster("Suitability_Prob_v35.tif")
r_36 <- raster("Suitability_Prob_v36.tif")
r_37 <- raster("Suitability_Prob_v37.tif")
r_38 <- raster("Suitability_Prob_v38.tif")
r_39 <- raster("Suitability_Prob_v39.tif")
r_310 <- raster("Suitability_Prob_v310.tif")
r_311 <- raster("Suitability_Prob_v311.tif")
r_312 <- raster("Suitability_Prob_v31.tif")
r_313 <- raster("Suitability_Prob_v313.tif")
r_314 <- raster("Suitability_Prob_v314.tif")
r_315 <- raster("Suitability_Prob_v315.tif")

TABLE2 <- read.csv(file=“XXXXX”, header=TRUE, sep=",") #Read in cvs that has  farm #'s

Here performing ONLY on the first of the 16 raster files
Set NAs to 0

r[==-3.402823e+38] <- NA
r[is.na(r)]=0

Calculating number of cells

ncells<-dim(r)[1]*dim(r)[2]
ncells

Calling in the # of farms for this state, stored in CSV

TotalFarm<-sum(TABLE2[,'X30'])
TotalFarm

For point  selection, set number of cells (ncells) and the number of points to be selected, (specific to each state)
Weighted by the value in the cells:

ptscell = sample(ncells, TotalFarm, prob=r[], replace=TRUE)

Distribute the points throughout the grid cells
Get the cell half-width:

hs = res(r)/2

Now find the center of those points:

centers = xyFromCell(r,ptscell)

Generate random uniform points in the cell by using the center and the
half-width/height

pts = cbind(runif(nrow(centers),centers[,1]-hs[1],centers[,1]+hs[1]),
            runif(nrow(centers),centers[,2]-hs[2],centers[,2]+hs[2]))

Print and finish:

plot(r)
points(pts)
print(pts)
write.csv(pts, file = "Points_5.txt",row.names=FALSE)

My attempts to apply this code over the series of 16 states has failed, in part, because I cannot get my 16 raster files into a “raster  object” that allows me to use the rasters individually with those changes. Specifically, I made a mosaic image out of the 16 rasters , however the changes that I made on the mosaic were not made on the original raster files (only on the new mosaic object), so when I need to go back and seed the files individually with the number of farms specific to each state, the modifications that I do to the mosaic, i.e. NA’s, are not there.

All of the tutorials that I have seen perform loops on rasters by building a brick or stack. But since these are not bands which have similar extents, rather, they are neighboring states, this approach is not appropriate here.

So, how do I perform a loop on rasters, that are in the “raster” class or otherwise apply the code written above, to each raster automatically?


texts: How do I find the overlapping area between polygons on 2+ vector layers in QGIS when the polygons on all layers are categorized?

texts:
Context

I am essentially trying to find the overlapping area between various categories of circles and various categories of polygons.

I looked for a good while at lots of posts on stack exchange and other websites but can't seem to find one that helps with the caveat I have below and the fact that I have categories.

And just as an FYI, I am extremely new to GIS in general and have been going at this problem for a couple days now. I am using QGIS 2.14.15.

Details

First, I have a vector layer with all the counties in Texas. Each county polygon has been categorized into 4 Tiers.

Second, I have 3 vector layers (could be combined later) that are essentially point coordinates that I buffered. Each resulting circle is categorized into 3 Coverage Types.

Just in case I did not explain that well, here is a picture:

I would like to end up with the data to answer the following questions:

How much Tier 1 area is covered by Coverage Type 1?
How much Tier 2 area is covered by Coverage Type 1?
How much Tier 3 area is covered by Coverage Type 1?
How much Tier 1 area is covered by Coverage Type 2?
Etc... for all Tiers and Coverage Types

And also:

What percentage of Coverage Type 1 Area (total) is Tier 1?
What percentage of Coverage Type 1 Area (total) is Tier 2?
Etc... for all Tiers and Coverage Types

I imagine the output will be in some sort of matrix/table format?

CAVEAT - If circles of the same type overlap, I want to make sure I do not double count (maybe there is a way to combine them temporarily?)


texts: Snapping end points of line to point using PostGIS

texts:
I have a problem with over 10k lines from linetable that do not end/start in a point from pointtable.
I need to get these lines snapped so I tried to intersect them by using st_snap
update linetable c set geometry = st_snap(c.geometry, (select geometry from pointtable d where st_dwithin(c.geometry,d.geometry,0.1)  limit 1) , 0.1)
where  not (st_intersects(ST_StartPoint(c.geometry), (select st_union(geometry) from pointtable ) )
or st_intersects(ST_endPoint(c.geometry), (select st_union(geometry) from pointtable ) )   )

The problem is that it just moves the line to intersect the pointtable, instead of just moving the start/endpoint.
How do I make only the start/end point intersect, rather than adding a vertex or moving a random node?


texts:
I'm trying to use the OSTN15 transformation to convert ETRS (EPSG:4258) to OSGB (EPSG:27700).
I have downloaded the NTv2 data from the OS website and extracted into the proj data directory. I have tried a couple of different proj4text values and none give the expected results.
Using the input point 51.40078220140, -3.55128349240 from the User Guide, the correct result should be 292184.87 168003.47 ±0.01.
SELECT st_astext(st_transform('SRID=4258;POINT(-3.55128349240 51.40078220140)', 27700));

The builtin 27700 uses this one:
+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +datum=OSGB36 +units=m +no_defs
Result: 292183.85 168003.44
Just adding +nadgrids to that and it's completely wrong:
+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +datum=OSGB36 +units=m +no_defs +nadgrids=@OSTN15_NTv2_ETRStoOSGB.gsb
Result: 292100.78 168057.91
Based on this blog post I tried this and it gives the same result:
+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +units=m +no_defs +nadgrids=@OSTN15_NTv2_ETRStoOSGB.gsb
Result: 292100.78 168057.91
I'm afraid I have no idea what I'm doing wrong. Are additional parameters or changes needed to the proj instructions, or have I missed something elsewhere?


texts:
I am trying to use v.vol.rst in Grass 7.0.4 to interpolate temperature values considering elevation.
The answer of Grass to my attempts is: WARNING: Unable to create cellout raster map without cellinp
It is first time with Grass, so I have no idea what might be wrong.

Here is the code from the command.

(Mon Jun 06 21:06:18 2016)
v.vol.rst --overwrite --verbose input=weceeurope_current_gsc@wce wcolumn=temperatur maskmap=wecereurope@wce elevation=wece_europe_current_3D
WARNING: Vector is not 3D
16 records selected from table
WARNING: Strip exists with insufficient data
Processing all selected output files will require 10500 bytes of disk space for temp files
WARNING: Points are more dense than specified 'DMIN'--ignored 1 points (remain 15)
WARNING: 15 points given for interpolation (after thinning) is less than given NPMIN=200
WARNING: There is less than 700 points for interpolation, no segmentation is necessary, to run the program faster, set segmax=700 (see manual)
Bitmap mask created
WARNING: Unable to create cellout raster map without cellinp
Finished interpolating
3D raster map <wece_europe_current_3D> created
The number of points in vector map is 16
The number of points outside of 2D/3D region 0
The number of points used (after reduction) is 15
(Mon Jun 06 21:06:18 2016) Command finished (0 sec)

I tried to use all output options, but result was the same or I got warning about more setting for cross-validation


texts: Deleting lines after manual union in ArcMap?

texts:
I have a question regarding the editing tools. I have made this selections of lines then used the "union" tool to make them into one entity but the lines that were merged together are still on my shapefile "under" my bigger entity. (see picture)

Is there a way to delete all of the lines I used to create the bigger entity without having to select them one by one ?

I can't select in my attribute table because I have other lines that I need to keep.


texts: Search Points within Buffered Polygon with PostGIS

texts:
I have a table of boundaries as polygons and another table of peaks as points. I want to find the points that are within the polygons and the polygons are extended with a buffer. What would be the fastest way to do that?

At the moment I am trying to work with ST_Within which gives me good results, but adding ST_Buffer does not seem to work:

SELECT p.name, p.elev, a.name FROM eu_peaks p, eu_admin_valid a where
    ST_Within(p.geom, ST_Buffer(a.the_geom, 0.02))
    order by a.name;

Is this a good approach? What should I change?


texts: Why does order of commands in raster processing matter?

texts:
My company has recently changed to ArcGIS 10.5 (from 10.1) and now some arcpy scripts that worked fine before don't run anymore. I am not completetly convinced this is actually the reasons of the problem, or whether the data/projection is somehow special but here it goes:

MWE:

import arcpy

arcpy.CheckOutExtension('Spatial')

def create_polyoutline(raster):
    print "creating outline for " + str(raster)
    ras = arcpy.CreateUniqueName("ras_int")
    arcpy.gp.Reclassify_sa(raster,"Value","0 9999999 1",ras,"DATA")
    arcpy.Delete_management(ras)

arcpy.env.workspace  = r"Z:\\"

in_ = r"Z:\ref\d_max_a.tif"
poly = create_polyoutline(in_)
in_ = r"Z:\ref\d_max_b.tif"
poly = create_polyoutline(in_)

this works if i first process grid "d_max_a.tif" and then "d_max_b.tif". But not the other way around.

The error is the following:

Traceback (most recent call last):
  File "<module2>", line 17, in <module>
  File "<module2>", line 8, in create_polyoutline
  File "c:\program files (x86)\arcgis\desktop10.5\ArcPy\arcpy\geoprocessing\_base.py", line 510, in <lambda>
    return lambda *args: val(*gp_fixargs(args, True))
ExecuteError: ERROR 999999: Error executing function.
This spatial reference object cannot be defined from the available information.
Failed to execute (Reclassify).

Both rasters have the same coordinate system assigned (CH1903+_LV95).

In other, more complex scripts the tool sometimes just fails with generic error 999999:

ExecuteError: ERROR 999999: Error executing function.
Failed to execute (Reclassify).

The 2 grids are basically the same, except for slightly different extents.
The problem seems to always arise after processing several grids (not necesarilly 2). However the grid that causes the script to fail can always be processed individually, after restarting the script.
Is there some kind of cache or in_memory magic at work, or why does the order of these 2 reclassify commands matter?


texts: Limiting search radius and improving performance of ArcSDE for Oracle?

texts:
I have a feature class with around 1 million polygon features. I'm using the following query to return features from the polygon feature class that intersect a specified point location (x, y not a point feature).

select * from polygon_feature_class where st_within(sde.st_geometry (x-value, y-value, null, null, srid), shape-field) = 1;

The problem is this query is very slow. It takes 10+ minutes to run through the table. Is there a way I can limit the search radius to say 1000 meters?

In my mind, the process would be "give me the features this point intersects but stop searching after you exceed a 1000 meter radius from the initial point".

In lieu of a way to limit the search radius, how can I improve this query?

I'm stuck the feature class and all its features and I need to start with a point specified as x,y. How I get the polygon(s) that contain the point is open but it needs to be much faster than what I'm currently using.

The query will run against ArcSDE running on Oracle, features stored as st_geometry. Eventually, this SQL will run from an ArcPy script but I'm currently using TOAD.


texts:
I am having trouble re-classing a field. I have two field one is called Rowid and one is Value.

Rowid is being reclass from 0 to 1

Value is being reclass from 1 to 0

I was able to reclass value, however, i am having trouble re-classing Rowid. Below is the script i used for Value. i am using the same for Rowid of course changing the field name and numbers.

import arcpy

#set env. settings

rasterPath = "C:/Users/Yulissa Arancibia/Desktop/Python Coding/Lesson_7/Lesson_7_Data/reclassriv"
inDEM = arcpy.Raster(rasterPath)

arcpy.CheckOutExtension("Spatial")

#set local variables

remapList = []

remapList.append ([1,0]) <---- here for Rowid i put ([0,1])

remapRangeValues = arcpy.sa.RemapRange(remapList)

elevReclass = arcpy.sa.Reclassify(inDEM, "VALUE", remapRangeValues, "NODATA") <--- changed the field to "Rowid"

elevReclass.save("C:/Users/Yulissa Arancibia/Desktop/Python Coding/Lesson_7/Lesson_7_Data/reclassriv2")
# Check out the ArcGIS Spatial Analyst extension license

arcpy.CheckInExtension("Spatial")


texts: Extracting substring from field in ArcMap using Python Parser of Field Calculator?

texts:
I have some parcel data where I need to extract a subdivision name from a long string. The format is always "Subdivision: ____ ______ _____" etc. BUT, there is no uniformity to what comes before or after "Subdivision: " or the actual name of the subdivision. In my example below it shows that "Block: " follows "Subdivision: " but that's not always the case.

I'd like to learn how to solve this issue using python, but VB can also be used. I was reading about re (edit: Regex) in python, but without some further explanation I'm a little lost. Here is a screen shot showing what the data looks like.

Any tips on where I should try to go with this?


texts: how do I reestablish all feature class references in mxd after their names change?

texts:
I have 200 shape files which have been named using an integer at the beginning of the name...i.e 01.03.shp 01.04.shp etc. They are all symbolized within an mxd document.

I need to import these shapes into a file-GDB which is not permitted due to the fact that the feature class starts with a number - I can iterate through the shapes and rename them using a prefix, but then i´ll need to change the mxd reference for every single feature class because the names are different. Is there a clever way to get around having to fix every single reference in the mxd?


texts: Creating copies of feature point geometries based on value in that feature's attributes using ArcGIS Desktop

texts:
I want to take a point feature and create copies of the geometries based on a value in the attribute table associated with that point.

Let's say you have 3 points; A, B and C.

Point A has a value of 2
Point B has a value of 5
Point C has a value of 0

The output from my desired result would have 7 points total.

2 records of Point A
5 records of Point B
0 records of Point C

I'm working in ArcGIS Desktop 10.7.1.

I have some experience in ArcPy as well if that would be a preferred approach.

Ideally I'd like to find a way to automate this; essentially selecting a feature and copying and pasting it by the number of times appropriate based on it's field value.

I have been thinking about utilizing search/update cursors, but not sure how efficient that would be.

For reference, there are about 32,000 points with a sum value of 76,000 in the desired field. So the final output should have 76,000 points.