texts:
I have a Virtual Table which is a result set for a particular process.

select
  "RECORD_IDENTIFIER",
  CAST ("UPRN" AS textstring) AS UPRN,
  "KEY",
  "WARNING_CODE",
  "EXCEPTION"
from Exceptions
where sent = 0

The whole process should be run by co-workers who have no programming knowledge but just need to run the process. The result of the query needs to be output to a text file .txt without field names with every field comma-separated and enclosed in double-quotes.

I have no idea how to do this last part i.e. export the virtual query to a .txt file. Any ideas?


texts: Iterating through feature classes in geodatabase using ModelBuilder

texts:
I’m trying to create a model so I can iterate through some feature classes in my gdb. I remember doing it one time while working on wetlands with the US states, where I used the state boundary shape file and the select by attributes tool.

How do I connect the Select by attributes tool to the Iterate Model to enable me select the feature classes I want to use a geoprocessing tool on?

I created a shapefile (utilities) with the list of selections with the names matching the feature classes in the gdb.


texts: How can I install Geomesa on Google Cloud using BigTable. I have data with lat , long and user identification. I need to create heat maps, path taken by most users, most frequented shops etc.I am looking for a solution on google cloud.


texts: Querying intersecting areas areas - Oracle Locator or Spatial?

texts:
My spatial database is Oracle. We are wondering if we can use basic Oracle Locator or need to get the more expensive, full, Oracle Spatial.

My application's requirements include having to be able to query "which geometries in the spatial db intersect this geometry X".

The whole locator / spatial issue is a bit confusing since it looks like Oracle initially installs all Oracle Spatial features even if you are just using basic Oracle db and do not have license for full Oracle Spatial. And you may end up using the Spatial features by mistake: http://www.spatialdbadvisor.com/oracle_spatial_tips_tricks/55/oracle-locator-vs-enterprise-licensing-the-sdogeom-package

What I understand of Oracle documentation http://docs.oracle.com/cd/E11882_01/appdev.112/e11830/sdo_locator.htm#SPATL340 is that geometry function SDO_GEOM.SDO_INTERSECTION is not supported by Oracle Locator.

So does this mean that we really need full Oracle Spatial license to fulfill the requirements? Or is there some other way then "SDO_GEOM.SDO_INTERSECTION" of finding about whether areas intersect?

We are using GeoServer on top of Oracle (and accessing the geometries with WFS-T). Does this change anything, will the intersection be somehow calculated in GeoServer (doubt it?) or will the intesection-calculation have to be done in Oracle (probably?)?


texts: Style by categories for items with no value in QGIS

texts:
I need to style a shapefile by graduated, but some items has no value for that items. (eg. the limit for each categories are 0/99 - 100/199 - and so on, but some items has no value, that means that the real value is zero).
I would like to represent all items with no value in the first step (0/99) but I do not want to create records with value "0" for those items.

How can I solve it?


texts: Special characters changing to "?" using ArcMap

texts:
I am using shapefiles in Spanish that I'm using to build a database and the special characters are displayed as normal,  but as soon as I ran an operation on the layer these characters change to ?.
I tried using field calculator to replace ? with them but as soon as I save the edits it reverts back.
I'm using ArcMap. When I added the shapefile the attribute tables had the special characters without a problem byt when I preformed an operation such as Merge or Dissolve or any other analysis that would deal with attributes, it would change every special charterer to ?. What I tried doing was using field calculator and replaced the ? with ñ for example and it worked but when I saved the edits the ? would return back as if I didn't replace them.


texts: Open source tools for creating and publishing animated web maps?

texts:
I am interested in creating animations from observations or model predictions and then displaying these animations in a web-based viewer such as OpenLayers.

The best example I have found of this functionality is at:

GODIVA2

In my opinion, the user interface is very unintuitive.  The steps required to get an animated map are:

Select a layer from the tree browser on the left.
Click on a day in the calendar on the right.
Click on the link labeled "first frame" which is to the left of the calendar.
Choose another day from the calendar and click on the link labeled "last frame".
A "create animation" link will pop up- click on that to generate the animation.

What open source toolkits are available with support for visualizing and animating time series?  On the front-end I am looking to display the results in a Google Maps-style viewer as that is what most people are familiar with.  The back-end is wide open.


texts:
I don't have much experience with QGIS and am working on a map with lines connecting locations created with the mmqgis plugin. Is there an easy way to make those direct lines somewhat curvy for aesthetic purposes?

I'm trying to achieve something like the map below but with several hubs. I have read suggestions on how to do this with PostGIS and other advanced tools, so I've done my research (including other posts here) and tried using several plugins, and I'm specifically asking for an easy solution because I don't have much experience.


texts: Couldn't load plugin sextantelwgeomprovider, QGIS 2.0.1 on Ubuntu 14.04.2

texts:
I'm getting this error message when I start QGIS

  Couldn't load plugin sextantelwgeomprovider due an error when calling its classFactory() method

  Traceback (most recent call last):
    File "/usr/lib/python2.7/dist-packages/qgis/utils.py", line 204, in startPlugin
      plugins[packageName] = package.classFactory(iface)
    File "/home/sara/.qgis2/python/plugins/sextantelwgeomprovider/init.py", line 50, in classFactory
      from sextantelwgeomprovider.SextanteLwgeomProviderPlugin import SextanteLwgeomProviderPlugin
    File "/usr/lib/python2.7/dist-packages/qgis/utils.py", line 453, in _import
      mod = _builtin_import(name, globals, locals, fromlist, level)
    File "/home/sara/.qgis2/python/plugins/sextantelwgeomprovider/SextanteLwgeomProviderPlugin.py", line 28, in
      from sextante.core.Sextante import Sextante
    File "/usr/lib/python2.7/dist-packages/qgis/utils.py", line 453, in _import
      mod = _builtin_import(name, globals, locals, fromlist, level)
  ImportError: No module named sextante.core.Sextante
  Python version:
  2.7.6 (default, Jun 22 2015, 18:01:27)
  [GCC 4.8.2]

QGIS version:

  2.0.1-Dufour Dufour, exported

Python path:

  ['/usr/share/qgis/python', '/home/sara/.qgis2/python', '/home/sara/.qgis2/python/plugins', '/usr/share/qgis/python/plugins', '/usr/lib/python2.7', '/usr/lib/python2.7/plat-x86_64-linux-gnu', '/usr/lib/python2.7/lib-tk', '/usr/lib/python2.7/lib-old', '/usr/lib/python2.7/lib-dynload', '/usr/local/lib/python2.7/dist-packages', '/usr/lib/python2.7/dist-packages', '/usr/lib/python2.7/dist-packages/PILcompat', '/usr/lib/python2.7/dist-packages/gtk-2.0', '/usr/lib/pymodules/python2.7', '/usr/lib/python2.7/dist-packages/ubuntu-sso-client', '/usr/lib/python2.7/dist-packages/wx-2.8-gtk2-unicode']

Because of that I can't use any geoalgorithms. Is there a way of solving the bug: "No module named sextante.core.Sextante"?

I tried to fix it by using Plugin > Manage and Install Plugins... but I still receive and error message "The plugin requires a missing module sextante.core.Sextante". How can I install this module?


texts: I'm new using ArcMap and I have a GDB of schools (and its enrollment data) at block level. What I need to do is to find the places with the most concentration of enrollment in a radius of 5km. Since I'm new using this tool I'm completely lost.


texts:
I will preface this by saying I am new to this field. I am using a Trimble geo7x with a zephyr 3 to collect elevation points, but my data are not what I expected to see: some are negative and do not come close when I compare uncorrected or corrected data to DEM data in gis.

I used com3 for my port setting in the field, and I’m thinking this is my error. I’m set for msl, geoid 12a, and I have the appropriate antenna height.

My question is: is there any way to correct this in office? Running a differential correction in pathfinder does not get me where I need to be.


texts: Iterating through all geometries using ArcPy search cursor?

texts:
I would like to use a python add-in which uses a button click to iterate through each selected polygon one by one.

The following code manages to zoom only to the first polygon. How can I adapt it to zoom to each selected polygon?

mxd = arcpy.mapping.MapDocument('CURRENT')
df = arcpy.mapping.ListDataFrames(mxd, "Layers") [0]
Envelopes = [] # store extents here

# find the selection set
SelLayer = arcpy.mapping.ListLayers(mxd,data_frame=df)[0] # first layer
fidSet   = arcpy.Describe(SelLayer).FIDSet

if len(fidSet) == 0:
    arcpy.AddMessage("Nothing selected")
else:
    # now cursor through it an get the geometries
    # storing their extents into a list
    with arcpy.da.SearchCursor(SelLayer,"SHAPE@") as SCur:
        for feat in SCur:
            # I'm going to assume polygon/polyline
            Envelopes.append(feat[0].extent) # grab the envelope

    df.extent = Envelopes[0] # first extent
    arcpy.RefreshActiveView()


texts:
I have two geodatabases and I import a large number of tables from one to another by using ArcCatalog (ArcGIS 10.2). However the tables and indexes created on them are imported with different srid.  Lets say, for the same table in geodatabase-1 srid that I am using is 6 but in geodatabase-2 it is 30002. I need to set srid for all tables to 6. The reason why I need same srid's is that two geodatabases must act like clones to each other.

Can anyone help me?


texts: Making map with gradation polygon using QGIS

texts:
I am making a history map on my website these days. (http://eastmaps.com/2019/09/10/%ed%85%8c%ec%8a%a4%ed%8a%b8/ < it's my test version)

But like this site(https://www.ancient.eu/map/), I want to give gradation to the lines outside the polygon on the map.

I am using QGIS.

Do you know how to give a polygon gradation?


texts: Spplot SpatialPolygonDataFrame and SpatialGridDataFrame

texts:
I have to plot in the same graph something like this:

Then, I have elements from SpatialPolygonDataFrame class and others elements from SpatialGridDataFrame class. When I try to plot them in the same graph only get this error:

Error in `[[<-.data.frame`(`*tmp*`, i, value = c(52.6195687909829, 52.6845097109409,  :
  replacement has 64 rows, data has 672.

However, I converted all elements to SpatialPolygonDataFrame class, but it does not work.

I've watched the data and I think that the problem is the size of the elements (row, columns) but I'm not sure. I do not how to do it.

This is a example of my routine.

# Interpolation kriging
k_066_1617H <- krige(data066_1617H$Vtec~1,locations = data066_1617H,newdata = gridpolygon2,model = Ajuste2)

class(k_066_1617H)

#[1] "SpatialPolygonsDataFrame"
#attr(,"package")
#[1] "sp"

Then, read a TIFF file

land <- readGDAL("image.tif")
class(land)

#[1] "SpatialGridDataFrame"
#attr(,"package")
#[1] "sp"

land <- Grid2Polygons(land) # Transform into SpatialPolygonDataFrame

class(land)

#[1] "SpatialPolygonsDataFrame"
#attr(,"package")
#[1] "sp"

Construct a vector with land and k_066_1617H for plot with spplot

zn.k_066_1516H = land

zn.k_066_1617H = k_066_1617H
zn = land
zn[["a"]] = land[["z"]]
zn[["b"]] = zn.k_066_1617H[["var1.pred"]]

Here, is the problem. Maybe the size of elements is the problem.

land has 28 rows and 24 columns and k_066_1617H has 3752 rows.

spplot(zn, c("a","b"),
            main = "Santander, VTEC día 066",
            #names.attr = c("Prueba","18 - 19", "16 - 17"),
            lty=0,
            lwd=7,
            sp.layout = c('sp.points', santos.coord, col='red', pch=19),
            contour=TRUE,
            col='black',
            pretty=TRUE,
            scales=list(draw = TRUE),
            labels=TRUE)


texts: change Symbol color on click, ArcGIS runtime 100 .net

texts:
I have a simple GraphicsOverlay symbol, how do I change the color of it on event triggered?
My sample code shown below:

    public MainPageViewModel()
    {
        // Create new Map without basemap
        Map = new Map(SpatialReference.Create(102100));

        // Add the tile basemap only after the map has loaded so it doesn't block the map if offline
        Map.Loaded += (sender, e) => { Map.Basemap = Basemap.CreateImagery(); };

        // Setup the graphics layers
        GraphicsOverlay editGraphics = new GraphicsOverlay();
        editGraphics.Id = "MyGraphics";
        GraphicsOverlays.Add(editGraphics);

        // Add some test graphics

        SimpleMarkerSymbol pointSymbol = new SimpleMarkerSymbol(SimpleMarkerSymbolStyle.Circle, SystemColor.Red, 25);
        var point = new MapPoint(0, 0);
        Graphic mapPoint = new Graphic(point, pointSymbol);
        GraphicsOverlays["MyGraphics"].Graphics.Add(mapPoint);
    }

    public ICommand ColorButtonCommand { get; set; }
    public void ColorButtonTapped()
    {
        // Handle color code goes in here
        // pointSymbol.Color = SystemColor.Gray;

    }
    public bool CanColorButtonTapped()
    {
        return true;
    }

I thought of declaring my SimpleMarkerSymbol pointSymbol variable as a global variable, and simply change color like what I did in my commented out code pointSymbol.Color = SystemColor.Gray, But I think it is a bad approach, there should be a easy and clean way to handle it.


texts: Multiplying Raster Value with Pixel size using QGIS?

texts:
i'm working on a map with population densities. My raster data contains values as [People/square meter] for each cell. I have per capita values of energy demand etc. which I want to multiply with the raster to get a spatial distribution of the demand, but this requires the number of people in one cell. Consequently i'd have to multiply the population density values with the area of each cell, but here's the problem:

Each cell is 0.0083 degrees but depending on the distance from the equator has a different area.

How can I extract the area values from this raster/ directly calculate population numbers of the cells?

My first try was to vectorize the raster, calculate the area for each polygon and rasterize again, but vectorizing it, QGIS consolidates cells with similar values to one polygon. From 660,000 raster cells I get 60,000 polygons.


texts: polygon dissolve or remove duplicate records with smaller geometry

texts:
I have a polygon table in postgres that has important columns of: address,type_(i know confusing wasnt me),subtype and geom.

this table is parcels buffered 1000 feet. some records like the picture below contains both the parcel geometry and buffer geometry with the same exact records address,type_,subtype. I use this layer to perform various distance calculations to other parcels by the b_uniqueid and the records with just the parcel geometries are giving me a boat load of trouble.

things I have tried:

drop table if exists test2;
create table test2 as
select st_union(geom) geom,address,type_,subtype from buffer_unions group by address,type_,subtype

gives me a

ERROR:  GEOSUnaryUnion: TopologyException: Input geom 1 is invalid: Self-intersection at or near point 20021031.841838982 -322552.71008741739 at 20021031.841838982 -322552.71008741739
********** Error **********

ERROR: GEOSUnaryUnion: TopologyException: Input geom 1 is invalid: Self-intersection at or near point 20021031.841838982 -322552.71008741739 at 20021031.841838982 -322552.71008741739
SQL state: XX000

and in arcgis the dissolve tool does not give me the correct results.

I have also tried:

select st_multi(st_collectionextract(st_collect(geom),3))::geometry(multipolygon,102689) geom,address,type_,subtype from buffer_unions group by address,type_,subtype

but because the 700 records of subtype OutdoorRecreationFacility has no address field those columns get dissolved down to one record.

What I am looking for is a query which goes through the table and analyzes
1. if a smaller geometry is inside a larger geometry and has the same address and subtype, dissolve them or alternatively maybe if have the same address and subtype compare there st_area's and just take the larger one (not sure how to do that)
2. for the OutdoorRecreationFacility records, because they have no addresses I need to either dissolve if a geometry falls completely within another geometry or just select the larger one

I tried this query but got the the topology error like above

select case
when st_contains(t1.geom,t2.geom) then st_union(t1.geom,t2.geom)
else t1.geom end geom
from buffer_unions t1 cross join buffer_unions t2
where t1.subtype = 'OutdoorRecreationalFacility' and t2.subtype  ='OutdoorRecreationalFacility'


texts:
I am using GRASS GIS 7.2.0. The intent is to clip my DEM raster. I have searched available raster tools but didn't found anything allow me to clip my raster using polygon vector or interactively.

How would I clip my raster in GRASS GIS? Is there any plugin or tool to do the job?


texts:
I have two coordinates.

Lat:53.307718   Long:-6.230189

Lat:53.307708   Long:-6.229988

.

I then use the distance formula

distance = SQRT((Y2-Y1)*(Y2-Y1) +(X2-X1)*(X2-X1))

.

I put it points in.

distance = SQRT((53.307708-53.307718)*(53.307708-53.307718) +(-6.229988--6.230189)*(-6.229988--6.230189))

and it get a value of 0.000201

What is the factor to convert that point to convert to meters?

I but the two points into this website http://www.movable-type.co.uk/scripts/latlong.html and found the distance to be 12.92 meters. But when i do it then for the next point the factor isnt the same.


texts:
I have a quite large high-resolution map in a format called .rik from the Swedish land survey agency.
I want that whole high resolution file in different format that can be used to print a map from a large area.

When i choose "Save as picture" in the Archive menu now, it only exports the view seen on the screen.
If I zoom-out to see the full extent of the map and then export it the resolution is poor.

I want to export the whole file in the highest possible resolution.

Is there a way to do this?


texts: Deliminated layer cuts off imported spreadsheet at row 20

texts:
Platform: Mac, OS: Mavericks, Qgis: 2.14 Essen

Exporting spread sheet from Excel as Windows .csv. No problem publishing points to map except that import will only process 20 of 23 row entries on spreadsheet. Spreadsheet is consistent. What is the resolution for this?  I can work around by splitting the spread sheet into two separate files but that seems a bit much for something I'm certain I must be missing, a simple oversight?.


texts: Rasterize Vector MBTiles - vector mbtiles to raster mbtiles?

texts:
I have a mbtiles which is a vector format; when I opened this in QGIS, dont see any symbology associated, but with Maptiler I can see all styles correspondingly. But my problem, I need to rasterize or convert vector mbtiles to raster mbtiles?
I did tried using following github forum https://github.com/bjornreppen/rasterize-vector-tiles/issues/6 but no luck to get the exact result!
Styled mbtiles from MapTiler

Non-style mbtiles in QGIS


texts:
everyone.

I created a table in oracle with a sde.st_geometry column. I followed the instructions of ArcGIS help to register the table into geodatabase. However, when I try to load data into this table, I got this error:

ora-01747: user.table.column, table.column or column description is invalid

My table creation statement is:

create table tbl_Comm_Cove
(
ID                   NUMBER(10,0)         not null,
Name                 CHAR(200),
Devi_Num             CHAR(8),
Code                 CHAR(200),
Unit_Grid            CHAR(10),
GID                  CHAR(200),
Belo_Road            CHAR(200),
Inst_Time            DATE,
Mana_Depa_ID         NUMBER(10,0),
Owne_Depa            CHAR(200),
Main_Depa            CHAR(200),
Mate_ID              NUMBER(10,0),
Diam                 NUMBER(10,2),
Devi_Stat_ID         NUMBER(10,0),
"Desc"               CHAR(200),
Shap_ID              NUMBER(10,0),
Init_Time            DATE,
Chan_Time            DATE,
Note                 CLOB,
constraint PK_TBL_COMM_COVE primary key (ID)
);
/
alter table tbl_comm_cove add loc sde.st_geometry;
/
CREATE INDEX tbl_comm_cove_idx ON tbl_comm_cove(loc)
INDEXTYPE IS sde.st_spatial_index
PARAMETERS('st_grids=1,3,0 st_srid=0');
/

The id is self—incremental.
My oracle version is 11.2.0, and ArcSDE version is 10.0

What should I do?


texts: Automatically reducing polygon size in QGIS?

texts:
I am working with two polygon layers (counties & states) which I will 'intersect' to find which counties are within which states.
Due to data inaccuracies I have to change all county polygons and make them a bit smaller. At this stage it doesn't really matter by how much, as long as I can do this change automatically.

I understand that polygons which are currently perfectly aligned now, will afterwards have a 'gap' between them, but that's fine.

I've tried 'Simpliyfy Geometries' but this only 'smoothens' edges and doesn't help in this case.

Is there any option in QGIS to reduce the size of polygons automatically? (i.e. move all polygon points 200 meters closer to the 'polygon mid-point').


texts:
A client of ours is having trouble with some raster tiles in ArcMap. They gave me a sample, and I was able to get it looking right with the following workflow:

All you have to do is go into each of your loaded raster images, right-click, Properties, Symbology tab. Under Stretch, change the drop-down to None, and unselect Apply Gamma Stretch. Then it should look normal.

They tried this, and got it to work on partial tiles, but not entire ones. Is there any better solution? Attached is a tile sample.  I'm using ArcGIS 10.2.2 for Desktop.


texts: Get feature information from WFS layer in OpenLayers 3

texts:
In case of OpenLayers 2 I am able to get feature information but how to get feature information from vector or WFS layer in OpenLayers 3

Below code how to extract information of the vector layer if multiple vectors are there?

  var layerWFS = new ol.layer.Vector({
  source: new ol.source.Vector({
    loader: function(extent) {
      $.ajax('http://localhost:8080/geoserver/wfs', {
        type: 'GET',
        data: {
          service: 'WFS',
          version: '1.1.0',
          request: 'GetFeature',
          typename: 'dgm:all_block_boundary_point',
          srsname: 'EPSG:3857',
          bbox: extent.join(',') + ',EPSG:3857'
        }
      }).done(function(response) {
        layerWFS
        .getSource()
        .addFeatures(new ol.format.WFS()
          .readFeatures(response));
         // console.log(response);
      });
    },
    strategy: ol.loadingstrategy.bbox,
    projection: 'EPSG:3857'
  })
});
map.addLayer(layerWFS);


texts:
I want to calculate the RMSE of two grids based on their z-value: a calculated and a reference grid.
I do have the statistics: Cells, Range, Min, Max, Sum, Mean, Variance, St.dev.

How could I solve the problem using e.g. ArcGIS or SAGA Gis?


texts: How to display related records in popup using ArcGIS Javascript API?

texts:
I am attempting to display related fields in a popup window.  I've created a GraphicsLayer, consisting of features from the results of a QueryTask.  Upon adding the features to the GraphicsLayer, I am assigning a popupTemplate to each feature:

var popupTemplate = {
  title: "{LONGNAME}",
  content: [{
    type: "fields",
    fieldInfos: [{
      fieldName: "SPACEID",
      label: "Space Identifier",
      visible: true
    },{
      fieldName: "relationships/1/VERTICALORDER",
      label: "Vertical Order",
      visible: true
    }]
  }]
};

However, this technique does not seem to be able to reference the related table. The [SPACEID] field value is visible, but the related field value is not. Am I doing something wrong, or is there a different way to accomplish this?

EDIT: Perhaps this functionality is not possible yet?  I've tested this with a simple FeatureLayer, and the relationship is not honored.

var layer = new FeatureLayer({
    url: baseURL + "/2",
    popupTemplate: {  // autocasts as new PopupTemplate()
        title: "{LONGNAME}",
        content: [{
            type: "fields",
            fieldInfos: [{
                fieldName: "SPACEID",
                label: "Space Identifier",
                visible: true
            },{
                fieldName: "relationships/3/FIRSTNAME",
                label: "Employee Name",
                visible: true
            }]
        }]
    }
});

I have verified the relationship is functioning, using the REST queryRelatedRecords window.  However the FeatureLayer does not see the table.  Is there something further that needs to be done?  I should note that the service is a MapService, in case that matters.

This issue is discussed here and elsewhere.


texts: Inserting rectangle in order to surround certain area (few raster cells) using ArcGIS for Desktop?

texts:
I´ve done a suitability analysis and now I need to point to this certain area.

Is it possible to select and frame certain raster cells to enlighten those?
Or can somebody explain how to insert a simple rectangle that cover an area somewhat precise?


texts: Load a GeoJSON layer with QGIS's Python module

texts:
I have been able to load a GeoJSON file into QGIS by simply clicking "Add Vector Layer", and selecting my file to load.

Now I'm using the qgis.core Python package, and I'm trying to load the GeoJSON file as a layer. I'm following instructions from the QGIS Python Cookbook - Loading Layers section; but I can't find how to specify that I'm loading a GeoJSON file.

I've tried:

QgsVectorLayer("my_geo.json","mygeojson","geojson")

and

QgsVectorLayer("my_geo.json","mygeojson","json")

to no avail. Any advice on what I could do?


texts: Sum of attributes of a subset of the features

texts:
This is my attribute table of a line layer:

How can I get the length of all segments where Rohr_1 equals SNRVe 7x12x2 ?


texts:
I currently have a DEM raster stored in an ESRI geodatabase. I've been asked to create a purely python (e.g., no arcpy functions) program that works with the values. Therefore, I need to convert the raster into something that can be read by python as an array. What file format would that be?

(Alternatively, is this a terrible idea and should I argue for sticking with arcpy functions?)


texts: Tilemill/PostgreSQL Connection refused Error

texts:
I have been trying to set up a MBTiles Export using PostgreSQL 9.2 and the Latest version of Tilemill. According to the postgreSQL Logs, I am receiving an error, and I cannot figure out why. I have them both setup on the same machine (localhost:5432).

Here is the error below (Connection Refused, and Client lost connection). Sometimes it loads into TileMill, but very intermittent and slowly. It doesn't allow me to zoom or pick the BBox I want to export, It never refreshes.

    2015-07-28 16:17:47 EDT STATEMENT:  SELECT ST_AsBinary("way") AS geom,"way_pixels" FROM (SELECT
        way,
        name,
        boundary,
        way_area/NULLIF(1222.99::real*1222.99::real,0) AS way_pixels
      FROM planet_osm_polygon
      WHERE (boundary = 'national_park' OR leisure = 'nature_reserve')
        AND building IS NULL
    ) AS national_park_boundaries WHERE "way" && ST_SetSRID('BOX3D(-9549125.069610503 4226661.916057105,-8609866.866042253 5165920.119625353)'::box3d, 900913)
2015-07-28 16:18:14 EDT LOG:  could not send data to client: No connection could be made because the target machine actively refused it.

2015-07-28 16:18:14 EDT STATEMENT:  SELECT ST_AsBinary("way") AS geom,"landuse","natural","waterway","way_pixels" FROM (SELECT
        way, "natural", waterway, landuse, name, way_area/NULLIF(2445.98::real*2445.98::real,0) AS way_pixels
      FROM planet_osm_polygon
      WHERE
        (waterway IN ('dock', 'riverbank', 'canal')
          OR landuse IN ('reservoir', 'basin')
          OR "natural" IN ('water', 'glacier', 'mud'))
        AND building IS NULL
        AND way_area/NULLIF(2445.98::real*2445.98::real,0) > 0.01
      ORDER BY z_order, way_area DESC
    ) AS water_areas WHERE "way" && ST_SetSRID('BOX3D(-9079495.967826378 3443946.746416901,-7200979.560689883 5322463.153553395)'::box3d, 900913)
2015-07-28 16:18:14 EDT FATAL:  connection to client lost


texts: How to Load SQL Geometry Data on OpenLayers3?

texts:
I am New in OpenLayers 3. I have to load Geometry data from sql to my OpenLayers map.

This is my code :

var tileLayer = new ol.layer.Tile({
source: new ol.source.MapQuest({ layer: 'osm' })
})

var vectorLayer = new ol.layer.Vector({
source: vectorSource,
style: new ol.style.Style({
    fill: new ol.style.Fill({
        color: '#0099FF'
    }),
    stroke: new ol.style.Stroke({
        color: '#0033CC',
        width: 2
    }),
    image: new ol.style.Circle({
        radius: 7,
        fill: new ol.style.Fill({
            color: '#0099FF'
        }),
        stroke: new ol.style.Stroke({
            color: '#0033CC',
            width: 2
            })
        })
    }),
});

map = new ol.Map({
target: 'map',  // The DOM element that will contains the map
renderer: 'canvas', // Force the renderer to be used
layers: [tileLayer, vectorLayer],
overlays: [overlay],
view: new ol.View({
    center: ol.proj.transform([2, 41], 'EPSG:4326', 'EPSG:3857'),
    zoom: 3
    })
});

This my code for a OpenLayers3 map.
I have data with geometry field in my table. So, i want to draw on my openlayers map with this geometry fied.

This is Data in my geometry Field.
0xE6100000010CDDB5847CD0FF56C00EBE30992A583D40

so, how can i draw with this type of data?


texts: Iteration over features in PyQGIS processing modules

texts:
I have an issue with processing.run commands where I would like to do an iteration over features in a layer with polygons.
zsj = QgsVectorLayer("/data/grassmaster/testVlado/zsj_ba_test.gpkg", "zsj", "or")

si = QgsProject.instance().mapLayersByName("Vectorized")[0]

it = zsj.getFeatures()

for i in it:

    clip = processing.run("qgis:clip", {'INPUT':si, 'OVERLAY':it[i], 'OUTPUT':'memory:'})
    print(clip['OUTPUT'])

Here is the error:
TypeError: 'QgsFeatureIterator' object is not subscriptable

How can I use this processing command for individual features in the layer? Respectively, how can I save a feature in the layer as a new (memory) layer?


texts: Setting the_geom to coordinate system EPSG5514 in CARTO?

texts:
How can I set my geometry column to coordinate system EPSG5514 Krovak in CARTO? I would like to publish maps on this system.

If I enter:

UPDATE
   komunikace_jtsk
SET
   the_geom = ST_Transform (ST_SetSRID (the_geom, 5514), 4326);

So the data will be displayed outside of the Czech Republic somewhere on the North Pole.
Also the data in the_geom layer contains geographic coordinates and not equal - EPSG5514.
the_geom:

{"type":"MultiPolygon","coordinates":[[[[24.83035935,59.75626015],[24.83035935,59.75626015],[24.83035935,59.75626015],[24.83035935,59.75626015],[24.83035935,59.75626015]]]]}


texts: Is it scientifically correct to pansharpen landsat reflectance product with pan band

texts:
I need to know if it is scientifically valid to fuse/pan-sharpen Landsat 8 surface reflectance products with pan band of that respective band? Landsat reflectance product details can be found at here. It needs to be mentioned that one needs to order surface reflectance product separately to get this product. This product contains only 7 band (30m) not IR and Pan band. So, again, my question is it valid to fuse 7 bands(30m) of surface reflectance product with normal(not surface reflectance) pan band(15m).
I want to use this pansharpened image for segmentation and following land cover mapping.So I need to know that is there any established practice of this type of pan sharpening  in academia with reference, if yes please cite.


texts: How can I map data points to OpenStreetMap (OSM) line data?

texts:
my goal is to filter measuring points of vehicles by mapping them to streets.
So I could filter out all the false data located in fields and rivers and so on.

I imported my data into a PostGIS database and can work fine with my Point data.
I also imported the part of OSM that i need into a PostGIS table.

Is using ST_DWithin the right way? I did use that to only show points in a certain distance to another point.

How can I now get only these points that are located, let's say 20 metres around the OSM lines?

EDIT1:

I have been trying to first filter the data and the OSM line data to a circle area of 5km and then using ST_DWithin to get the mapping:

WITH
    Data as
    ( SELECT gid,geom FROM schema1.data AS Mp
        WHERE Mp.id=14
        AND EXISTS (
            SELECT 1 FROM schema1.base As Base
            WHERE Base."ID"=14
            AND st_covers(st_makeenvelope(-180, -90, 180, 90, 4326), Mp.geom)
            AND ST_DWithin(Mp.geom::geography, Base.geom::geography, 5000)
            )
    ),
    Line As
    ( SELECT * FROM public.planet_osm_line As Line
        WHERE EXISTS (
            SELECT 1 FROM schema1.base As Base
            WHERE Base."ID"=14
            AND ST_DWithin(Line.way::geography, Base.geom::geography, 5000)
            AND Line.highway='motorway'
            )
    )
SELECT
    Data.*
FROM
    Data, Line
WHERE
    st_covers(st_makeenvelope(-180, -90, 180, 90, 4326), Data.geom)
    AND ST_DWithin(Data.geom::geography, Line.way::geography, 20);

Unfortunately "Data" still consists of 5014309 and "Line" of  5686 rows.
What can I improve in my query to speed up things here?

If i understand the query execution correctly, postresql is trying to build a temporary table and thus joining Data and Line to 5014309*5686 rows...

I would appreciate any help.


texts: Ionic Leaflet delay/latency retrieving images

texts:
I have a very basic mapping application using Ionic 4 / Leaflet and functionally it is working. However I notice a significant delay on loading basemap tiles when in app compared to in my desktop dev environment. The images load almost instantly on desktop versus there being noticable delay in app.

The exact same embedded image e.g.

Takes 55ms on a desktop browser vs 420ms in app.

In Dev/Desktop

In App (taken with chromes remote devices dev tool)

What explains the 230ms “connection start” in the app? Is this related to IONIC or the webview in any way and is there a way to circumvent this? It is very noticeable and makes the performance of the app look much worse than it actually is.

The app itself is just a very basic Leaflet application

this.map = Leaflet.map("map", {
    center: [this.latitude, this.longitude],
    zoom: 6,
    minZoom: 3,
    maxZoom: 17,
    inertiaDeceleration: 1000,
    attributionControl: false,
    zoomControl: false,
    maxBounds: Leaflet.latLngBounds([-90, -180],[90, 180])
});

const basemap = Leaflet.tileLayer('https://{s}.tile.opentopomap.org/{z}/{x}/{y}.png',  {
    maxZoom: 17,
    minZoom: 0
});

basemap.addTo(this.map);

Using "leaflet": "^1.5.1"

Ionic:

Ionic CLI : 5.4.2 (/usr/local/lib/node_modules/ionic)
Ionic Framework : @ionic/angular 4.7.4
@angular-devkit/build-angular : 0.801.3
@angular-devkit/schematics : 8.1.3
@angular/cli : 8.1.3
@ionic/angular-toolkit : 2.0.0

Cordova:

Cordova CLI : 9.0.0 (cordova-lib@9.0.1)
Cordova Platforms : android 8.1.0
Cordova Plugins : cordova-plugin-ionic-keyboard 2.2.0, cordova-plugin-ionic-webview 4.1.2, (and 6 other plugins)

Utility:

cordova-res : not installed
native-run : 0.2.8

System:

Android SDK Tools : 26.1.1 (/home/anthony/Android/Sdk)
NodeJS : v12.11.1 (/home/anthony/.nvm/versions/node/v12.11.1/bin/node)
npm : 6.11.3
OS : Linux 5.0


texts: About spatialite views and Attribute table

texts:
I try to use views with spatial column in spatialite. I can connect to spatialite database and load the layer without problems but when I try to view data in Attribute Table I can't get eny data; all rows and field show ERROR.

Data render well and information in point show correct data.

Has someone get the same error?.
Is this a bug?

I get this results in different versions of qgis (osg4w in windows, master version in ubuntu 12.04 and qgis version included in osgeolive 6.0)


texts: Splitting river channel polygon by fishnet overlaying it using ArcGIS Desktop?

texts:
I have to split river channel polygon by fishnet in ArcGIS Desktop.  There are 500 squares in fishnet of 1km sq area which overlay above river channel polygon.

Is there any tool in ArcGIS Desktop by which I can directly split the river channel polygon into 500 squares?


texts:
I want to download the feature layer which I input as a string parameter to my GP service. When I run the script I got the following errors.

First time when I run it:

job-submitted  gp_service_run2.php:110
job-failed

When I reload the page:

MapView.js:381 [Violation] Added non-passive event listener to a scroll-blocking 'wheel' event. Consider marking event handler as  'passive' to make the page more responsive. See https://www.chromestatus.com/feature/5745543795965952
(anonymous) @ MapView.js:381 set  @ MapView.js:380 m.updateDependencies @ MapView.js:436 updateDependencies @ MapView.js:432 onResume @ MapView.js:349
b.resume @ MapView.js:351 (anonymous) @ MapView.js:362 m @ dojo.js:1079
p.onInstall @ MapView.js:362 b @ MapView.js:432 m.installHandlers @ MapView.js:433 g._connect @ MapView.js:342 (anonymous) @ MapView.js:341
(anonymous) @ dojo.js:1077 a @ dojo.js:363 m @ dojo.js:364 p @ dojo.js:370
(anonymous) @ dojo.js:374 gp_service_run2.php:110

job-submitted gp_service_run2.php:110
job-executinggp_service_run2.php:110
job-failed

In all other occasions when reloading same error and:

job-submitted gp_service_run2.php:110
job-failed

Where I have made the mistake in my code?

<script>
    require([
        "esri/Map",
        "esri/views/MapView",
        "esri/layers/FeatureLayer",
        "esri/tasks/Geoprocessor",
        //"esri/tasks/FeatureSet",
        "esri/views/MapView",
        "esri/layers/support/ImageParameters"
      ],
      function(
        Map,MapView, FeatureLayer,Geoprocessor,ImageParameters,FeatureSet ) {
        var gpUrl="http://www.xxxxxxxxxxxxxxxxxxxxxxx/GPServer/Script";
        var map = new Map({
          basemap: "gray"
        });

        var view = new MapView({
          container: "viewDiv",
          map: map,
          center: [80.9,7.8],
          zoom: 8
        });
        var featureLayer2 = new FeatureLayer({
             url:"http://www.xxxxxxxxxxxxxxxxxxxxxxxxxxxx/FeatureServer",
          outFields:["*"],
            popupTemplate:{
          title:"Parcel Deatils",
         content:"Lot No:"+"{sht_pg_id}"
          }
            });

        map.add(featureLayer2);
        var gp =new Geoprocessor(gpUrl);
        //gp.outSpatialReference = map.spatialReference;
       gp.outSpatialReference = {
        wkid: 5235
       };

        var params = {
         "sheet": featureLayer2,
            "Plan_No":"28006103S001"
             };
        gp.submitJob(params).then(drawResultData, errBack, progTest);

         function drawResultData(result) {
         //var resultLayer = gp.getResultData(result.jobId);
          //map.add(resultLayer);
          //document.getElementById("demo").innerHTML = ""+result.jobId;

          if (result.jobStatus == "job-succeeded") {
                    gp.getResultData(result.jobId,"save",downloadFile);
                      }     };

        function progTest(value) {
          //message.innerText = "Job status: " + "'" + value.jobStatus + "'";
          console.log(value.jobStatus);
        }
        function errBack(error) {
          //message.innerText = "";
          console.log("gp error: ", error);
        }

        });

  </script>


texts: Use Python Fuzzy string matching library to match string between a list with 30 value and a given value and get the closest match in field calculator

texts:
I am trying to get the closest match between different strings but
The Function is limited to only 5 values to compare and I want to increase it to 30 value at least, any ideas how to work around this from ArcMap field calculator?

from fuzzywuzzy import fuzz
from fuzzywuzzy import process
def SetMatchPercenage(query,v1,v2,v3,v4):
    choices = [v1,v2,v3,v4]
    result= process.extractOne(query, choices)
    return result

The output is the top match from the given values like this

  ('Barack H Obama', 95)


texts: Define "geospatial data" for a non-GIS professional

texts:
I realize this is a sort of a subjective question, but I've seen other subjective questions on the site, so I thought it I'd put it out there.
A colleague recently had a request for "geospatial data" with a grant that they were submitting, and were confused by exactly what was being asked for.  My explanation was that they were asking for any GIS data related to the project area. My explanation seemed to work, but I wanted to put it to the community at large for a more defined definition.

How would you define "geospatial data" to a non-GIS professional?


texts: How can I convert .tfw/.tif files for use in ArcGIS?

texts:
I started a project in QGIS using .tfw files, and was able to simply open and manipulate the images by running them through the program. I tried to restart that same process, except with ArcGIS, but didn't expect the files to be unreadable. ArcGIS doesn't recognize these files as compatible, but I haven't been able to find a conversion program .tfw > .mxd or similar.

Any suggestions?


texts: How to change projection and show latitude / longitude in QGIS?

texts:
I am new to QGIS and pretty ignorant with the subject, but I want to be able to see the  latitude / longitude of a point (with a mouse click) in Quantum GIS.
I imported a shapefile (roads of finland from here : http://www.maanmittauslaitos.fi/en/node/6567 ), and there is a box a the bottom that show the current coordinate in unknow units to me : 521636,6947168, same the coordinate capture don't show lat/long. The coordinate system seems to be WGS84 - 4326.


texts:
I use Geoserver 2.8.1, Openlayers 2.12, GeoExt 1.1 on a Windows machine. I can successfully print a map using the Print Provider below:

var printProvider = new GeoExt.data.PrintProvider({
        method: "GET", // "POST" recommended for production use
        capabilities: printCapabilities,
        customParams: {
            mapTitle: mapTitle,
            comment: mapNote
        }
    });

In the previous example the "printCapabilities" variable points to this location and everything works fine;

http://99.99.99.99:8080/geoserver/pdf/info.json?var=printCapabilities

I have large complicated maps with many layers and using the "GET" option my print map requests frequently exceed the maximum request length and return an error.

I know I can change "GET" to "POST" in order to send longer map print requests. If I change my Print Provider to the following;

   var printProvider = new GeoExt.data.PrintProvider({
        url: "/geoserver/pdf",
        autoLoad: true,
        method: 'POST',
        customParams: {
            mapTitle: mapTitle,
            comment: mapNote
        }
    });

my print fails and FireBug gives me the following error explanation;

"XMLHttpRequest cannot load http://99.99.99.99:8080/geoserver/pdf/info.json. Response to preflight request doesn't pass access control check: No 'Access-Control-Allow-Origin' header is present on the requested resource. Origin 'http://99.99.99.99' is therefore not allowed access."

I already have an OpenLayers proxy setup as follows:

OpenLayers.ProxyHost = "/proxy.ashx?url=";

and I use the "URL-Rewrite" tool in Microsofts Internet Information Services (IIS) manager to redirect any URL's and avoid the same origin problem. It works for everything else but for some reason the POST option makes the server think that I am violating the same origin policy.

What do I have to do to make my system recognize that the print request is coming from the same server as everything else?


texts: I need to load a BAG file (Bathymetry attributed grid) in QGIS (3.4) according to the supported file page https://www.gdal.org/formats_list.html it is supported, but is not compiled by default and requires "libhdf5" I have no idea what this means, is there a plugin I need to install or do I need to compile something my self? I have no programming experience, so compiling software myself is not an option.


texts: How exactly is an SLD_body object (when extracted from its XML href) structured when using a cElementTree to parse it?

texts: I'm trying to export data from a PostGIS database into a geoserver and I've gotten it to where I can get the PostGIS database to point to the geoserver (I'm doing this as a python script for a GIS program) with gsconfig. So essentially I need to establish the services so that I can see the data in the geoserver data store/namespace that I have created. Therefore I'm using the xml.etree from cElementTree package in Python to parse the data from my database using gsconfig's catalog object. However, I don't know that much about XML and how its structured, but I have the gist of the element tree structure.


texts: Create common vertex where line was snapped to another using "snap to segment"

texts:
I have a schematic diagram within QGIS consisting of a single layer of line segments.  Lines have been snapped to others using "snap to segment".  When using the snap to segment option the two lines that intersect do not share a common vertex.

I am attempting to use Road Graph to navigate the schematic, find paths and distances.  The issue is that the shortest path algorithm requires two lines that share a common point.

Is there a plugin, script, or otherwise known procedure to add a common vertex/point on two lines that have been snapped using "snap to segment"?


texts: Updating the TEXT attribute on an OpenLayers point

texts:
I have the following map and vector layer declaration that is working:
var map;
var vectorLayer = new ol.layer.Vector({
  target: "points",
  source: new ol.source.Vector(),
  style: new ol.style.Style({
    image: new ol.style.Icon({
      anchor: [0.5, 0.5],
      anchorXUnits: "fraction",
      anchorYUnits: "fraction",
      src: "RedDot.svg",
      }),
            label: "My Junk",
            text: new ol.style.Text({
                font: '15px Calibri,sans-serif',
                fill: new ol.style.Fill({ color: '#000' }),
                stroke: new ol.style.Stroke({
                    color: '#fff', width: 2
                    }),
                text: "My Truck",
                offsetY: 18,
                })
    })
  });

(key info - Note the default "My Truck" label for each point)
I have the following function that adds a point to the map (I don't honestly know if it's creating a new layer for each point, or adding all the points to the default layer):
function add_map_point(lat, lng, txt) {
  var source = vectorLayer.getSource();
  var pointFeature = new ol.Feature({
    geometry: new ol.geom.Point(ol.proj.transform([parseFloat(lng), parseFloat(lat)], 'EPSG:4326', 'EPSG:3857')),
    });
  source.addFeature(pointFeature);
  }

This all works - no issues.
What I would like to do is use the TXT variable passed to that function to update the TEXT label ("My Truck" by default).  It seems like I should be able to do something simple like source.getStyle().setStyle(getText().setText(text: txt)); (over simplified I assume).
Could someone point out how this could be done?
I've been searching all morning, and all I can find is redeclaring the style again, which seems unnecessary.


texts: Setting fill colors using hex value in ArcGIS Online (arcgis.com) webmap?

texts:
Has anybody been able to use hex values (e.g. #FF0000) for setting polygon fill symbols in an ArcGIS Online (arcgis.com) 'webmap'?

I see predefined colors to select from, but no hex input box.


texts: How to obtain overlapping semi-transparent circles effect in ArcGIS for Desktop?

texts:
I have been trying to obtain the following markers effect.

How can I do the same in ArcGIS?

This image is from a book.


texts: Overlaying and exporting rasters in QGIS?

texts:
I have in QGIS a higher quality GeoTIFF overlayed on a larger, lower quality GeoTIFF. I reprojected one raster into the coordinate system of the other.

How can I export these to a single GeoTIFF while retaining the resolution of the high quality raster?

I've tried using the Merge tool (with no data set to 255 to remove whitespace) but get an error at 50%: "Illegal band #." The resulting layer is low quality and much brighter.

I also tried using Build Virtual Raster, but the resulting layer turns red and green (with "Place each input file into separate band" checked). Without it checked, the high quality raster doesn't show up in the new layer.

The high quality raster has 3 bands and the lower quality raster has 4 bands. There is no information by the 4th band though


texts: How to create accurate heat maps and hot spot maps in QGIS?

texts:
I am mapping 5 heritage sites using DGPS point data. Each site has 100s, and in two  cases 1000s, of points representing various attributes. I want to show the density of attributes and then demonstrate a threshold (no. of attributes) in hot spots for each heritage site.

I have no problems creating heat and hot spot maps following the QGIS instruction manual and Ujaval Gandhi’s tutorials, but I am uncertain about how to use the advanced raster heatmap plugin settings to set the x and y cell size, kernel shape etc.  I seem to get different results every time I create a map.  So, is it ok to leave the advanced settings off? and just choose the radius? I set the radius to 10m, 5m and 1m depending on the attribute I want to display.

I need the heat mapping to be representative and the hot spot mapping to be accurate.


texts: I can't find the color sampler (to select a color from anywhere on the desktop and have this copied into QGIS) in QGIS 2.6. I attach a screenshot on where QGIS says it is meant to be.


texts: type mismatch: 'esri_8' when using field calculator

texts:
I am trying to run the VBScript 'round' function using the field calculator on a column for a numeric field that contains a large number of rows (> 500k).  When I run the code on a subset of the field (~100 rows), it runs perfectly fine, but when I try to run it on the entire table I get the subject error. It says 'there was an failure during processing, check the geoprocessing results window for details'.  Does anyone know what could be causing this?

This is an annotation feature class that I would like to keep in the ArcMap environment while editing to avoid having to import it back.

my code is in the field calculator is: Round([TextString], 1)

The annotation feature class is currently in a file geodatabase for testing purposes and then I will run it in SDE, ArcMap 10.2.2, Sql-Server, Windows 7 64-bit SP1.


texts: PyQGIS: Load, process & overwrite layers from GPKG

texts:
I have 3 layers in one geopackage, all with the same name but different geometry: Point, Line, Polygon (originally coming from a MapInfo .TAB, which supports multi-geometry in one file). I already have loaded them in my project. Now i want to run native:centroids over the polygon-layer to get the centroids and save the result back into the Geopackage (and overwrite the existing point layer + delete if possible the line and polygon layer).
My script so far:

from qgis.utils import iface
from qgis.core import *
import os

layers = qgis.core.QgsProject.instance().layerTreeRoot().layerOrder()

for lyr in layers:
    if lyr.type() == 0 and lyr.name() == 'lightPosts':
        basename = lyr.name()
        geom = lyr.wkbType()
        layer_path = lyr.source()
        layer_path = layer_path.replace("Polygon","Point")

        input_layer = 'lightPosts'
        if geom == 3:
            params = {'ALL_PARTS' : False, 'INPUT' : input_layer, 'OUTPUT' : layer_path}
            processing.runAndLoadResults("native:centroids",params)

Following error message it prints:

_core.QgsProcessingException: Could not create layer H:/Test_Data.gpkg|layername=lightPosts|geometrytype=Point.gpkg:
Creation of data source failed (OGR error:
sqlite3_open(H:/Test_Data.gpkg|layername=lightPosts|geometrytype=Point.gpkg)
failed: unable to open database file)

It seems it can not acces or overwrite the layer already contained in the geopackage. How can i circumvent this problem?


texts: How to recover deleted data from an edited point layer?

texts:
As a newbie in GIS, I didn't fail to make a newbie mistake. While editing a point layer I accidentaly deleted all the introduced data. I duplicated a layer, as it was more convenient (same style and attributes), but forgot to save it as a new layer before editing it and saving it. It resulted in the deletion of all the previously introduced points.

It appears nonetheless that the datas haven't totaly disapear. If I open for instance the related .dbf file with notepad, it still show the deleted lines (not the case in excell though). Same is true for the shapefile that weighs much more than the empty one I saved later.

Is there a way to retrieve the lost data and to make them reappear on the layer?

Thanks for your help.


texts: Creating a polygon with four coordiates?

texts:
I need to draw polygon shape with the following coordiantes

private List<IGeometry> geos2 = new List<IGeometry>();
  var gf = new GeometryFactory(new PrecisionModel(), 3857);
                cor1 = new Coordinate(16144349.4032217,-4291402.04717672);
                corner2 = new Coordinate(16145538.9946882,-4291401.77806892);
                corner3 = new Coordinate(16145967.1248121,-4292444.80890257);
                corner4 = new Coordinate(16145352.7123475,-4292584.54329502);
                corner5 = new Coordinate(1500000, 1500000);

I tried with following line of code to draw

 geos2.Add(gf.CreatePolygon(new Coordinate[] { cor1, corner2, corner3, corner4 }));

But noting happening . but i could create Point with the following code

geos2.Add(gf.CreatePoint( corner2));

Could you please tell me how to make a polygon.


texts:
I classified this image in ArcGIS 10.3, but my Professor said the pixels are big and I should make it smaller.
This is my problem: is there any tool that can use to make the picture grains finer. Or are there any parameters i can adjust in order to make the grain size smaller.


texts:
Doubled question (https://stackoverflow.com/questions/44726971/changing-a-raster-stack-projection-in-r?noredirect=1#comment76436491_44726971)

I have a NetCDF (nc) file I'm using to create a Raster stack in R.  The coord. ref. is listed as "NA" when I load it in, but plots as Polar Stereographic.

When I try set its CRS, it lists that it has changed in the Raster info, but still plots incorrectly.  How do I change the CRS to WGS1984 and plot it as Equirectangular?

I need to crop it using lat/long coordinates etc, so can't really leave it as stereographic projection.

Data is available here https://drive.google.com/file/d/0Bz0W7Ut_SNfjeE9ObXpySzJ5UWs/view

Code used is:

  library(ncdf4)
  library (raster)
  library(sp)
  nc<-nc_open('NOAA_Snow_JanJune2016.nc')
  print(nc)
  varname<-names(nc$var)
  r<-stack('NOAA_Snow_JanJune2016.nc',varname='snow_cover_extent')

  crs(r) <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"
  plot(r[[1]])

(I can provide images of plots etc if you don't want to run the above)


texts: How can I design a similar Low Pass Filter as given below which takes into account only 4-connectivity and ignores the pixels that complete 8-connectivity?


texts: Turn certain layers on and off in ArcGIS Web AppBuilder Developer Edition

texts:
Does anyone know a way in ArcGIS Web AppBuilder Developer edition that I can create a button or a link that turns certain layers on or off?

For example, if I have three layers

[1] Homes
[2] Rivers
[3] Zip Codes

and I want to create a button that turns on layers [1] and [3], so that the user only has to click once to do this instead of twice.

(In reality, there are more like 20 layers, so the time savings would be significant)


texts: How to prevent cross-world lines in ggplot world map?

texts:
I'm having issues displaying 'great circles' paths and polygons in ggmap, an extension of ggplot2. Please enter the following code to replicate my issues:

https://gist.github.com/Robinlovelace/379e0a3db23b53965c73

# load packages we'll be using
x <- c("ggmap", "geosphere", "sp")
lapply(x, library, character.only = TRUE)

download.file("https://dl.dropboxusercontent.com/u/15008199/tmp/origins.csv", "origins.csv", method = "wget")

origins <- read.csv("origins.csv")
os <- SpatialPoints(coords = origins,
proj4string = CRS("+init=epsg:4326"))
  dest <- geocode("Leura")
destp <- SpatialPoints(coords = matrix(c(dest$lon, dest$lat), nrow = 1),
  proj4string = CRS("+init=epsg:4326"))

path <- NULL
for(i in 1:nrow(coordinates(os))){
  path_tmp <- gcIntermediate(coordinates(os[i, ])[2:3], p2 = dest)
  path_tmp <- data.frame(path_tmp)
  path_tmp$person <- i
  path <- rbind(path, path_tmp)
}

wrld <- borders("world", colour="gray50", fill="gray50")

# plot the result
ggplot() + wrld +
  geom_point(data = origins, aes(lon, lat)) +
  geom_path(data = path, aes(lon, lat, group = person)) +
  coord_fixed(ratio = 1)

As shown, all lines crossing 180 degree longitude jump across the whole world, instead of 'going round the back' as it should. Any ideas of a quick fix?

Even when the destination is moved to the centre of map there are still problems:

On a wider note, are there any prospects of the bug in which projected geom_polygons fail? - This issue is mentioned on the moribund ggmap GitHub issue tracker and described on in a blog post.

Beginning to wonder if it would be easier to do this in base graphics...

Edit: Example II

Download data

x <- c("ggmap", "geosphere", "sp")
lapply(x, library, character.only = TRUE)

download.file("https://dl.dropboxusercontent.com/u/15008199/tmp/origins.csv", "origins.csv", method = "wget")

origins <- read.csv("origins.csv")
os <- SpatialPoints(coords = origins,
proj4string = CRS("+init=epsg:4326"))
dest <- geocode("Leura")
destp <- SpatialPoints(coords = matrix(c(dest$lon, dest$lat), nrow = 1),
proj4string = CRS("+init=epsg:4326"))

path <- NULL
# for(i in 1:2){ # for testing
for(i in 1:nrow(coordinates(os))){
pathl <- gcIntermediate(coordinates(os[i, ])[2:3], p2 = dest, breakAtDateLine = T)
cl <- class(pathl)
path_tmp <- matrix(unlist(pathl, use.names = T ), ncol = 2)
path_tmp <- data.frame(path_tmp)
path_tmp$person <- i
# ifelse(cl == "matrix", path_tmp$person <- i , no = path_tmp$person <- paste(i, c(rep(1, nrow(pathl[[1]])), rep(2, nrow(pathl[[2]]))), sep = "."))
path <- rbind(path, path_tmp)
}

wrld <- borders("world", colour="gray50", fill="gray50")

# plot the result
ggplot() + wrld +
geom_point(data = origins, aes(lon, lat)) +
geom_path(data = path, aes(X1, X2, group = person)) +
coord_fixed(ratio = 1)

I've got it working by preventing lines from crossing the longitude extremes and following code described here but still not found a solution to the original question.


texts: Rasters changing pixel size and range of values when imported into QGIS?

texts:
I have observed that when I import rasters into QGIS their pixel size changes and also, the full range of the raster's value does not appear. For instance, I may have a raster with values ranging from 0 to 2890, and QGIS shows my values to range from 10 to 2850.

Does anyone know why this happens and how can it be fixed?


texts: Problems with Saga integration in QGIS 2.14.4 LTR

texts:
I've installed QGIS 2.14.4 LTR (32 bit) on my WIN7 64bit system.
But I believe something is wrong with the "Processing" plugin as many of the features were missing from what I've used earlier.
My QGIS window didn't have any Simple/Advance interface option for switching between various providers in the processing toolbox.

Also the option for setting path was not seen in processing>>Options>>Providers>>SAGA
The toolbox showed 235 SAGA Algorithms(with options such as "Grid Tools" still missing) but running them resulted in following error msg

 *****************************************************************
 Uncaught error while executing algorithm
        Traceback (most recent call last):
          File "C:/PROGRA~2/QGIS21~1.4/apps/qgis-ltr/./python/plugins\processing\core\GeoAlgorithm.py", line 203, in execute
            self.processAlgorithm(progress)
          File "C:/PROGRA~2/QGIS21~1.4/apps/qgis-ltr/./python/plugins\processing\algs\saga\SagaAlgorithm212.py", line 134, in processAlgorithm
            filename = dataobjects.exportVectorLayer(layer)
          File "C:/PROGRA~2/QGIS21~1.4/apps/qgis-ltr/./python/plugins\processing\tools\dataobjects.py", line 321, in exportVectorLayer
            if not os.path.splitext()[2] in supported or not isASCII:
        TypeError: splitext() takes exactly 1 argument (0 given)

The issues are presented/highlighted in the attached image.

I've tried the OSGEO Network installer as well as the standalone installer. I've uninstalled qgis several times and subsequently cleaned the registry and recycled the system (deleted the .qgis2 folder) for the next installation. But the result is the same.


texts: Drawing the boundaries on coordinate file and classifying the points

texts:
I have imported the data on coordinates in Excel which are classified into 3 types: Class 1, Class 2 and Class 3. When I Add delimited text layer, I get the following output:

I have two basic questions:

How can I draw boundaries around these points?
How can I classify these points into Class 1, Class 2, Class 3 so that each class is distinctly represented with a different color. Right now, everything looks the same. I just need each point to be distinctly represented and not boundary around each point!

I am using QGIS 2.4


texts:
Possible Duplicate:
  Why does QGIS scalebar show a wrong scale?

I want the scale-bar in QGIS to be in meters/kilometers. Usually it's in degrees, but if I change it to meters it is displayed in cm. What's wrong? (Version is 1.7.4, coordinate system is WGS84)
The units are wrong in the scale bar on the normal view and when I want to add a scale bar in the print composer, but the "measure"-tool gives correct numbers. It seems that one "cm" is 1km in reality.
I tried a range of different coordinate systems (projected and unprojected) but it's still the same. Part of the data comes from Open Street Maps part I collected with my GPS, so WGS84 should be correct, no?

Can anybody help? Thanks a lot!

This question has already been asked by an other user, but with no answer that solved my problem (Why does QGIS scalebar show a wrong scale?)


texts: How to update feature class by iterating through values with branching?

texts:
I'm really struggling with getting this to work.  Basically, I want my script to get the "newID" value and use that as a look up value.  If the "newID" is found in the OldID list then do one thing, else do something else.  The problem I'm having is it only grabs the first ID value and then does loops the rest of the code based on whatever category the first ID fell under.  I need it to get the value of each field in the new data table, and write it to the corresponding field in the old data table.  It does this just fine, but obviously it isn't checking the ID's after the first one, or updating/appending the geometry.

Here's what I've got so far.  Just showing applicable code....

rows = gp.SearchCursor(new_data)
row = rows.next()

while row:
    newID = row.getValue('ID')
    newShape = row.SHAPE
    y = 0

    while newID in OldID_List:
        gp.AddMessage('')
        print 'Starting the updating for ratios'
        where = 'ID' + ' = ' + str(newID)
        new_rows = gp.UpdateCursor(old_data, where)
        new_row = new_rows.next()

        while y <= len(Field_List):  # checks to see if ID exists -
                                     # this means just update
            for each in Field_List:
                newValue = row.getValue(each)

                if newValue is None:
                    new_row.setValue(str(each), 0)
                else:
                    new_row.setValue(str(each), str(newValue))

                y += 1

                if y > len(Field_List):
                    new_row.Shape = newShape
                    new_rows.UpdateRow(new_row)
                    new_row = new_rows.next()
                    row = rows.next()
                    newID = row.getValue('ID')
                    newShape = row.SHAPE
                    y = 0

    while newID not in OldID_List:  # these are new objects to be appended
        flag = 1
        nsert_Rows = gp.InsertCursor(old_data)
        nsert_Row = nsert_Rows.NewRow()

        while y <= len(Field_List):
            for each in Field_List:
                newValue = row.getValue(each)
                print 'Starting appends'

                if newValue is None:
                    nsert_Row.setValue(str(each), 0)
                else:
                    nsert_Row.setValue(str(each), str(newValue))

                y += 1
                if y > len(Field_List):
                    nsert_Row.ID = newID
                    nsert_Row.Shape = newShape
                    nsert_Rows.insertrow(nsert_Row)

                    row = rows.next()
                    newID = row.getValue('ID')
                    newShape = row.SHAPE

                    appends += 1
                    y = 0
del rows
del new_rows
if flag == 1:
    del nsert_Rows

Any help or suggestions would be appreciated.

Thanks.


texts: Create foldable section in QGIS tool using R script

texts: I am developping a GUI tool in QGIS using an R script and I am wondering how to add this kind of foldable section in the graphical interface.


texts: Circles created by points inside bounding box query - PostGIS

texts:
The table that i want to perform this query has 3 columns: id, the_geom and seen_count.

The query i've been using to get points inside a bounding box is as follows:

SELECT id, ST_AsGeoJSON(the_geom) FROM table WHERE ST_Contains(ST_MakeEnvelope({0}, {1}, {2}, {3}, 4326), table.the_geom)

What i want to do is: selecting points, which are transformed to circles using their seen_count as radius, intersect with bounding box. Is it possible? If so, would it be an expensive query? Thanks in advance.


texts: Creating quick coarse scale preview (JPEG/PNG) of raster mosaic using GDAL?

texts:
The script I developed stacks and mosaic satellite imagery, and the resultant imagery is in some cases "50GBs" in size. Therefore opening the image in a GIS environment to check the mosaic is sluggish and slow.

Is there anyway I can quickly create a JPEG/PNG overview/preview (it can be very coarse) of the big mosaic so that the users can quickly check the imagery using the preview file instead of opening/loading the 50GB image mosaic?


texts: Obtaining Postcode Boundaries for India?

texts:
I was wondering if there are shapefiles for postcode regions in India.

I did find some shapefiles that give the outline of the administrative units in the country but not with postcodes. Right now I am using Google Map API that highlights the region under pin-code.

It looks like the information is not present in public domain. It looks like either you can crowdsource the information or buy it from a thirdparty called MapMyIndia.

I found a blog by Anand.S where he wrote about his crowdsourcing idea. Other than this I couldn't really find anything interesting: mapping-pin-codes


texts: Removing a group from layers panel using PyQGIS?

texts:
How can I remove a group that I created previously using :

groupName="machin"

root = QgsProject.instance().layerTreeRoot()

group = root.addGroup(groupName)


texts: How to increase resolution of 3D image in GRASS NVIZ

texts: I am trying to make 3D visualization of terrain data and orthophoto using GRASS NVIZ  via QGIS. I have been able to import both the DTM and the orthophoto in NVIZ window and successfully draped the image over DTM. But the resolution looks very poor. Is there any way to increase the resolution ? It seems like NVIZ GUI is not so user-friendly. For example, in Arcscene, there is a option to optimize the rendering that  can increase the resolution of the 3D draped image. Here is the screeshot from GRASS NVIZ:


texts: Polygon scratch file is converting to Multipolygon shapefile in QGIS

texts:
I digitized somewhat complex (many vertices) boundaries of neighborhoods in QGIS as a Polygon shapefile.  Turns out QGIS thinks they are Multipolygons and Tableau doesn't like them.  They are not multipolygons as per the definition of a multipolygon.

I copied the features into a Polygon scratch layer and confirmed that the wkbType is Polygon.   I then saved the scratch layer using both the Export - Save Features as and the Make Permanent options.  In both cases, it converted the Polygon scratch data to a Multipolygon shapefile.

I have many uses of the data, mostly in QGIS and Python code which don't care, but Tableau is treating the polygons all as one polygon.


texts: Creating query for set of polygons which would select adjacent polygons for which combined area is greater than number?

texts:
I have a set of parcels, represented as polygons shapefiles and I need to build a query, to select all the polygons, which are adjacent, with the sum of the areas of 2 polygons,greater than 1000 acres.

Is this even possible to query in arcmap/pro, or it requires scripting?

Google was so not so keen in providing an answer to such a complex question. I know how to do simple queries in Arcmap, and calculate areas and joining adjacent polygons, but this is a little bit more than what I know.

I could do it manually, but this is going to be a recursive task, so need some automation.


texts: openlayers - how to render vector layer in openlayers with WebGl

texts:
I was used this tutorial and copy code of tutorial step by step but don't work for me.

As mentioned in the tutorial, I was used ?renderer=webgl as a query string.

here is jsfiddle of my code: https://jsfiddle.net/4gvdsjwj/2/

What is wrong?


texts: importing files from grass mapset to QGIS Layer Panel

texts:
I am performing an NDVI calculation against 2 sets of rasters. One produced by LANDSAT5 the other by LANDSAT8.

r.mapcalc "L8_NDVI = float(L8Band5 - L8Band4) / (L8Band5 + L8Band4)"

The L8_NDVI file gets created successfully in my GRASS mapset.  I see it in all of the subdirectories (cell, cellhd, hist,etc.), but It does not appear in the QGIS Layers Panel. Is there a GRASS "Show Output On Map" command that I need to execute as follwoup to my ndvi command?

I tried the simple r.in.gdal command to import the layer form the mapset to the QGIS map, but it prompts me to overwrite the output file which also happens to be the input file.


texts:
In the past I have just manually matched the color of the border to the color of the fill, but this can be time consuming.

I am wondering, is there any way to have the program automatically match the border color to that of the fill?

The fill color is being assigned through the styling tab in layer properties.

Using QGIS 2.8.6


texts: Oracle SDO_UTIL.GETVERTICES returns more than 1 vertex

texts:
I'm using Oracle 12.1.0.2.0 and I am getting extrange results using sdo_util.getvertices.

This query returns 2 "vertex" instead of one (as expected)

 SELECT T.x, T.y
    FROM TABLE(
      SDO_UTIL.GETVERTICES(
        MDSYS.SDO_GEOMETRY(
          2001,23030,NULL,
          MDSYS.SDO_ELEM_INFO_ARRAY(1,1,1,3,1,0),
          MDSYS.SDO_ORDINATE_ARRAY(662301.48,4203898.34,1,0)
        )
      )
   ) T;

         X          Y
---------- ----------
 662301,48 4203898,34
         1          0

The expected result would be only 1 vertex (the first one)

I'm afraid that I did not configure something well.

Does anyone know what could happen?

Maybe some incorrect database/session parameter.


texts: How to correctly set the Independent Cities of Virginia in Open Street Map

texts:
I would like to use the data from Open Street Map to identify which county I am in, but I noticed something unusual with the Independent Cities of Virginia. Essentially, there are a few cities that are independent of counties, even if they are completely contained within a county. When I do a query (Via the Map Quest Geocoding API), it returns that I am in the surrounding county, and not the city. For instance, here is Augusta County, which contains two independent cities, Staunton and Waynesboro. I've seen this for a half dozen other independent cities, BTW.

It's a bit tricky to see, but Staunton has a county boundary approximating the city boundaries. When I call Staunton, it correctly shows Staunton City as the county. When I look at Waynesboro, it stays it is in Augusta, not Waynesboro City. What I'm trying to do is to figure out how to make the appropriate correction. I can see a couple of possibilities, including:

Somehow trace an interior boundary that matches the city of Waynesboro, and identify that county as Waynesboro City.
Somehow make the existing Waynesboro think it's a county.

I'm not quite sure how either of those could be done, or which is the best, or even if there's a better solution all together. How can I make this correction?


texts: "All the inputs are not current" after Precondition = false in Model Builder

texts:
I am using an Iterator to Iterate through a list of Strings for Feature Selection/Aggregation/Count. There is the possibility that no Features are aggregated. In this case, I want the iterator to go to the next %String%. However, in such a situation, the model runs as follows:

Executing (Iterate Field Values): IterateFieldValues ...
Succeeded ...
Executing (Select): Select ... Data.shp in_memory\TagFeaturesSelect ... LIKE '%;seo;%' - ...
Executing (Get Count): GetCount in_memory\TagFeaturesSelect
Row Count = 11
Executing (Calculate Value): CalculateValue check(11) "def check(rowcount)... Boolean
Value = true
...
Executing (Aggregate Points): AggregatePoints in_memory\TagFeaturesSelect ...
Finding point clusters...
Constructing polygons...
...
Executing (Get Count 2): GetCount C:\01_TOOLS\Intermediate\SelTaglocationsAgg.shp
Row Count = 0
The process did not execute because the precondition is false.
Succeeded at ...
Executing (Spatial Join): SpatialJoin ...\SelTaglocationsAgg.shp in_memory\TagFeaturesSelect ...\output_seo.shp ...
Start Time: ...
All the inputs are not current.
Succeeded at ...
Executing (Spatial Join (3)): ...
Start Time: ...
All the inputs are not current.
Succeeded at ...
Executing (Summary Statistics): ...
Start Time: ...
All the inputs are not current.
Succeeded at ...
...

This is the Model I am talking about:

I also tried using Calculate Value with result = false for Row Count = 0, but the same result. The Model also ignores whether I set [Row Count 2] as a precondition for [Spatial Join (3)] or not.
My model still continues, but I get all the errors and I believe, this is not by intention?


texts: Problems using 2 maps with multiple layers and loadend event

texts:
I have 2 or more maps with one ore n layers on each map.
Now I want zoom to extent for each map, but the getDataExtent returns only for the last layer in the last map a bounds.

Can you help me?

Here is the code

var jekyllMapping = (function () {
    return {
        mappingInitialize: function() {
            var maps = document.getElementsByClassName("jekyll-mapping");
            for ( var i = 0; i < maps.length; i++ ) {
                var zoom      = maps[i].getAttribute("data-zoom"),
                    lat       = maps[i].getAttribute("data-latitude"),
                    lon       = maps[i].getAttribute("data-longitude"),
                    layers    = maps[i].getAttribute("data-layers"),
                    title     = maps[i].getAttribute("data-title"),
                    map, markers, center;

                maps[i].setAttribute('id', 'jekyll-mapping' + i);

                markers = new OpenLayers.Layer.Markers("Markers")
                map = new OpenLayers.Map('jekyll-mapping' + i);
                map.addLayer(new OpenLayers.Layer.OSM());
                map.addLayer(markers);
                if (lat && lon) {
                    center = new OpenLayers.LonLat(lon, lat).transform(
                        new OpenLayers.Projection("EPSG:4326"),
                        new OpenLayers.Projection("EPSG:900913"));
                    map.setCenter(center, zoom);
                    markers.addMarker(new OpenLayers.Marker(center));
                }

                if (layers) {
                    // layers are url's to kml files
                    layers = layers.split(' ');
                    while (layers.length > 0){
                        var m = new OpenLayers.Layer.Vector("KML", {
                                strategies: [new OpenLayers.Strategy.Fixed()],
                                protocol: new OpenLayers.Protocol.HTTP({
                                            url: layers.pop(),
                                            format: new OpenLayers.Format.KML({extractStyles: true, extractAttributes: true, maxDepth: 2 })
                                })
                        });
                        if (m) {
                            m.events.register("loadend", m, function () {
                                // raises for all layers but only the last layer on second map gives me a bound, why???
                                var bounds = m.getDataExtent();
                                if (bounds) {
                                    map.zoomToExtent(bounds);
                                }
                            });
                            map.addLayer(m);
                        }
                    }
                }
            }
        }
    };
}());

window.onload = function() { jekyllMapping.mappingInitialize(); }


texts: Shp2pgsql encoding issue / error. Special latin /utf characters (áéíóú) are wrong

texts:
I am not able to load accented vowels coming from a Shapefile, to a PostGis geodatabase.

This is my information:
Operative System: Windows 10.
PostgreSql Database: postgresql-10.5-2-windows-x64.exe.
Database encoding:  UTF8.
Database collation: English United States 1252.
Database character type: English United States 1252.
Sql client in Windows: psql.
Windows CMD default code page: 437.
Tool to create SQL from Shapefile: shp2pgsl.

1) Commands that I have ran to create the Sql:

shp2pgsql -s 21892 -c -g geom -D -k -i -I -W "UTF-8" municipios2.shp public.municipios2 > mun1.sql

shp2pgsql -s 21892 -c -g geom -D -k -i -I -W "LATIN1" municipios2.shp public.municipios2 > mun2.sql

On top of mun1.sql and mun2.sql (first line) reads:

SET CLIENT_ENCODING TO UTF8;

2) This is the encoding that Notepad++ shows for those sql files:

UCS-2 LE BOM

3) In Notepadd ++, next sample words, are inside those sql files:

CI├ëNAGA  (must be CIÉNAGA)
UNI├ôN (must be UNIÓN)

4) I tried in Notepad ++ to open those sql files with a different encoding or convert to a different encoding (UTF-8) but it does nothing or shows worst results:

CI├â┬ëNAGA (must be CIÉNAGA).
UNI├â┬ôN (must be UNIÓN)

5) This is the procedure that I did to load some of those files inside Postgis
- In windows command prompt (before use shp2psql of psql) I changed the code page:

chcp 1252

Inside psql, before load, I changed the client encoding to one of those options:

SET client_encoding = 'UTF8';
SET client_encoding = 'LATIN1';

I load (with psql) those files to the geodatabase:

\i 'D:/mun1.sql'

with this error
ERROR:  invalid byte sequence for encoding "UTF8": 0xff
psql end.

Note: If I convert those sql files to UTF8, they will load to PostGis, but with wrong words (CI├â┬ëNAGA, UNI├â┬ôN)

Of course if I have tested "Postgis 2.0 Shapefile Loader", and it works very well with UTF8. But I need shp2psql to use bath load of many shapefiles with specific sql sentence to load.

What should I do to solve the problem?

I have ogr2ogr option (I have not work with that), but I have been using shp2psql and I want to know what is the problem.


texts:
My goal is to install the latest build version of QGIS in order to test and identify bugs.  I'm on OS X (10.10.4) and am attempting to build from the QGIS nightly source.

I'm following the official instructions here.  I went to the Github QGIS page and downloaded the .zip for the Master version, then installed all the requirements.  When I run make, I quickly get the following error:

[  0%] Generating qgsversion.h
fatal: Not a git repository (or any of the parent directories): .git
make[2]: *** [qgsversion.h] Error 128

The build then progresses to 7% complete, then ends with:

[  7%] Built target ui
make: *** [all] Error 2

Any ideas on how to troubleshoot the above error?

I also am open to an alternative way to install and be able to test the QGIS nightly builds.  (Note: The qgis.dakotacarto.com nightly build appears to not be updated.  As of October 15, 2015, qgis.dakotacarto.com QGIS Master 2.11 appears to be updated.  See update 4 below for details.)

Update 1:
I cloned the repository by running the following in Terminal app:

clone https://github.com/qgis/QGIS.git

The CMakeError.logreturns the following errors:

Determining if the Q_WS_X11 exist failed with the following output:
Change Dir:     /Users/mbp/Desktop/qgis_git_clone/QGIS/build/CMakeFiles/CMakeTmp
...

Determining if the Q_WS_WIN exist failed with the following output:
Change Dir:     /Users/mbp/Desktop/qgis_git_clone/QGIS/build/CMakeFiles/CMakeTmp
...

Determining if the Q_WS_QWS exist failed with the following output:
Change Dir:     /Users/mbp/Desktop/qgis_git_clone/QGIS/build/CMakeFiles/CMakeTmp

Update 2:
Got same result when getting QGIS master via the following Terminal commands:

git init
git remote add -f -t master -m master qgisupstream git://github.com/qgis/QGIS.git
git merge qgisupstream

Update 3:
I updated to 10.11 El Capitan for unrelated reasons.  I thought it worth mentioning though that I got the same error as noted in update 1.

Update 4: qgis.dakotacarto.com site has updated version QGIS Master 2.11.  I did need to update my GDAL via the Kyngchaos GDAL Complete 1.11 framework package.


texts:
I'm struggling whit AGOL and related tables.
I've made a gdb in Arcmap and related a child table "test" to its parent feature class "ca" with the tool "Create relationship Class" and a double field as primary and foreign key.
Everything works perfectly into ArcMap:

When I publish it as service in AGOL, the related record exist but it's empty and says me "No results found". The attached table is empty as well... Any idea? Thanks !


texts: How to apply changes in QSettings immediately?

texts:
With PyQGIS, you can change various settings such as making the groups in the Layers Panel (or table of contents) bold:

from PyQt4.QtCore import QSettings
QSettings().setValue("/Qgis/legendGroupsBold", True)

The settings are written which can be seen from the options or from typing:

QSettings().value("/Qgis/legendGroupsBold")
>>> u'true'

But it does not get applied. I thought using

QSettings().sync()

might help but unfortunately not. The only way for it to apply would be to either go to the GUI (Settings > Options) then click OK or restart QGIS.

Is there a method to doing this through PyQGIS?

Context:

I want to add this as a macro in a project file to ensure (for aesthetic purposes) that only group names are bold and layer names are not. Therefore, it would be ideal if the QSettings are applied when the project has loaded (the initial settings would be reapplied when the project closed).


texts: How to copy attribute values of fields in selected features of one layer to another Qgis programmatically?

texts:
I am trying to get the values of attributes of the selected features of a layer field by field and copy to another layer that I created in qgis with pyqgis.

Is here a example for this. Every example I found leaded me to the dead end.


texts: ArcGIS standalone application deployment

texts:
I'm trying to build a deployment for a ArcGIS standalone application. I'm using ArcGIS 10 Engine controls in a WPF application with desktop and spatial analysis license.
I'm using classes that only belong to desktop (IGxDialog, ...) and also engine controls (AxTOCControl, AxMapControl, ...)
I followed the example of ESRI for creating deployment for engine applications
The problem is installing runtime is not enough, but if i install the application on a pc that include an installation of ArcGIS desktop it works fine.

Now I need to build a complete installation that includes both desktop and engine DLLs without installing ArcGIS desktop, and I don't know how that can be accomplished
any suggestions???


texts: Select IDLE 2.7 (not 3.4) as default script tool editor for ArcPy with ArcGIS for Desktop?

texts: I'm trying to set my default editor in ArcMap (10.3) to IDLE 2.7, but despite using the 2.7 path, when I right-click on a script tool and select Edit, it opens IDLE 3.4:


texts:
I have just downloaded and installed Google Earth Pro, and I cannot get place names to show up at all.

Can you suggest a course of action?


texts: Calculating area of land cover within 30m of streams for whole U.S.?

texts:
I'm trying to do a very simple analysis where I need to find all land cover of a certain type within 30 meters of a stream. I have a raster from the National Land Cover Database (NLCD) from which I've pulled out the relevant pixels and the National Hydrography Dataset stream network.

All I need to do is buffer the streams by 30m and either:

Vectorize my NLCD land cover pixels and clip by the stream buffers, or ...
Rasterize my stream buffers and pull out the NLCD pixels that overlap them.

(I know these aren't quite equivalent methods in terms of results, but that's OK for the purpose of this question)

The issue is that I need to do this for the whole U.S. Even if I break up the streams and the land cover by county (i.e. into ~3,100 pieces) to avoid out-of-memory errors from buffering/clipping/etc. such large datasets, and it takes 1 minute to process each county (WAY optimistic, it seems to me), the analysis would still need to run for >50 hours!

Is there any way optimize this so I can get the runtime down to perhaps 24 hrs?

I am working in ArcGIS Desktop (ArcMap 10.3.1). My data is a line feature class in a file geodatabase (for the streams) and an 8-bit unsigned integer TIFF that I've imported into that same geodatabase. Everything is projected to USA Contiguous Albers Equal Area Conic.


texts:
I'm getting an error when trying to switch basemaps in my application(right now only works in IE9).

TypeError: Unable to get value of the property '_basemapGalleryLayerType': object is null or undefined

I'm using standard code provided in the ESRI samples. The full code is a mashup of various samples, so my error may be in the order I've placed things. However, I've tried reordering the code and haven't had any luck.

When I remove mapReady from the dojo.connect(map, "onLoad", mapReady); listener, the basemap gallery works, however it causes my initial basemap to fail and the mapReady function to fail as well (obviously). There is no link that I can determine between the initial basemaps and my mapReady function so I'm not sure what the deal is there...

Also, if I clear the browser cache the entire app will work, 1 time only, then on refresh, it fails as described above. Very inconsistent behavior.

Any help would be greatly appreciated.


texts:
I have shapefile named `band-tailed_pigeon.shp which is fine and no problem on displaying or geoprocessing on Desktop but in ArcPy when I try to merge some layes to it I am getting this error

ExecuteError: ERROR 000354: The name contains invalid characters
Failed to execute (Merge).

I tried to do the proccess with a simple name like sample.shp which it works but I HAVE TO keep the same name as we are using them in several projects.
Here is the code which I have
FinaloutLocation = "E:\\GIS\\Data\\Final\\band-tailed_pigeon.shp"
inFeaturesToMerge = [dissolve_1, dissolve_2, dissolve_3, dissolve_4, dissolve_5]
arcpy.Merge_management(inFeaturesToMerge, FinaloutLocation)


texts: Joining Fields from Multiple Tables to One Table in ModelBuilder?

texts:
I have 40 tables with about 2000 row with a shared ID that I want to merge into 1 table with 2000 rows. I only want to retain a distance field calculated within each table. The distance field is uniquely named Dist_miles_[name] with "[name]" as a unique value for each table. This way in the final table I will have: Key value, Dist_miles_[nameofAdvisor1], Dist_miles_[nameofAdvisor2], Dist_miles_[nameofAdvisor3]... Dist_miles_[nameofAdvisor40].

I would like to use iterate table and join field, and delete field to remove all the extraneous fields. However, the specificity of the variable I want to keep poses problems because I can't select it in the join field tool. I was hoping the join field tool would take a wildcard and be able to pull only the Dist_miles* into the final output table, but it doesn't.

Is there a workaround this?

Used Polygeo's tip to create modelbuilder. Works fine for single table. Tried to iterate on tables but doing that created an error at Calculate tool: "A field name was not found or there were unbalanced quotation marks. ERROR 999999: Error executing function.
Failed to execute (Calculate Field)."


texts: GeoServer problems on OpenGeoSuite 4.0 upgrade

texts:
We've tried to upgrade OpenGeoSuite 3.1 on an Ubuntu server to version 4.0 following http://suite.opengeo.org/opengeo-docs/installation/ubuntu/upgrade.html#installation-ubuntu-upgrade
However we are experiencing several problems afterwards:

WPS extention (which doesn't seem to be included in OGS 4.0 anymore) doesn't work after dropping the according jars in the lib directory (GetCapabilities & DescribeProcess do work, Execute not)
GML 3.x is not served

Has anyone experienced similar problems or can give us a hint? Can it be an issue in the Snaphot version?

A general question: Why does OpenGeoSuite use Snapshot versions of Geoserver and not the latest stable release?


texts: Accessing Mosaic Dataset Properties to Change Default Mosaic Operator using ArcPy?

texts:
I would like to change the default mosaic operator on mosaic datasets my script creates. However, I can only manually change the mosaic operator from mosaic dataset properties from ArcCatalog by right clicking the mosaic dataset. Is there a way this could be access and changed from the script itself?

This is the properties window that I am trying to access using arcpy:

I'm trying to accomplish something similar to this:

for mosaic_dataset in mosaic_datasets:
    arcpy.CreateMosaicDataset_management(gdb_path, mosaic_dataset, coord, "1","8_BIT_UNSIGNED")

which works fine, but the next step is to change the mosaic operator from python rather than manually changing it and I cannot find any documentation.


texts: Create GeoTIFF from JPG and JFW files with gdal

texts:
I have a JPG image and the corresponding geospatial information on a JFW wolrd file in a ETRS89 geodesical system and UTM coordinates (EPSG:25830)

5.0000000000
0.0000000000
0.0000000000
-5.0000000000
484752.5000000000
4811767.5000000000

I want to create a GeoTIFF on a SRS EPSG:4326 and extract the edge lat/lon coordinates [ulx uly lrx lry]. I am using the folowing command without success which maybe requires to do a warp before (although I don't think so because line 2 and 3 are null on JFW file), and I do not know how to extract the edge values on that projection.

gdal_translate -of GTiff -a_srs EPSG:4326 -projwin 0 0 [width_jpg] [height_jpg] [from_file_jpg] [to_file_tif]

Does someone know how to do it? It is necessary to define width and height of image if I am using the whole image?

EDIT

The code used to display the tile is the following one:

//WMS overlay
var format = 'image/png8';
var tiled = new OpenLayers.Layer.WMS(
           "Show Layer", "http://<IPaddress>:8080/geoserver/" + workSpace + "/wms",
           {
               LAYERS: layerName,
               STYLES: '',
               TRANSPARENT: true,
               VERSION: '1.3.0',
               format: format,
               tiled: true
           },
           {
               sphericalMercator: true,
               buffer: 0,
               displayOutsideMaxExtent: true,
               isBaseLayer: false,
               opacity: 1.0,
               yx: { 'EPSG:4326': true }
           }
);
map.addLayer(tiled);


texts:
I'm working with Leaflet and I'm also working with JSTS to do the spatial operations. In most of the cases I want to create Polygons inside another Polygon Geometry.
In this case the yellow shape is the selection. The light blue line is the polygon that will intersect with the selection (yellow shape) and create the new Polygon Geometry:

The intersect of both  will create the new blue polygon :

So far so good. But if we repeat the process and we need to create another Polygon (Note: overlaps are forbidden between blue polygons, if there is any overlap, it will generate the difference):

Up to this point it seems all is ok, both geometry are created correctly.

The problem is when I select the bottom blue shape. I can see that geometry didn't generate correctly.

I have this code to transform between WKT to geojson. Probably the issue it's here in the precision of my geojson. (I can't paste the code right)

    var precisionModel = new jsts.geom.PrecisionModel(1000000000000000000000000)
    var geometryPrecisionReducer = new jsts.precision.GeometryPrecisionReducer(precisionModel)
    var geomFactoryJSTS = new jsts.geom.GeometryFactory(precisionModel)
    var geoJSONReader = new jsts.io.GeoJSONReader()
    var geoJSONWritter = new jsts.io.GeoJSONWriter()
    var readerWKT = new jsts.io.WKTReader(geomFactoryJSTS)
    var isValid = new jsts.operation.valid.IsValidOp()
    readerWKT.read(WKT.GeoJSONtoWKT(geojson[i]))

Other posibility is that my yellow shape is in EPSG25831 and Leaflet transform it 'on the fly'. It is posible that Leaflet forget some precision here?

Anyone has idea how can I fix that?


texts:
I simply want to make a georeferenced imageOverlay draggable in leaflet.
I found several unanswered posts (e.g. here and here).

Since it is easy to make markers draggable, someone tried to adapt the code for imageOverlay here. The issue was closed and refered to the interactive option for events on imageOverlays.

According to the last two links, it shouldn't be so difficult but I have no idea how to do it - do you?


texts: Create a continuous raster buffer for polyline

texts:
I would like to create a continuous raster that is a buffer around a polyline in ArcGIS.  For example, lets use a maximum buffer distance of 20m.

At a 20m distance from the line, the raster value would be 0;
At a 0m distance from the line (i.e. right on top) the raster value would be 1;
At a distance of 10m from the line, the raster value would be 0.5;
At a distance of 5.25m from the line, the raster value would be 0.2625;

I realize I could create multiple buffer polygons around the line, then convert to raster and normalize, but is there any easier way to do this that also allows for a more continuous raster?


texts: Road network comparisons - how to extract roads that do not completely fall within the buffer

texts:
I am trying to compare how well one road network aligns with another - I have created a buffer around the first network and would like to select only the roads from the second network that fall COMPLETELY within that buffer, OR only have small extensions that go beyond.

I have done both a clip and erase of the second network by the buffer and joined those sets to find the set that fall both within and outside the buffer to remove them, but need to find a way to not remove those that "mostly match"

i.e. how can I remove the roads that do not match the buffer but are included because part of the road is within the buffer, as in a T intersection, without also removing those roads that have smaller parts not falling within the buffer?

It doesn't have to be a 100% solution, there is room for a reasonable amount of error


texts: Symbolizing Layers from a PostGIS Database to many Client Computers

texts:
I was hoping someone could help with a symbolizing issue.  I have a PostGIS database on a server with a number of layers.  The layers of a series of polygons, lines and points that all share the same symbology.  I want to symbolize the layers with the same symbol schema on computers in lab that are connected to the database.  I have used the MultiQML tool to symbolize the layers and saved them on the server, which did not load the correct symboled features on the connected lab computers.  Do I have to have the QML files on each lab computer, or is the a way to have a QML file on the server that each lab computer will be able to read via the database connection and symbolize the layers correctly without needing to store the QML locally?

Thanks,
Ryan


texts: How to choose a projection for area measurement on building to continental scales?

texts:
I have many geographical objects described as polylines, polygons, multipolygons. All objects are in WGS84. Object could be: a building, stadium, park, street, highway, city, state, even country and continent

I need:

calculate accurate area/length for each object.
calculate intersection are for objects.

I've read some topics here:
How to select projection

and there was a comment from @AndreJ:

  You are right for the points, but the whole continent is just too
  large for UTM. For the intersection, try EPSG:3857 or 4326. –  AndreJ
  Jan 29 '14 at 11:35

I don't understand it. 4326 is in degrees, I can't use it to measure area of smth, I need meter unit. 3857 is for "web", 3857 should give significant errors and distortions if used for computation...

I want to automate the process of course. I have billions of objects. Is there any java library for that? I used Geotool/JTS but it's only for planar geometry. I need to correctly project before using it :)

I see several approaches:

Dirty: Get object bounding box and select proper UTM zone for projection and use geo-tools java lib to do it.
If object is big enough (I have category for each object, I know if it's state/country/small lake in the forest)

2.1. Use some special projection for huge objects

2.2. Or split it (continent for example) by UTM zones, correctly project parts of object and them do calculation on top of these pieces.

I don't want to invent one more bicycle, probably we already have something ready for that?

EDIT:
please see that I need to calc area and area of intersection. I need to preserve shape to calc correct in intersection area.

Coordinate system for accurately calculating areas of polygons that cross UTM Zones? is about area, but not intersections...