texts:
I'm looking for a method to shift the poles vertically in Mollweide projection using QGIS 3.

I have come across this website and I want to do something similar.

This is the default view of the Mollweide projection by Jason Davies.

This is my take on Mollweide projection in QGIS 3 with minor alterations.

This is Mollweide projection by Jason Davies with the poles shifted vertically

How can I achieve the same result in QGIS 3?


texts: Selecting table records by attribute if in Python list

texts:
I would like to create a subset of a table based on values in a Python list that I created from a search cursor. I have not been able to get the where_clause correct and keep getting errors. Here is my code:

pu = 'PlanningUnits'
hab = 'lu_Habitat'

with arcpy.da.SearchCursor(pu, "unique_id") as cursor:
    selected_pu = sorted({row[0] for row in cursor})

where_clause = """{} IN selected_pu""".format(arcpy.AddFieldDelimiters(hab, "unique_id"))
selected_hab = arcpy.TableToTable_conversion(hab, r'in_memory', "selected_hab", where_clause)

I have also tried the following as a where_clause:

where_clause = """{} IN """.format(arcpy.AddFieldDelimiters(hab, "unique_id")) + str(tuple(selected_pu))

In both cases, I am getting the following error:

ExecuteError: ERROR 999999: Error executing function. An invalid SQL statement was used. An invalid SQL statement was used.


texts: instal pgrouting on postgres 9.3 + postgis 2.1

texts:
Just installed postgres 9.3 + postgis 2.1 on debian 7.

Postgis works fine, but installation of PGrouting fails in psql.

(no errors while installing sudo apt-get install postgresql-9.3-pgrouting)

But this SQL:

CREATE EXTENSION pgrouting;

gives

ERROR: could not load library "/usr/lib/postgresql/9.3/lib/postgis-2.1.so": /lib/x86_64-linux-gnu/libc.so.6: version `GLIBC_2.14' not found (required by /usr/lib/libproj.so.0)
SQL state: 58P01
Context: compilation of PL/pgSQL function "pgr_startpoint" near line 1

Any suggestions?

I also tried the combination postgresql 9.2 with postgis 2.0, same problem.

hans


texts: Assuming I start ArcGIS 8 times on my PC. One for each processor on my PC. Now lets say I run similar but slightly different models in each ArcGIS. Lets say reclass, zonal stats, or math. Now assume all 8 have the same inputs and outputs but the process is different. Could I get cross-pollinated results with all 8 using the same scratch directories or are the temp files tagged with the processor ID or similar to keep them distinct?


texts: Cut feature from the same layer with the buffer in QGIS

texts: Is it possible in QGIS to cut with a buffer a feature from the same layer, like it is in the editing mode for ArcGIS (see picture below)?


texts: QGIS Grid rasters have no resolution setting

texts:
In earlier versions of QGIS you could create rasters like heat map and interpolation, and define the resolution of the raster.

In 3.10 I see these dialogs are new, and raster size fields disappeared. So no matter what settings I use, the resulting raster is about 280 px in the biggest dimension.

How do I change the resolution?


texts: Getting information about what road is at certain latitude and longitude

texts:
I am using ArcGIS Desktop and I have a bunch of points that are part of a road network. I have created a shapefile using the starting and end point for each segment.

Is there a way to add a field to this data that corresponds to the road that it lies on?


texts: How can I create a RasterLayer object layers from a list of files stored as data frame?

texts:
I am trying to load raster files stored in AWS s3 bucket and I've managed to read them as data frame. How can I change them to rasterlayer objects?

library(aws.3)
library(raster)
library(rgdal)

#set aws s3 credentials
Sys.setenv("AWS_ACCESS_KEY_ID" ="XXXXXXXXXXXXXXXX",
       "AWS_SECRET_ACCESS_KEY" ="XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX",
       "AWS_DEFAULT_REGION"="us-east-1")

#read objects
r_files<-get_object_df(bucket="my-bucket")

#check dataset
str(r_files)

 > 'data.frame':    8 obs. of  8 variables:
   $ Key              : chr  "niger_LST_day_2004.tif"
  "niger_LST_night_2004.tif" "niger_black_albedo_2004.tif"
   "niger_builtup_2000.tif" ...
   $ LastModified     : chr  "2018-06-05T21:07:37.000Z" "2018-06-
     05T21:07:37.000Z" "2018-06-05T21:07:37.000Z" "2018-06-05T21:07:37.000Z" ...
   $ ETag             : chr  "\"f0659cc4d1f606acb571a3c510cff3ce\""
 "\"cee34ed55e0293957455e1576d0cba3f\"" "\"bdb26dfaac36b394c29aee55dde1834e\"" "\"5b7df07a486678b76daf5e1f08aeea0a\"" ...
   $ Size             : chr  "4615511" "4504193" "5581624" "490477" ...

str(r_files)

I followed the suggestion posted on Creating raster layer from data frame in R? changed the data frame to matrix and I get this.

# change to matrix
r_files_2 <-as.matrix(r_files)

str(r_files_2)
> chr [1:8, 1:8] "niger_LST_day_2004.tif" "niger_LST_night_2004.tif" ...
- attr(*, "dimnames")=List of 2
  ..$ : chr [1:8] "1" "2" "3" "4" ...
  ..$ : chr [1:8] "Key" "LastModified" "ETag" "Size" ...

and to create raster objects I will go with:

for (i in r_files_2){assign(unlist(strsplit(i, "[.]"))[1],raster(i))}

But, the following error shows up

> Error in .local(.Object, ...) :
  `/home/rstudio/aws_s3_testing/niger_LST_day_2004.tif' does not exist in the file system,
 and is not recognised as a supported dataset name.

What would be the best way to save each file to my home directory("aw_s3_testing") as tiff file?


texts: Cannot install QGIS 2.18 from repository on Debian Stretch; package "gdal-abi-2-1-1" not installable

texts:
I've been using QGIS from the repository (http://qgis.org/debian/ stretch main). I tried to do an update today (it had been a while) but it could not be updated because qgis : Depends: gdal-abi-2-1-1 but it is not installable.

I uninstalled QGIS, switched the repository to http://qgis.org/debian-nightly-release/, and tried to reinstall, but it gave me the same error.

 qgis : Depends: gdal-abi-2-1-1 but it is not installable
        Depends: libqgis-analysis2.18.0 but it is not going to be installed
        Depends: libqgis-app2.18.0 but it is not going to be installed
        Depends: python-qgis (= 1:2.18.0+git20161104+b23e29f+14stretch) but it is not going to be installed
        Recommends: qgis-plugin-globe but it is not going to be installed
        Recommends: qgis-plugin-grass but it is not going to be installed
        Recommends: qgis-provider-grass but it is not going to be installed

A bit of searching shows that this package is available in sid, but not stretch.


texts: Creating Shape by Density Criteria using QGIS?

texts:
I have a map with points plotted.

Is there a way that I can make a circle of a certain radius (let's assume 10 miles) if the density of points in any given 10 mile radius is over a certain number?


texts: Value at a particular spot on Geoserver Heatmap

texts:
Geoserver is drawing from a PostGIS database to generate a heatmap using this SLD extension. Leaflet then overlays the heatmap over the base layer.

Question: How can you determine the value at a particular spot on the heatmap? This value to display will be color ramp value X * constant C. In other words, how can we determine the color ramp value at that particular spot on the heatmap?

There's probably other methods of finding this value without using the heatmap's color ramp mapping.


texts:
I looked at this question, but it only answers part of my problem.
I want to create Thiessen polygons based off of the green points within the blue polygon boundary ONLY.  I've set the Processing extent to the exchange boundary, and have set the raster analysis mask to the same exchange boundary.

When I run the Thiessen polygons tool, it will still use the largest extent of the polygon, not the actual borders of the polygon.  You can see that polygons are created for points outside of the blue border.  Is there a step I'm missing?


texts:
I'm having hard time doing the following analysis in QGIS:

I have watershed polygons with individual IDs.

Then I have four different shape files, each having different land use polygons.

I would like to know the areas of each land use inside a individual watershed polygon.


texts: Stopping leader lines overlapping each other?

texts:
I am using ArcGIS Desktop 10.3.

I often have a problem where my leader lines are overlapping other leader lines in the same layer. I do not want to convert to graphics. I understand how to stop labels overlapping other features, it is the leader lines that are the problem. I use Maplex a lot and have altered all the settings I can think of. Sometimes changing the leader tolerance makes a slight improvement but not all the time.

Example below.

Here I have a long offset because some areas are very busy. All the external zones are allowed. These are currently advanced callouts but I'm having the same problem with simple ones.

Is there a setting I'm missing somewhere?


texts: For loop to calculate percentages using update cursor

texts:
I'm trying to run a for loop to calculate acreage percentages for a shapefile. I've converted field1 and field2 to float. I'm having trouble executing the CalculateField tool using the updateCursor. I want values where the area is greater than 0.000. I would also, eventually, like to run a more complex SQL query than only calculates percentages for a specified date range, values greater than 0.000 and values that don't have a null area. Is it possible to run multiple queries using the UpdateCursor?

Here is the part of the code I'm using to execute the calculatefield tool.

import arcpy
from arcpy import env
env.workspace = "X:/UserFolders/AnneGilligan/WC_74_Testing"
fc= "Export_Output_3.shp"
rows= arcpy.da.UpdateCursor(fc, '"Area" > 0.000')
newfield="ac_std"
field1="!acres!"
field2="!Area!"
expression1="(!acres!/!Area!) * 100" #standard calculation
expression2="(!Area!)/!acres!) * 100" #standard calculation
expression_type="PYTHON_9.3"
code_block1=""
code_block2=""
for row in rows:
    if field1 < field2:
        arcpy.CalculateField_management(fc,
                            newfield,
                            expression1,
                            expression_type,
                            code_block1)
    else:
        arcpy.CalculateField_management(fc,
                            newfield,
                            expression2,
                            expression_type,
                            code_block2)
    rows.UpdateCursor(row)
del row
del rows


texts: PostGIS error while adding feature: no result buffer

texts:
Can anyone explain to me the following error and possible fixes?  I am using QGIS 2.12.3 to connect to a PostgreSQL 9.5 database.

What confuses me is that I was able to add a feature and save, but then the second feature I tried to add threw this error.  I've had this happen on previous days, only with QGIS 2.12.  When I use 2.8 I do not get this problem (but some other issues with default values arise with 2.8 so I would prefer to use 2.12)

Any clues as to what's going on?

Update for more detail:
I have a postgis table of streets, and another table of signs (mostly freeway navigation signs) which are used for a navigation app.  The signs table does not have any geometry and is in the format (from_link, to_link, navigation instructions, sign_ID).  The from/to links are both ID numbers from the streets layer.

I never have any issue with adding features to the streets layer, only on the signs layer.  This is consistent for more than one layer (I have a separate table for each city) which leads me to believe there may be something weird going on with the lack of geometry?

Only constraint on the sign table is a primary key on the sign_ID

And here is the SQL used for creating the table.  You can ignore most of the fields, they are just various instructions for our navigation code:

CREATE TABLE public.austin_signs
(
  "Source_LinkID" bigint,
  "Destination_LinkID" bigint,
  "Sign_ID" bigint DEFAULT 0,
  "Seq_Number" character varying(255),
  "Exit_Number" character varying(255),
  "Exit_LangCd" character varying(255),
  "ExitNum_Tr" character varying(255),
  "Alt_Ex_Num" character varying(255),
  "Branch_RouteID" character varying(255),
  "BrRte_Tr" character varying(255),
  "Branch_RouteDir" character varying(255),
  "Sign_TextType" character varying(255),
  "Sign_Text" character varying(255),
  "Lang_Code" character varying(255),
  "SignTxt_Tr" character varying(255),
  "Trans_Type" character varying(255),
  "Toward_RouteID" character varying(255),
  "TowRte_Tr" character varying(255),
  "Straight_On" character varying(255),
  id integer NOT NULL DEFAULT nextval('austin_signs_id_seq'::regclass),
  CONSTRAINT austin_signs_pkey PRIMARY KEY (id)
)
WITH (
  OIDS=FALSE
);

I am the owner of the database so I don't think it is a permissions issue


texts: GEE - Reprojecting before exporting causes "User memory exceeded" error

texts:
I am studying ecosystem distributions in Norway using satellite data and figured that the Google Earth Engine is a very viable tool for this purpose. I am well aware that processing data with a high spatial resolution for a large (or should I rather say: elongated) country like Norway is a computationally demanding task, but I really hope there is room for improvement in my inexperienced workflow!
Details:

Calculate a number of spectral indices using a multi-year cloud-free composite from SENTINEL-2 surface reflectance and export the data (both indices + SR) to Google Drive
This is building upon an existing project which has a large number of prepared environmental rasters in a UTM CRS, plus, the training data has some legal restrictions which is why I cannot upload it directly to GEE
--> Need to reproject the satellite data to match the CRS, mask the exact spatial extent of the data, and merge it into one big GEO TIFF.
Exporting in the original maps Mercator CRS works like a charm and takes less than an hour
However, merging tiles, reprojecting, and masking on my local machine causes a lot of troubles and errors, have been experimenting with SAGA GIS, QGIS / GDAL (gdalwarp), and the raster package in RStudio, all of those solutions take many days of processing or crash for various reasons
--> Would really like to do the masking and reprojecting directly in the GEE code editor
Uploaded one of the environmental rasters to use as a mask after reprojecting using the GEE reproject() function

Here is an excerpt of the code:
function maskS2clouds(image) {...}

function addNDVI(image) {
  var ndvi = image.normalizedDifference(['B8', 'B4']).rename('NDVI');
  return image.addBands(ndvi);
}
...
var collection = ee.ImageCollection('COPERNICUS/S2_SR')   // Retrieve Sentinel 2 collection
    .filter(date)                                         // Filter date
    .filterBounds(geometry)                               // Filter tiles
    .filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))  // Prefilter clouds
    .map(maskS2clouds)                                    // Apply additional cloud filtering
    .map(addNDVI)                                         // Add NDVI
...

// Retrieve median values
var median = collection.median();

// Subset indices
var indices = median.select(['NDVI',...]); // Total of 5 indices

// Reproject to UTM
var proj = ee.Projection('EPSG:32633');
var indices_utm = indices.reproject(proj);

// Use raster mask
var indices_utm_nor = indices_utm.updateMask(ImportedRasterTif);

// Export to Google Drive
Export.image.toDrive({
  image: indices_utm_nor,
  description: 'Median_indices',
  scale: 30,
  folder: 'Median_indices',
  region: indices_utm_nor,
  maxPixels: 1e12
});

However, when applying the reproject() function before exporting, this results in a User memory limit exceeded. error. Any tips on how I could circumvent this in my case or other suggestions for doing everything locally?
Thanks in advance!


texts: KML layer missing metadata in ArcGIS Pro

texts:
Using ArcGIS Pro 2.7
I created a KML file from a raster with the Layer to KML tool in ArcGIS Pro. The KML layer does not show any metadata description in the Catalog View. The input raster was created in the same map and does show a complete metadata description, so I am not sure why the corresponding KML file does not have its own metadata description. (The raster layer metadata was auto-populated with a spatial reference and geoprocessing history, and the metadata for the KML layer did not). Why could this be?
I am wondering if the issue is to do with the KML as a xml file format.
Blank metadata of the KML layer from the Catalog View:


texts: I am carrying out some work, and my client requires isochrones generated between 1 and 9 hours (daily driving limits in EU). I'm scratching my head at the moment trying to find something that could generate this larger drive time, is this something I can do in QGIS? Esri Drive time analysis maxes out at 5 hours. There are limitations to most services. Am I missing something obvious?


texts:
I am on Ubuntu 14.04. When I am installing PostGIS-2.1, it's asking for libgdal1 library. But when I am installing libgdal1, it's removing grass and gdal.

Can someone please guide me how to install grass, gdal and postgis?

Thanks

Zia


texts: Call Toolbox operators of ENVI from IDL code

texts: Is it possible to call the clustering/classification algorithms included in the toolbox of ENVI (gui) from IDL (source code) and how?


texts: Check if point is near a linestring Spatialite

texts:
I have a LINESTRING made of two points start point and end point in spatialite.

UPDATE points SET road_geometry = GeomFromText('LINESTRING(24.889142 67.196188,
24.889673 67.196579)',4326) where name='name';

and Now, I want to check that if any other point is within that line string or how far the point is from the LINESTRING.


texts: gdal_translate and openweathermap xyz url

texts:
I'm trying to fetch OpenWeatherMap GeoTIFF layers using gdal_translate from the website

https://tile.openweathermap.org/map/temp_new/{z}/{x}/{y}.png?appid=API_KEY

I used this command:

gdal_translate https://tile.openweathermap.org/map/temp_new/{z}/{x}/{y}.png?appid=API_KEY geo.tif

but I got this answer:

ERROR 1: An error occurred while creating a virtual connection to the
DAP server: Error while reading the URL:
https://tile.openweathermap.org/map/temp_new/{z}/{x}/{y}.png.ver.
The OPeNDAP server returned the following message:
Unauthorized: Contact the server administrator.`

How do I fetch it?


texts: EE Spatial Join-- Generate columns in primary collection

texts:
I have two feature collections- a feature collection of different polygons (primary collection) and a .csv with various location data (secondary collection). I would like to append each column in the secondary collection to the each feature in the primary collection.

As it is now, I have a list of features (columns in secondary collection) appended to each feature in the primary collection, but when I export the data to a csv, all that shows up is a column with <Feature> in each column.

Here is my "secondary" feature collection. I would like to have these columns in the output csv:

Here is a snippet of the output CSV. See that "points" is the name of the secondary feature collection, and information such as Species, and No_Fish is lost.

Here is my pseudo-code solution, extracting properties from an already joined feature collection:

for feature in joined:

    feat_points = feature.get("Points")

    //Iterate over each point in secondary collection that has been joined
    for point in feat_points:

        //Set each property name in secondary collection to primary collection
        for property in point:
            feature.set("Property_name",property)

Here is a link to my code:
https://code.earthengine.google.com/a3ff4b90ac87ccc9199ad4861a9d4f54


texts: Python crashes after reading a certain number of mxd files

texts:
In one of my work, I need to go through about a few thousand mxd files, get some information from each MapDocument object (scale, printing paper size, map title, etc.), and write those information to a record (a csv file, for example).

My Python script run perfectly fine for the first 30 mxd files. When it reached the 31st mxd file, I got a Windows error saying "python.exe has stopped working. A problem caused the program to stop working correctly..." which gave me only 1 option to close the program. No error from the script itself. To be particular, the script failed at the line elms = arcpy.mapping.ListLayoutElements(...).

I run the script the 2nd time, this time I started from that 31st mxd file that I failed. The script run perfectly again for the next 97 mxd files. The same Windows error popped up, and again, Python crashed at the line elms = arcpy.mapping.ListLayoutElements(...).

I run the script the 3rd time, starting from the mxd file that the script failed to read. This time Python crashed after 49 mxd files, still at the line elms = arcpy.mapping.ListLayoutElements(...).

I repeated the above 3 run, and got exactly the same result, i.e. Python crashed after reading the exact number of mxd files, at exactly the same line.

I guess it's not the problem of the mxd file. It's probably arcpy itself or the way Windows and Python (or arcpy?) work together.

I then wrote a simple test script with just one for loop that repeatedly opens one of my mxd, and calls that line elms = arcpy.mapping.ListLayoutElements(...). Python crashed after 73-77 loops. I also tried with arcpy.mapping.ListDataFrames(mxd)[0] and mxd.pageSize, and they all made Python crash after 73-77 loops. Interestingly, some arcpy methods (e.g. mxd.dateSaved, mxd.description) did not cause the crash at all.

Anyone else has experienced this? I'm using Standard ArcGIS 10.5.

PART OF MY SCRIPT:

The last thing printed out before every Python crashes is the line print count, "mxd opened", so my best guess is that the next line elms = arcpy.mapping.ListLayoutElements(mxd, "TEXT_ELEMENT", "title") is what causing the crash.

def get_record(my_mxd):
    global count
    try:
        mxd = arcpy.mapping.MapDocument(my_mxd)
        count += 1
        print count, "mxd opened"
        elms = arcpy.mapping.ListLayoutElements(mxd, "TEXT_ELEMENT", "title")
        print "Element listed"
        if len(elms) == 0:
            fail_writer.writerow([my_mxd])
            return
        ...
    except Exception as e:
        print e
        fail_writer.writerow([my_mxd])

UPDATE 2

After several crashes and re-run the script starting from the crashed MXD, I finished reading ~900 MXDs and got another weird behavior. Instead of crashing, now Python gave me the error Object: CreateObject cannot open map document which was caught by try except, all the way to ~5000 MXDs.

But again, when I re-run the script, made it start from the first MXD that got this error, it run perfectly for the next ~100 MXDs, then throw the error Object: CreateObject cannot open map document all the way to ~5000 MXDs.

Looking forward to what I will get after this :)).

UPDATE 3

After ~2000 MXDs, the crash came back. This time I got some consistent crashes with several MXDs, i.e. re-run the script from the failed MXD kept crashing Python.

The fun part: some consistent MXDs crashed Python after pushing memory usage to 1GB, others crashed Python with just 170MB memory usage. At least Python struggled with those MXDs a few minutes before being knocked out. There's only one MXD that beat Python in 3s with just 170MB. All these MXDs can be opened in ArcMap without any problem.


texts: gdalsrsinfo conflict between command line and python console

texts:
I am working on a PyQgis application which loads WRF model output as a raster to overlay on a map which is in the Pseudo-Mercator projection.  The WRF data is projected to Lambert-Conformal and am using PyProj to retrieve the projection metadata and GDAL to transform the original NetCDF file to GeoTiff.  Initially I was getting issues with "gdal_translate" not being able to decipher the Proj4 string I produce, so I decided to use the GDAL API directly to debug; this required converting the Proj4 string to an OGC WKT.  I use the QProcess class in the application to run the following:
gdalsrsinfo -o wkt ${MY_PROJ_STR}

where "${MY_PROJ_STR}" is:
"+proj=lcc +lat_1=13.342499732971191 +lat_2=47.317501068115234 +lat_0=30.330001831054688 +lon_0=-86.16000366210938 +x_0=0+y_0=0 +a=6370000 +datum=WGS84 +no_defs"

When I run the "gdalsrsinfo" command from the application I get the following stderr output:
ERROR 1: ERROR - failed to load SRS definition from "+proj=lcc +lat_1=13.342499732971191 +lat_2=47.317501068115234 +lat_0=30.330001831054688 +lon_0=-86.16000366210938 +x_0=0+y_0=0 +a=6370000 +datum=WGS84 +no_defs"

yet when testing it from the command line, it get the WKT translation:
PROJCS["unnamed",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Lambert_Conformal_Conic_2SP"],PARAMETER["standard_parallel_1",13.34249973297119],PARAMETER["standard_parallel_2",47.31750106811523],PARAMETER["latitude_of_origin",30.33000183105469],PARAMETER["central_meridian",-86.16000366210938],PARAMETER["false_easting",0],PARAMETER["false_northing",0]]

I'm not sure if it's the PROJCS="unnamed" part that's throwing it off.  I checked to make sure the versions of GDAL used in Python and from the command line weren't in conflict, and they seem to be the same:
print(gdal.VersionInfo())
>> 2030200

gdalinfo --version
]$ GDAL 2.3.2, released 2018/09/21

My other system specs for reference:  am using Python 3.7, PyQt5 5.13.2, pyproj 3.0.0.post1, PROJ 5.2.0, and QGIS 5.10.  The operating system is Fedora 31. I'm not sure if the syntax of the PROJ string is correct or if I am missing a step in producing the correct reference system. I've also added "GDAL_DATA" and "PROJ_LIB" environmental variables to my ~/.bash_profile.


texts: Open multiple files at once by QGIS actions buttons

texts:
I'm planning to put into my table a column with folder path, so it will be easy to open a folder with needed files in attributes table. Now I'm wondering if I could open not the folder, but all files in it.

It would be great to just get all pdf in directory, no matter how they are named (so I don't have to store filenames anywhere). But the schema is simple: all files have prefix and number from _1 to _7, eg. a_1.pdf, a_2.pdf, a_3.pdf... But sometimes one of those files is missing.

I tried to put all paths to varchar, as string like "C:\a_1.pdf;C:\a_2.pdf(...)" or as array but "open" doesn't work with more than one path. Any suggestions?


texts:
I have gis data which i can access over Arcgis Rest with Arcgis js api and i have non-gis data which i get from different db.I want to mashup this two like in this link

http://www.cartographicperspectives.org/index.php/journal/article/viewFile/cp76-hunt/1311  (look page 6)

this is my code

DemMin = new Number; DemMax = new Number;
jQuery("#testMenu").click(function () {
// sample -- > is feature layer already defined somewhere else in the code

    $.getJSON('/home/getresult', null,  // this returns array of json like [{grid_no= 2,avg_heat=3},{grid_no= 1,avg_heat=5}]
function (data) {
   $.each(data, function (i, C) {

       if (C.avg_heat < DemMin) { DemMin = C.avg_heat; }
       if (C.avg_heat > DemMax) { DemMax = C.avg_heat; }

       $.each(sample.graphics, function (k, G) {

           if (G.attributes.GRID_NO == C.grid_no) {
               G.attributes.AVG_HEAT = C.avg_heat;  // (1) adding featurlayer.graphics.attributes not featurelayer field i know ::)),  look for (2)
               alert(G.attributes.AVG_HEAT);   // this alert works and shows G.attributes.AVG_HEAT which i added
               return false;
           }
       });
   });

   var testrenderer = new SimpleRenderer(new SimpleFillSymbol().setOutline(new SimpleLineSymbol().setWidth(0.5)));

   testrenderer.setColorInfo({//set classification values and symbology color information
       field: "AVG_HEAT",   // (2)  This my adding, why i getting error ,this is field not in feature layer i know:))
       minDataValue: DemMin,
       maxDataValue: DemMax,
       colors: [
           new Color([255, 0, 0]),
           new Color([255, 255, 0])]
   });

   sample.setRenderer(testrenderer);

   sample.redraw();

});

the error i get of course

--  esri.layers.FeatureLayer: unable to find 'AVG_HEAT' field in the layer 'fields' information  "feature_layer_url"

is there anyway to this code working?

Note : I know Arcgis js api  "Gas price data by state mashup"  sample is similar mashup but please stay on this sample !!


texts: Problems with QGIS (v 1.8.0) layer labels

texts:
I have a PostGIS (v 2.0) layer representing lakes. The attributes "chinchar" and "chintrans" contain Chinese characters and Pinyin (Roman) characters in UTF8 format. When I use Layer -> Labeling in QGIS, this works very well to name the lakes in about 90% of the cases, but the remaining 10% refuse to display a label. I am using the following string for "Label this layer with...":

"chinchar" || '\n' || "chintrans"

The new labeling system seems to be very good and I would hate to go back to the "old" labeling techniques of defining a separate column in the table or defining a view in PostGIS.
Thanks for any ideas...


texts: Automatic water extent mapping from Sentinel-1 SAR using Otsu's method

texts:
var imgVV = ee.ImageCollection('COPERNICUS/S1_GRD')
        .filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VV'))
        .filter(ee.Filter.eq('instrumentMode', 'IW'))
        .select('VV')
        .map(function(image) {
          var edge = image.lt(-30.0);
          var maskedImage = image.mask().and(edge.not());
          return image.updateMask(maskedImage);
        });

var desc = imgVV.filter(ee.Filter.eq('orbitProperties_pass', 'DESCENDING'));
var asc = imgVV.filter(ee.Filter.eq('orbitProperties_pass', 'ASCENDING'));

var spring = ee.Filter.date('2019-03-01', '2019-04-20');
var lateSpring = ee.Filter.date('2019-04-21', '2019-06-10');
var summer = ee.Filter.date('2019-06-11', '2019-08-31');

var descChange = ee.Image.cat(
        desc.filter(spring).mean(),
        desc.filter(lateSpring).mean(),
        desc.filter(summer).mean());

var ascChange = ee.Image.cat(
        asc.filter(spring).mean(),
        asc.filter(lateSpring).mean(),
        asc.filter(summer).mean());

Map.setCenter(83.77193293042244,21.595416478393343, 12)
Map.addLayer(ascChange, {min: -25, max: 5}, 'Multi-T Mean ASC', true);
Map.addLayer(descChange, {min: -25, max: 5}, 'Multi-T Mean DESC', true);
print(imgVV);

Got following error
ImageCollection (Error)
Collection query aborted after accumulating over 5000 elements.


texts: Seeking Strava users routes on city scale?

texts:
I'm a GIS student and I would like to access all the routes used by the users of Strava at the scale of a city (Rennes, France) in SHP, GEOJSON or JSON format.

Is it possible to make requests by city or geographic coordinates on the Strava API?

Is there an application "for dummies" where you could download all the data inside a polygon?

I use QGIS if it's relevant and I have almost no knowledge of using APIs in general. I don't even know if the data is free, I try to contact Strava METRO but they didn't answered.


texts: Meaning of 'bnds' and dimensions in OCO2 v9 L3 data product

texts:
I cannot seem to understand what, for example, 'lat_bnds' means? Same go for 'lon_bnds', lat/lon_dimensions. These values come from the structurally flat NCDF4, sourced from the JPL as a custom product.

Reading all available documentation does not hint of any definitions. Here's a visual explanation of my data:


texts:
I would like to explain my situation.

I have a GNSS receiver, which outputs me NMEA data (read below).
I receive RTCM data via NTRIP from a caster.
I would like to correct my GPS position with centimeter-accuracy.
I do not want to use factory-made GNSS devices which perform this correction on their own.

I know (from other questions on GIS_STACKEXCHANGE) that with only NMEA data I cannot perform any such correction. So I need to buy a GNSS receiver which gives me raw data (is it right?).
Now some questions:

Do I really need these raw data?
Which kind of data need I own? Ephemeris? Pseudoranges? What else?
Do I need to get these same data from the caster? I studied the format and I got that there are several message codes each of one refers to a particular kind of data. Which message need I look to?


texts: Join two tables completely on geometry (postgis)

texts:
I have two tables with the same street network as geometry (but not a same id).
What I want to do is to join both tables with all attributes to one for all geometries which overlap exactly.

Unfortunately my query just gives me emtpy lines back. I am working with left join, but all other joins won't work as well.

DROP TABLE IF EXISTS test;
CREATE TABLE test AS
SELECT *
FROM  table1
RIGHT JOIN table2 ON ST_Equals(table1.the:geom, table2.the_geom);

Does anybody know an answer for this?


texts: How do you obtain the execution times of processes in qgis?

texts: I am working on QGIS. I need to determine how much time QGIS needed to finish a specific task. I have been unable to get execution times for the processes. What is the proper way to get processing time under QGIS?


texts: Count number of days by year and by month where temperature fall below 0 celcius degree

texts:
I want to generate raster for each month , where the value of the pixel refers to the number of days per month and per year where temperature fall under 0 degree celsius. In this example the year range is [2010,2012], which means that i will get 36 images. But what i get , is 36 images with no band.
Any solution ?
var years = ee.List.sequence(2010,2012);
var months = ee.List.sequence(1,12);
var imcol = ee.ImageCollection("MODIS/006/MOD11A1");
var opla = function(image){ return image.select('LST_Day_1km').multiply(0.02).subtract(273.15).lte(0).clip(table)};
var filtred = imcol.map(opla);

var sommation_frost = function(y){
  var sommar = filtred.filter(ee.Filter.calendarRange(y,y,'year'));
  var montar = function(m){
    var x = sommar.filter(ee.Filter.calendarRange(m,m,'month')).reduce(ee.Reducer.sum())
    return x
  }
  return months.map(montar)
}

var sortie = years.map(sommation_frost)


texts:
I'm trying to do something really simple, plot a rectangle using psxy on a .cdf grid file. Something like:

grdimage grid.cdf -JM20c -R-70/0/20/60 -Bf20g10 -K > image.ps
psxy rectangle.dat -JM -R -B -W0.5p -O >> image.ps

My rectanle.dat file is a two column file with 5 lon/lat points. Once plotted, the rectangle has two straight vertical sides and two curved upwards horizontal sides and I cannot figure out why. Any ideas?


texts:
I am trying to create a raster from a polygon (ie polygon to raster conversion tool) in ArcGIS 10.3 but I receive a 999999 ERROR message. The raster is at 1km from a 1km fishnet/mesh polygon for the entire continent of Africa. When I try to create this same raster from a subset of the polygon (ie I clipped out a small section of the 1 km fishnet/mesh polygon and used the same methodology to create the raster), it worked well with no issues. However, when I try to create it for the larger area of all of Africa, it fails.

I have tried to trouble shoot - it is in a personal geodatabase (not a shapefile) so I would hope that file size would not be an issue, I have unchecked the "background processing" option in Geoprocessing > Geoprocessing options, and I have even tried restarting ArcGIS just in case it had some "hang up" from creating the fishnet (I have tried in catalog too).

Does anyone have any ideas as to why this might be happening?

Maybe it is just that I don't have a large enough processor to handle the 1km raster for the entire continent of Africa?


texts: Modifying default properties for HTML Annotation in QGIS?

texts:
How do I modify the default properties for the HTML annotation in QGIS 2.18?

For example, I don't want the red spot as a symbol or I want to change the frame width from 1.00 to 0.2.


texts: Grouping of layer with GeoServer and OpenLayers

texts:
I need help to display layers group/non group
1. with checkbox to select layer or group
2. icon/symbol for each layer.

I saw example codes from OpenLayers where it's adding the layers statically. I would like to populate it dynamically so in future if a layer is added/remove it will be displayed without changing any code.

Sample layout of what I want to achieve is posted below, please note that checkboxes are displayed as dots.

    Groping of Layer

                Layer

                        School

                        Road

                                Metal Road

                                Unmetal Road

                        Parecls

                        Water Feature

                        Cities

                        Disulict

                        India


texts:
I cannot find the AutoForm plugin. I am using QGIS 3.6, have this plugin perhaps given another name, or where do I find the functionality.

I am trying to create lookup tables following this tutorial
https://blog.sourcepole.ch/2017/02/10/autoform-plugin-for-qgis/


texts: Create new string field and replace numbers from other field with text

texts:
I would like to create a string field in a shapefile where the string expression is dependent on the numbers given in another field. For example if there is a '1' in the source field there should be the expression 'Kontrolle' in the string field; if there is a '2' in the source field there should be the expression 'Kontrolle Baumkolonne' and so on.
Can I realize this in the field calculator?


texts: Create contour line from point table in PostGIS

texts:
I have a point table containing values ranging from 0 to 100 and I would like to draw a contour line for the value 70. (See image below)

I know how to create contour lines in QGIS for raster data but now I would like to do it with vector data (my point table) using PostGIS so that it can be integrated in a python script.

Any idea on how to start?

EDIT:
This question is different from Generating contour lines and filled contours from points as it doesn't involve raster data (neither in the input nor in the output)


texts: I have county-wide parcels that don't match up. They are in the same feature dataset, have same projection etc. Is there a tool to get the borders/parcels to match? One layer of parcels seems to be more rough and not as accurate but I at least want the borders to match up.


texts: Issue on Loading Required ClusterLayer.js With dojo

texts:
I have js file like this and I have saved the ClusterLayer.js in all directories like
root --> ClusterLayer.js
root --> extras/ClusterLayer.js
root -->js/extras/ClusterLayer.js
root -->js/ClusterLayer.js

and the file is like
    require([
    "esri/map",
    "dojo/_base/array",
    "esri/layers/ArcGISDynamicMapServiceLayer",
    "esri/geometry/Geometry",
    "esri/geometry/Point",
    "esri/geometry/webMercatorUtils",
    "extras/ClusterLayer",
    "esri/graphic",
    "esri/symbols/SimpleMarkerSymbol",
    "esri/symbols/SimpleLineSymbol",
    "esri/symbols/SimpleFillSymbol",
    "esri/Color",
    "esri/InfoTemplate",
    "esri/dijit/HomeButton",
    "esri/dijit/OverviewMap",
    "dojo/parser",
    "esri/layers/GraphicsLayer",
    "esri/SpatialReference",
    "esri/dijit/PopupTemplate",

    "dijit/layout/BorderContainer",
    "dijit/layout/ContentPane",
    "dojo/domReady!"
    ], function(Map, arrayUtils, ArcGISDynamicMapServiceLayer, Geometry, Point, webMercatorUtils, ClusterLayer, Graphic, SimpleMarkerSymbol, SimpleLineSymbol, SimpleFillSymbol, Color, InfoTemplate, HomeButton, OverviewMap,   parser, GraphicsLayer, SpatialReference, PopupTemplate)

but still getting error

Can you please let me know why this is happening?


texts:
Trying to use CQL bounding-box filter, but get a strange result like this:

loadFeatures({"type":"FeatureCollection","totalFeatures":15640,"features":[]})

It showing a total 15640 features in that bbox, but why there are no list of them and "features" are just empty list?

This is my CQL filter:

http://loom-gis.geo.ut.ee:8080/geoserver/ermas/ows?service=WFS&version=1.1.0&request=GetFeature&typeName=ermas:testdata_geopnt&outputFormat=text/javascript&format_options=callback:loadFeatures&maxFeatures=50&srsname=EPSG:3301&EPSG:3301&CQL_FILTER=%28BBOX%28geometry,494981.10875,%206437855.3625,%20684824.8587499999,%206569105.3625%29%29


texts: Putting highcharts chart in OpenLayers popup?

texts:
I build a openlayers map on my website.
and there are many feature vector on the map.

now I want to put a highcharts charts in popup.
to show the temperature in every city
something Like this

If it is possible, then how do I do it?

I used contentHTML but it is failed

highcharts code

$('#chart').highcharts({
    chart: {
        type: 'area'
    },
    title: {
        text: 'Area chart with negative values'
    },
    xAxis: {
        categories: ['Apples', 'Oranges', 'Pears', 'Grapes', 'Bananas']
    },
    credits: {
        enabled: false
    },
    series: [{
        name: 'John',
        data: [5, 3, 4, 7, 2]
    }]
});

part of openlayers code

var content = "

" ;
        var popup = new OpenLayers.Popup.FramedCloud("popup",
            OpenLayers.LonLat.fromString(feature.geometry.toShortString()),
            null,
            content,
            null,
            true
        );


texts:
I am trying to use time manager to show patterns of movement for tagged wildlife.  I have a GPS point with time and date every two hours.

Then, I select the layer with time manager, and select the date and time column (date and t) in the example and confirm the time format mirrors the information layout of the feature.

I have selected the time frame size to be 2 hours, however, time manager picks up every feature on the same hour even if they are from different dates.  p.s. the same happens if time frame size is changed to date.

I hope the explanation is clear.  Looking forward to your comments/suggestions.  hopefully not missing something blatantly obvious.

Second attempt:

i have updated the attribute field format to reflect date and time format as follows:

I then add the layer to time manager and add this format to the time format:

the result is as follows:


texts:
I need to know the type of unit that Zoom, or maxZoom/minZoom has.
I've been finding in OpenLayers API but without success.

Do you have any idea?


texts: Open Trip Planner - Error no graph registered with routerID

texts:
I am new to using OTP. I followed the basic tutorials instructions to start OTP, with own data, as well as with the data indicated within the tutorial. I get it to run the grizzly server, but opening local host shows no results and in my Windows command line the ERROR message "no graph registered with router ID" occurs. Further when I want to switch to the analysis window in OTP another error message occurs saying something like "Not possible, out of bounding box". My folder structure is

  D:\Program\Portland\otp

for the jar-file and the inpput (gtfs.zip &osm.pbf) data is stored under

  D:\Program\Portland\otp\graphs

my command to launch otp reads as following:

D:\Program\Portland\otp> java -Xmx1G -jar otp-0.19.0-shaded.jar --build --inMemory --analyst

I checked all solutions offered found via google and the result kept the same. Anyone an idea what myproblem could be?


texts: Why does zoom not work for one of the openlayers base layers?

texts:
I have map with 2 base layers (osm, bing) and few overlays. I put layerswitcher outside of map in separate div.
I wrote script for toggle layers in that div and it works fine. But there is one problem with base layers. If I select bing as a base layer in my div it displays fine but I can't zoom in or out. Point features (ones from the overlay) are moving but base layer is static.
Here is my code:

BASE LAYERS

              <input type="radio" name="baselayer" value="osm" id="IDosm_layer" onclick="toggleControl(this);" autocomplete="off" checked>
              OpenStreetMap

              <input type="radio" name="baselayer" value="bing" id="IDbing_layer" onclick="toggleControl(this);" autocomplete="off">
              Bing Aerial

OVERLAYS

              <input type="checkbox" value="2" id="IDcurrent_position" onclick="toggleControl(this);" autocomplete="off">
              Overlay1

That's the html part of code for layer tree div and here is javascript code:

function toggleControl(element){

     if(element.value == "osm"){
      bing.setVisibility(false);
      osm.setVisibility(true);
     }
     else if(element.value == "bing"){
      osm.setVisibility(false);
      bing.setVisibility(true);
     }

     for(i=2; i<12; i++) {
        if(element.value == i) {
            map.layers[i].setVisibility(element.checked);
        }
    }

}

I hope you understood my problem and I'm looking forward to your answer.


texts:
Im trying to add StreetView functionality to my OpenLayers web app, and would like to add the overlay showing StreetView availability to the map itself.  I have seen an example of it being done successfully here, but i cannot get the code to work in my own solution.  Specifically :

var layer = new OpenLayers.Layer.Google("Google Street", {sphericalMercator: true});

layer.mapObject.addOverlay(new GStreetviewOverlay());

results in the error:

layer.mapObject is null

Has anyone got any idea what i might be doing wrong, and if so, how can I fix this?


texts: How to create proj4.js and proj4leaflet.js for coordinate transformation

texts:
I need to know how to create a proj4.js and proj4leaflet.js which requires me to convert coordinates from EPSG:31700 to wgs84.

Both scripts I will load in leaflet script.


texts: Tabulate Area 2 Tool not returning correct area

texts:
I am having trouble with the output of the "Tabulate Area 2" tool in the following update:

http://blogs.esri.com/esri/arcgis/2013/11/26/new-spatial-analyst-supplemental-tools-v1-3/

I have a layer of overlapping buffers (polygon) as the zones and then a layer of parks where I'm trying to find the area of parks (also polygon) that intersect with each individual buffer. The tabulate area 2 tool seems like it should be exactly what I want. However, the area of intersection it's returning is not what it should be and I have no idea why. I've gone ahead and projected both files so they're the same and in meters, and I've also converted one of my polygon files from a multi-part to a single-part. The numbers are off by thousands of meters. Does anyone have any ideas?

The highlighted FIDs at the top should have values of 87618, 23323, 21837, and 280335 meters, respectively.


texts: How to prepare and upload points/lines to GPS using QGIS 2.8

texts:
I would like to upload some points (e.g. locations of transmission towers) and lines (e.g. streams, roads) in .shp format to my Garmin Oregon 600, using QGIS 2.8.2.

Could anyone provide with step-by-step instructions?


texts: My First Custom Widget (WIP) : Rubber Band Zoom

texts:
I have created a new widget and called it "MagnifyingGlass" (at the request of the client). I have copied from the WAB (Web App Builder) the Sample Widget (the SimpleWidget folder structure for the 'Get Map Id' widget example).

I have renamed everything.
I have added it to the site config file so that it shows up.
It is now a visible button on one of my dev maps.
It is the duplicate BOOKMARK widget button to the right of the buttons on the bottom left.
If you hover over it is says 'Zoom.'

The problem is what happens when you click it!
It opens a window / popup...

The desired functionality is simply to activate this javascript command "this.map.enableRubberBandZoom();" and nothing more.

The whole goal for this widget is to simply add a button that users can click that lets them THEN click and drag to draw a Zoom window.

Apparently, this is functionality that SHOULD be (obviously) on all maps. SO...
In my Widget.js file (remember I copied the Simple Widget example files) I have commented out the templateString, rewrote _getMapID to now be _init simply because I don't know what I'm doing (and I'm not sure where to look).

The whole goal is the Rubber Ban Zoom Box creation after clicking the Widget button. Period. But I am uncertain where and how, in the Widget.js file, to make this happen.

I have asked this on ESRI's GeoNet site already: https://geonet.esri.com/message/438305


texts: PyQGIS - Style Features individually based on attribute value - custom renderer as option?

texts:
I am trying to colorize all features in all layers in my QGIS 2.8.2-Document based on a color-value, that is being determined for every Feature individually.
I do have a QtGui.QColor instance for every feature, but I still need to apply it to the Feature via a Renderer.

from PyQt4 import QtGui, QtCore
canvas = qgis.utils.iface.mapCanvas()
layers = canvas.layers()

for layer in layers:
    iter = layer.getFeatures()
    for feature in iter:
        if  feature['ELEM_NR'] != '' and feature['ELEM_NR'] != NULL:
            #
            color = QtGui.QColor()

            # ...
            # hexValue is being determined based on feature['ELEM_NR']
            # ...

            color = QtGui.QColor(hexValue )

            # how can I apply color to feature`?

Background (why I do not simply use graduated symbols):

Every Layer has an attribute that is kind of an id.
I want that every Feature (no matter which layer it is in), that has the same attribute-value, is colored the same way.

Because the the number of features differs between layers and the value range of the attribute is quit big, neither the catgorized nor the graduatedSymbol-Renderer descibed here would bring the desired result (at least in my understanding).

Using a custom renderer as described in the documentation looks promising, but I can't figure how to include this.

Any Idea?


texts: ArcGIS JavaScript Print task not working in Internet Explorer 10. error: access is denied

texts:
I have created a custom print task using my own custom templates and the Export Web Map tool in Server Tools of ArcGIS 10.2 to create a geoprocessing service similar to the print task made by Esri.

The geoprocessing service sits in the same directory as my map service for the mapping application I am working on.  In Internet Explorer 10, my map application works perfectly fine - but when I click on the Print button to print I get an error in the console stating "Access is denied".

I ended up adding a proxy to my application after getting the following error in the console when using the print button originally: "Error: esri.config.defaults.io.proxyUrl is not set.  IF making a request to a CORS enabled server, please push the domain into esri.confi.defaults.io.corsEnabledServers".

So I added a proxy into my application and then got the "Access is denied" error.  I tried the same functionality with Esri's out-of-the-box print task but had the same error.

Does anyone know why this is happening or how to fix it?  My print button works in both Firefox and Chrome with no problems, and I thought by adding the proxy the problem would be resolved.


texts: Reloading layer in Layers Panel with Python after modifying it in DB Manager

texts: I often have to reload layers in the layers panel after doing changes in the DB-Manager (add a column, delete a default value, ...). At the moment I use the changeDatasource plugin. Can this also be done using Python?


texts: Convert MIF and MID files into postgis tables

texts:
I have tables like this:

CREATE TABLE postcode_area_boundaries
(
    ogc_fid INTEGER,
    wkb_geometry GEOMETRY(MultiPolygon,3857),
    postarea VARCHAR(2),
    area_km DOUBLE PRECISION,
    area_mi DOUBLE PRECISION,
);

CREATE TABLE postcode_area_points
(
    ogc_fid INTEGER PRIMARY KEY NOT NULL,
    wkb_geometry GEOMETRY(Point,3857),
    postcodearea VARCHAR(254),
    postcodeareaname VARCHAR(254)
);

I have also files:

Postcode Area Boundaries.MID

"AB","Aberdeen"
POLYGON((383707 866131,383465 866240,383302 866387,383202 866505,383089 866586,.....))
"AL","St Albans"
POLYGON((526248 203348,....))
.......

Postcode Area Boundaries.MIF

VERSION 300
DELIMITER ","
COORDSYS EARTH PROJECTION 8,79,"m",-2,49,0.99960127,400000,-100000
COLUMNS 2
  POSTCODEAREA CHAR (0)
  POSTCODEAREANAME CHAR (0)
DATA
POINT 384222.000 819125.000
POINT 518661.000 209809.000
POINT 407763.000 285514.000
POINT 371106.000 147748.000
.... as many points as many was areas in MID file (about 126)

I probably should use ogr2ogr but have no idea how. I am trying with commands like this:

ogr2ogr -f PostgreSQL PG:"<dbconn>" "Postcode Area Boundaries.MIF" -t_srs EPSG:3857 -nlt MULTIPOLYGON -overwrite -nln postcode_area_boundaries

But I have errors like this:

  Warning 1: Geometry to be inserted is of type Point, whereas the layer
  geometry type is Multi Polygon. Insertion is likely to fail ERROR 1:
  INSERT command for new feature failed. ERROR:  Geometry type (Point)
  does not match column type (MultiPolygon)

Can anyone help me with this?


texts: realtime maps for coastal zone using arcmap10 with added layers

texts:
i would like to make the best mapping serices for clients for exact mapping layers for High tide line and low tide line ( which the data is being collected physically in x,y format) with this data should i make maps in arcmap10 with desired layers and desired scale for presenting the sylish maps.

if yes kindly be send each step by step procedure

-jawaharkumar
jawahar_ch@yahoo.com


texts: Getting MemoryError from Fiona when merging large amounts of buffered data (Python)

texts:
I am using fiona and shapely together to buffer points and then merge them and write them to a shapefile.  For almost all of my datasets, this is working fine.  However, I am currently getting a MemoryError on the mapping(merged) call.

Here's the error that I'm getting:

  Traceback (most recent call last):
    File "C:\Working\Projects\USACE_Data\Toolbox\V4.1\Tools\ImportEHydro.py", line 149, in execute
      helper.ExcludeNewerDataWithinArea(lyrUSACE, fileTemp, fcSurveyJob, long(sordat.strftime("%Y%m%d")))
    File "C:\Working\Projects\USACE_Data\Toolbox\V4.1\Tools\HelperClass.py", line 92, in ExcludeNewerDataWithinArea
      self.USACEMinimumBoundingPolygon(usaceLayer, newerDataPoly , "SORDAT")
    File "C:\Working\Projects\USACE_Data\Toolbox\V4.1\Tools\HelperClass.py", line 609, in USACEMinimumBoundingPolygon
      firstBufferFC = self.BufferFeatureClass(tempCopy, "new_buffer", mean + (4*std), mergeType, mergeField)
    File "C:\Working\Projects\USACE_Data\Toolbox\V4.1\Tools\HelperClass.py", line 540, in BufferFeatureClass
      'geometry': mapping(merged)
    File "C:\python\ArcGIS10.5\lib\site-packages\fiona\collection.py", line 353, in write
      self.writerecords([record])
    File "C:\python\ArcGIS10.5\lib\site-packages\fiona\collection.py", line 347, in writerecords
      self.session.writerecs(records, self)
    File "fiona\ogrext.pyx", line 1215, in fiona.ogrext.WritingSession.writerecs
  MemoryError

The relevant bits of code are below:

with collection(os.path.join(newShpDir, outputShp + ".shp"), "w", "ESRI Shapefile", schema) as output:
     # Merge all of the features into one feature
     for i in range(len(features)):
         try:
             merged = cascaded_union(features[i])
         except:
             arcpy.AddWarning("  Could not merge features for " + (str(featureList[i]) if len(featureList > 0) else "(unknown)") + "...skipping...")
             continue

         # write the output to a shapefile
         dictProperties = {}
         if mergeAll:
             dictProperties = {mergeField : str(bufferDist)}
         else:
             dictProperties = {mergeField : featureList[i]}

         output.write({
             'properties': dictProperties,
             'geometry': mapping(merged)
         })

The error is happening at the 'geometry':mapping(merged)  It's trying to merge about 350k shapes into one.  Anyone know if there is a limit to how many it can handle or if this is normal?  I can handle it with an exception and then union some, merge, union some, merge, etc...but I'd rather not if I don't have to...


texts: How to populate From Node and To Node fields for a polyline feature class? (ArcGIS 10)

texts:
I need to populate FNODE and TNODE fields for my roads dataset. I have an intersections point feature class that has unique ids, and I need each road segment to have an FNODE and TNODE value that corresponds to the intersection id. The roads are broken into segments at each intersection so that each segment only touches 2 intersection points. I thought that maybe Network Analyst generated these, but it doesn't look like it actually writes these values. Does anyone know how to generate these?

This dataset won't be used in Network Analyst - it will be used in a custom legacy add in that a client uses.


texts: Exporting geo-referenced collada from ArcGIS Pro to Sketchup?

texts:
Using the Multipatch to Collada conversion tool on ArcGIS Pro I have been able to export 3D models and open them in Sketchup maintaining their 3D features, however they lose their spatial reference and appear as individual entities floating around in space.
What can I do to keep the spatial reference of the exported 3D files in Sketchup?


texts:
I have data of one-dimensional points (locations in a straight line). I want to check if the cluster poisson process models or Cox models are a good fit for the data.

Since my data only has x-coordinate, I tried the linear network poisson process model. However, lppm supports only Poisson models.

The kppm methods, which have the cluster models, require two-dimensional data. So, I added a dummy y-column with zero values and range of (ymin=0,ymax=0.001). ymax=0 returns errors during computation. Now, I am able to fit Matern cluster, etc.

My question is what would be the best way to handle one-dimensional data?

Is adding a dummy column with a non-zero range the only solution? Or did I miss some details about point patterns or process models?

Suggestions for alternatives are also welcome.


texts: How do I show an aliased field in a CartoDB.js Infobox?

texts:
I used CartoDB.js to do a SQL query of a layer which also performs some calculations on the layer's fields. This part works fine, but I just can't figure out how to get the calculated field into my infobox template.

For example, in this sample code from the CartoDB github:
https://github.com/CartoDB/cartodb.js/blob/develop/examples/tutorials/tutorial-2.html
if I changed line 47 to be:

sql: "SELECT pop_other/10 AS testval FROM ..."

How would I show that "testval" in an infowindow? I've tried doing {{content.data.testval}} on line 20, but that doesn't work.


texts: Slow rendering of tin of large area in ArcMap?

texts:
i am using ArcGIS 10.3.1 and working on TIN generation for entire H.P(India).

I have downloaded  DEM & generated TIN but it is very slow and when I am doing anything the entire TIN starts forming again and it is taking a lot of time. I have tried many options like started working on geodatabase , done some system setting but still the problem exists.

My system is 8 GB RAM and also have dedicated & 2GB extended graphics. I have tried opening in ArcScene for 3D view but it is only processing and after sometime it stops.

I am unable to  do any analysis. I don't know how to proceed.

Does anyone have a solution for this issue?


texts:
For the last several weeks, I have been increasingly seeing this GET 500 error in the Chrome developer console of a published GEE app:
GET https://earthengine.googleapis.com/v1alpha/projects/earthengine-legacy/maps/[bigstring]/tiles/[2-dig-num]/[3-dig-num]/[3-dig-num] 500

Very often followed by:
Tile error: Earth Engine capacity exceeded.

The underlying operation often being performed when the error occurs is .mosaic()
Here is a screenshot:

Then, finally, an error often appears near the top of the map: "There was an error loading some parts of the map"
Followed by a fatal error in the middle of the map from which the only option is to close or refresh the browser: "Something went wrong while displaying this app. Please come back later."  Screenshot:

Of course I'd like to know what is causing the error to begin with.  But in the meantime, is there a way to trap / catch the error?  (I have tried putting try...catch blocks around the section of code, but the catch doesn't get triggered)
Anybody else having problems like this lately with GEE?


texts:
A few days ago I stumbled over a map which had beautifully rendered, detailed topographic tiles of the whole world. It was similar to google-physical with a strong contrast (shadow/light) and contour-lines. I do believe it was not from a commercial provider like ArcGIS but an open platform (Not sure though). I stumbled over it in relation to a search for openlayers.
Unfortunately I missed to bookmark the site, I can't find the link in my history anymore. Does anyone know what I'm mumbling about or could post examples of topographic tiles which fit the description?

This question is very similar to what I'm looking for (And much better worded), but it's 2 years old and partly outdated


texts: Estimating river Channel Width using Python/ArcGIS/MATLAB/R?

texts:
I am trying to estimate the river channel width from digitized polygon or poly-lines shapefiles. To do so, I converted poly-line shape files to points.shp at regular intervals (at 1 cm). Later, I calculated distance between corresponding points distance from both bank. However, points distribution are not uniform throughout both bank. Therefore, this method will yield much error in continuous width data set.
I have calculated manually at regular interval of 100 m of 20 km river length using distance tool in ArcGIS but I need this in automation mode which can provide continuously width of the river channel.

I am open to other method. Can we draw perpendicular lines from right bank to left bank or vice-versa? Or can we draw transect lines at user defined interval across the river? If this possible then we can convert these perpendicular lines into shapefile. Later we can use intersect tool to get channel width.
I posted snap of the points from polylines below for which I need to estimate width.This is just small reach of the river.
Has anyone done anything like this before?
If someone have better idea I would be happy to use that.


texts: QGIS Orfeo toolbox segmentation error "\otbcli_Segmentation" is not recognized

texts:
I have an old aerial image that I am trying to use Orfeo Toolbox's Segmentation tool in QGIS on. The image is grey-scale. I don't know if this is a problem or not, but I am experimenting at this point, I haven't used image analysis before. It would be great if I could get this to work, I have a few images like this that I would like to use.

I found a youtube video to run a meanshift segmentation to get me started. The video left the defaults in the options window of Segmentation, so I did the same. The log of the operation follows:

Algorithm Segmentation (meanshift) starting...
\otbcli_Segmentation
-in
"C:\Users\d\Documents\GIS\Musklutifs\Musk002.tif"
-filter
meanshift
-filter.meanshift.spatialr
5
-filter.meanshift.ranger
15
-filter.meanshift.thres
0.1
-filter.meanshift.maxiter
100
-filter.meanshift.minsize
100
-mode
vector
-mode.vector.outmode
ulco
-mode.vector.neighbor
true
-mode.vector.stitch
true
-mode.vector.minsize
1
-mode.vector.simplify
0.1
-mode.vector.layername
layer
-mode.vector.fieldname
DN
-mode.vector.tilesize
1024
-mode.vector.startlabel
1
-mode.vector.out
"C:\Users\d\AppData\Local\Temp\processing5cb9f70856de45bfa3be2c41f8dad49d\04296f5ee3534a49a59470503afd57aa\mode.vector.out.shp"
'"\otbcli_Segmentation"' is not recognized as an internal or external command, operable program or batch file.
Converting outputs
Loading resulting layers

The following layers were not correctly generated.
Output vector file
You can check the log messages to find more information about the execution of the algorithm

To see if it was just the grey-scale image, I also used the tool with a color raster. The operation failed with the same log message. I just updated to QGIS 2.18.3 with the OTB-__ library (forgot which one it's supposed to be) in the advanced install before trying to use this tool.


texts: Shifting all features in vector dataset using Bash/OGR

texts:
Let's say that I put together a Shapefile and all the features have their vertices shifted by a constant amount.

What's the easiest way of shifting all the features (hence the (x,y) position of their vertices) by an arbitrary shift?

I have lots of files that I would apply this correction to, so a Bash/OGR answer would be preferred.

I ended up using Spatialite for this, as it has the nice function ShiftCoords.


texts:
I am using GDAL to convert a .tif to .shp file.

gdal_polygonize.py input.tif -f “ESRI Shapefile” out.shp

However this does not work and results in:

Traceback (most recent call last):
  File "/usr/local/bin/gdal_polygonize.py", line 166, in <module>
    drv = ogr.GetDriverByName(format)
  File "/Library/Frameworks/GDAL.framework/Versions/1.11/Python/2.7/site-packages/osgeo/ogr.py", line 5255, in GetDriverByName
    return _ogr.GetDriverByName(*args)
TypeError: in method 'GetDriverByName', argument 1 of type 'char const *

I don't believe there is anything wrong with the .tif file or the command


texts:
I have MOD13Q1 NDVI images. Data are named in this format:

MOD13Q1.A2007049.h26v06.006.2015161224938_MODIS_Grid_16DAY_250m_500m_VI.tif

(YYYYDDD) Here 2007 is a year and 049 is day in year.

I have to rename it in this format: (YYYY.MM.DD) 2007.02.18

I am renaming it manually for every single file. How to perform a batch process for renaming all the files?


texts: Create in memory polygon with RasterToPolygon

texts:
I need to create a number of polygons in a loop from a raster. Essentially, the loop will call a function a number of times which needs to convert it to a polygon:

def myfunc(raster):

    #do some raster processing
    return arcpy.RasterToPolygon_conversion(raster, "in_memory\poly")

But the object it returns seemingly has nothing in it? The end result is I want to create a single shapefile from all the polygon features I will create. It seems silly to call RasterToPolygon, writing out a temporary shapefile for each polygon and then collecting them altogether as they will easily all fit in the RAM available.

Is there any way to output a Polygon geometry from RasterToPolygon ?. I.e, can I not do some thing like this:

poly = arcpy.RasterToPolygon_conversion(raster, arcpy.Polygon)


texts: Accurately tracking location in outback?

texts:
I am looking for a way to efficiently track my whereabouts/tracks whilst traveling through the Australian bushland. In the past, I've just recorded my longitudes and latitudes from a basic GPS onto a notebook and plotted dots to roughly estimate my tracks however the accuracy is low and the process time-consuming.

Was there any GPS device or software I could use that could track/record my location and export that as a trail onto a map or as list of coordinates after the work's finished?


texts: Finding lines of similar attribute with respect to nodes connected to its both Ends

texts:
I have line networks with node generated for Start and End Point . how to check if a line has correct attribute with respect to point created at its End. for instance FuID_Line = ABC_BDQ and Start Point having FuID_Po =ABC and End point having FuID = BDQ.

Update:

I think I am not clear enough previously in my question. actually I dont want to generate points. I have lines and with end points(Start/End) already given what I want to check is if line has
FuID = 560_960 attribute. it should have start point value pID=560 and End point value pID=960.

I want to join attributes to line feature something similar to this one

+--------+---------+-----+-------+---------+
| obj_id |  Fuid   | pid | pid_1 | remark  |
+--------+---------+-----+-------+---------+
|      1 | 560_960 | 560 |   960 | match   |
|      2 | 960_560 | 960 |   560 | match   |
|      3 | 261_690 | 260 |   690 | unmatch |
|      4 | 420_921 | 420 |   926 | unmatch |
+--------+---------+-----+-------+---------+

Note:Lines are sometimes overlapped to each other.


texts: Axis of rendered shapefiles not showing numbers that are latitude and longitude

texts:
I am using rgdal to render some shapefile. The first is a set of building zones from NYC Open Data. The second is the NYC coast line from OSM.

The Zone files seems to generally render an X axis between 850000 and 1150000. And a Y axis with numbers from 150000 and 250000. The OSM coastline renders an X-axis between -8400000 and -8100000 (note the difference in order of magnitude as well as being negative). The OSMs Y axis is between 4900000 and 5050000.

I'm using this as my example and I don't the author doing anything special to format the axes.

My code is here:

# render a plot based of shapes in nycgiszoningfeatures_201311_shp/
# shapes downloaded from https://nycopendata.socrata.com/api/geospatial/kdig-pewd?method=export&format=Shapefile
require(rgeos)
require(maptools)
require(rgdal)

# OSM Coastline of NYC Data from http://osm-extracted-metros.s3.amazonaws.com/new-york.coastline.zip
mapBase <- readOGR('new-york-city.coastline', 'new-york')

shapeDir <- 'nycgiszoningfeatures_201311_shp'
shapeFiles <- dir(shapeDir, '\\.shp$')
shapeLayers <- sub('\\.shp$', '', shapeFiles)
# to see the files sprintf("%s/%s", shapeDir, shapeFiles)
# TODO: I'm told for loops are evil in R. Find a better way
for (i in 1:length(shapeLayers) ) {
    map <- readOGR(shapeDir, shapeLayers[i])
    #plot(mapBase, col="khaki2",border="grey", axes=TRUE, xlim=c(-8300000,-8200000), ylim=c(4900000, 5050000))
    #plot(map, add=TRUE, col="gray",border="blue")
    plot(map, col="gray",border="blue", axes=TRUE)
}

How do I get these shapes to plot with longitude and latitude numbers for the axes instead of those numbers.


texts: Creating simple popup on click in OpenLayers

texts:
I have different figures created in my Django project.(icon with text, point,circle etc)
but I try to create a simple popup with the coordinates when I click the mouse but I have a error:

My code is from an example, seen here

var map = new ol.Map({
    target: 'map',
    layers: [
        new ol.layer.Tile({
            source: new ol.source.OSM()
        })
    ],
    view: new ol.View({
        center: ol.proj.transform([-0.92, 52.96], 'EPSG:4326', 'EPSG:3857'),
        zoom: 6
    })
});

var popup = new ol.Overlay.Popup();
map.addOverlay(popup);

map.on('click', function(evt) {
    var prettyCoord = ol.coordinate.toStringHDMS(ol.proj.transform(evt.coordinate, 'EPSG:3857', 'EPSG:4326'), 2);
    popup.show(evt.coordinate, '

Coordinates

' + prettyCoord + '

');
});

I get this error:

  Uncaught TypeError: ol.Overlay.Popup is not a constructor.

when trying to create --> var popup = new ol.Overlay.Popup();


texts: GeoServer 2.18.1 can't add Data Source to Oracle Database 10g Enterprise Edition Release 10.2.0.1.0

texts:
Target: Adding Data Source to Oracle 10g Enterprise Edition Release 10.2.0.1.0
Problem: When trying to add Data Source - Oracle NG GeoServer show the error message:

Error creating data store, check the parameters.

Error message:

Unable to obtain connection: Cannot create PoolableConnectionFactory
(ORA-00604: error occurred at recursive SQL level 1 ORA-12705: Cannot
access NLS data files or invalid environment specified )

Server configuration:
Apache Tomcat/9.0.13
Geoserver (war) 2.18.1.
Drivers Oracle placed in:
D:\Tomcat9\webapps\geoserver_new\WEB-INF\lib\gt-jdbc-oracle-24.1.jar
D:\Tomcat9\webapps\geoserver_new\WEB-INF\lib\ojdbc8-19.7.0.0.jar
LOG:
Caused by: java.sql.SQLException: ORA-00604: error occurred at recursive SQL level 1 ORA-12705: Cannot access NLS data files or invalid environment specified

  at oracle.jdbc.driver.T4CTTIoer11.processError(T4CTTIoer11.java:509)
  at oracle.jdbc.driver.T4CTTIoer11.processError(T4CTTIoer11.java:456)
  at oracle.jdbc.driver.T4CTTIoer11.processError(T4CTTIoer11.java:451)
  at oracle.jdbc.driver.T4CTTIfun.processError(T4CTTIfun.java:1123)
  at oracle.jdbc.driver.T4CTTIoauthenticate.processError(T4CTTIoauthenticate.java:552)
  at oracle.jdbc.driver.T4CTTIfun.receive(T4CTTIfun.java:553)
  at oracle.jdbc.driver.T4CTTIfun.doRPC(T4CTTIfun.java:269)
  at oracle.jdbc.driver.T4CTTIoauthenticate.doOAUTH(T4CTTIoauthenticate.java:501)
  at oracle.jdbc.driver.T4CTTIoauthenticate.doOAUTH(T4CTTIoauthenticate.java:1292)
  at oracle.jdbc.driver.T4CTTIoauthenticate.doOAUTH(T4CTTIoauthenticate.java:1025)
  at oracle.jdbc.driver.T4CConnection.logon(T4CConnection.java:747)
  at oracle.jdbc.driver.PhysicalConnection.connect(PhysicalConnection.java:793)
  at oracle.jdbc.driver.T4CDriverExtension.getConnection(T4CDriverExtension.java:57)
  at oracle.jdbc.driver.OracleDriver.connect(OracleDriver.java:747)
  at oracle.jdbc.driver.OracleDriver.connect(OracleDriver.java:562)
  at org.apache.commons.dbcp.DriverConnectionFactory.createConnection(DriverConnectionFactory.java:38)


texts:
I'm a GIS student and in my project I need to generate vegetation indexes. I've generated it succesfully yesterday using both python2.7 and grass gui .I'm using WinGRASS-6.5.svn-r54428-453 on Win7. Today -I don't know why- grass started to generate 1 pixel NDVI map. The formula is as simple as it could be: float(b4-b3)/(b4+b3) as most of you already know. My data is not corrupted, I'd saved original maps somewhere else and it was my first check that whether there is a corruption. Also using QGis 1.8.1 I can generate correct NDVI map with same data. You can check it here:

r.info map=ndvi4@PERMANENT
 +----------------------------------------------------------------------------+
 | Layer:    ndvi4@PERMANENT                Date: Tue Jan 01 04:06:35 2013    |
 | Mapset:   PERMANENT                      Login of Creator: murat           |
 | Location: newLocation                                                      |
 | DataBase: D:\gis\grassWS                                                   |
 | Title:     ( ndvi4 )                                                       |
 | Timestamp: none                                                            |
 |----------------------------------------------------------------------------|
 |                                                                            |
 |   Type of Map:  raster               Number of Categories: 255             |
 |   Data Type:    FCELL                                                      |
 |   Rows:         1                                                          |
 |   Columns:      1                                                          |
 |   Total Cells:  1                                                          |
 |        Projection: UTM (zone 36)                                           |
 |            N:          1    S:          0   Res:     1                     |
 |            E:          1    W:          0   Res:     1                     |
 |   Range of data:    min = nan  max = nan                                   |
 |                                                                            |
 |   Data Description:                                                        |
 |    generated by r.mapcalc                                                  |
 |                                                                            |
 |   Comments:                                                                |
 |    float(L71175034_03420050627_B40@PERMANENT -                             |
 |    L71175034_03420050627_B30@PERMANENT) /                                  |
 |    (L71175034_03420050627_B40@PERMANENT +                                  |
 |    L71175034_03420050627_B30@PERMANENT)                                    |
 |                                                                            |
 +----------------------------------------------------------------------------+
(Tue Jan 01 04:06:51 2013) Command finished (0 sec)

I uninstalled WinGRASS-6.5.svn-r54428-453  and installed WinGRASS-7.0.svn-r54453-453, I've deleted mapset and location and the recrated them also. I couldn't think anything else please give me some information.


texts: formatting issues with an attribute table in the QGIS Atlas

texts: I'm having issues formatting an attribute table in the Atlas of QGIS. All the data are linked correctly, but the values and words are not centered in the middle of the rows (see picture). Changing the cell margins and lie width does not fix the problem. Under appearance, Display Header = 'on first frame', Empty tables = 'draw headers only', oversized text = 'wrap text'. Under Frames, I have "use existing frames" selected. And I have tried different combinations on the appearance, which don't seem to help. Does anyone know of a solution?


texts: Importing ERDAS-processed GTOP030 dem data into ArcGIS Desktop?

texts:
I have two GTOP030 DEM images. After downloading them, I have imported into ERDAS and mosaiced them. Then I have exported it to .tif file in order to open it in ARCGIS. But when I tried to open it in ARCGIS, "the memory could not be 'read'" error message appeared on the screan and turned off ARCGIS. what's wrong with this.

Do you have any idea for this?

Another thing I ask from you is how to remove No data value from this DEM file because min value of this DEM image is -9999?


texts: Extracting data from Oracle Spatial and preserving spatial data

texts:
My table currently contains an SDO column that I want to export (along with the rest of the data) but all I have at my disposal is PL SQL Developer.

Is there anything I can do to get at that data and put it into a text or xml format?

I'm not accustomed to working with Oracle Spatial data.


texts: Getting AttributeError when exporting Heatmap using qgis2web?

texts:
I want to make a heatmap layer on my webgis.

using these options

but when i'm trying to export it using qgis2web i always get this error.

2019-06-21T10:15:53     WARNING    Traceback (most recent call last):
              File "C:/Users/Rizki Permana/AppData/Roaming/QGIS/QGIS3\profiles\default/python/plugins\qgis2web\maindialog.py", line 323, in previewMap
              preview_file = self.createPreview()
              File "C:/Users/Rizki Permana/AppData/Roaming/QGIS/QGIS3\profiles\default/python/plugins\qgis2web\maindialog.py", line 286, in createPreview
              dest_folder=utils.tempFolder()).index_file
              File "C:/Users/Rizki Permana/AppData/Roaming/QGIS/QGIS3\profiles\default/python/plugins\qgis2web\leafletWriter.py", line 100, in write
              folder=dest_folder)
              File "C:/Users/Rizki Permana/AppData/Roaming/QGIS/QGIS3\profiles\default/python/plugins\qgis2web\leafletWriter.py", line 268, in writeLeaflet
              useVT, useShapes, useOSMB)
              File "C:/Users/Rizki Permana/AppData/Roaming/QGIS/QGIS3\profiles\default/python/plugins\qgis2web\leafletLayerScripts.py", line 88, in writeVectorLayer
              new_obj = heatmapLayer(layer, safeLayerName, renderer)
              File "C:/Users/Rizki Permana/AppData/Roaming/QGIS/QGIS3\profiles\default/python/plugins\qgis2web\leafletLayerScripts.py", line 430, in heatmapLayer
              hmWeightId = layer.fieldNameIndex(hmWeight)
             AttributeError: 'QgsVectorLayer' object has no attribute 'fieldNameIndex'

How do I solve it?


texts:
I need to visualize in a map the best month to do something in an area.

I've thought of classifying months in 4 groups: (Dec, Jan, Feb // Mar, Apr, May // Jun, Jul, Aug // Sept, Oct, Nov) and displaying them by colors according the season temperature: Blue - cold, red - warm.

The problem that I find with this approach is that it is a global map, and does not make much sense to show January and February with the same color in Australia and Italy.

Any Idea on how I could display global month/season data in a map having in count this?


texts: ArcGIS spatial join - WITHIN A DISTANCE including distance

texts:
I am applying a spatial join in ArcGIS between points (target features) and lines (join features). I want to have a 1-to-many join, so I use the match option 'within a distance'. So far so good.
The annoying thing is that I need to have the distance between the point and the matched lines within the search radius. I do not find any solution for this. Is there someone who has a good tip for me that will do the trick?

The final goal is that I need to join the points to the line that is nearby, but also has some other characteristics (equal value on some attributes, but also a timestamp that needs to be higher in 1 class than the other which I don't succeed in working with in arcGIS anyway). These kind of data preparation I do afterwards in other tooling. That's why I can't just take the one that is closest...


texts: Grouping gps-points by location and time using QGIS?

texts:
I tracked my holiday-routes with some friends for many years. We tracked 3 cars each year, every car with another GPS-Tracker. The GPS-Tracker locate the position every 10 minutes. I'm using QGIS 3.2 and have at least about 1 million GPS-Points.

Now i want to route every year as a line, not the single cars.
I have Attributes about:
1. Identifier per way
2. timestamp
3. lat
4. lon

Would it work with a time-distance-matrix to group the points by location and attributes?


texts: How to change the port on which GeoServer 2.13 runs on when it is installed as platform independent binary on Linux

texts:
Note this question IS NOT about how to change Tomcat port or jetty port or Mac or Windows native installation of GeoServer. Each configuration for different platform is different. I tried many different ways to change port to 8585 from default 8080 but it ignores my changes and still gets launched on 8080. I tried changing startup.sh and added -Djetty.port=8585 in the exec part and tried jetty-http.xml but nothing works.
I am not surprised this is not documented anywhere on their manual.


texts: Seeking border of Netherlands to use in ArcGIS Pro?

texts:
I am a beginner in ArcGIS Pro. I asked a question before (Expanding IDW interpolation to borders of country in ArcGIS Pro?), but I feel like I do not understand the basics enough, so a more basic question now.

I have a dataset with concentrations measured at 20 locations in the Netherlands. I want to cover all of the Netherlands with IDW interpolation/extrapolation, using the Geostatistical Wizard.

So not like this, but everything up until the borders:

I understand I have to set Analysis/Environments/Processing extent to the border of the Netherlands.

How can I make a shapefile of the border of the Netherlands so that I can use it to define the processing extent?

I used this, but that doesn't work:

https://learn.arcgis.com/en/projects/cartographic-creations-in-arcgis-pro/lessons/add-and-symbolize-countries.htm


texts:
I have an odd shaped (not square) raster .tiff file that I converted to MrSid using Lizardtech's GeoExpress. It converts the nodata to black. This is applying a value to no data.

Eventually I want to overlay this MrSid image on top of other layers, so the nodata cells need to be transparent.

How do I do this?


texts: Need Help on Multiple External SLDs for GeoServer LayerGroup

texts:
I have a layer group called "basemap" in my GeoServer made up of 3 layers but their SLDs needs to be styled dynamically so I have decided to create a XML web service for the external SLDs but I am wondering whether I can pass multiple SLD files as shown below.

How do I go about it?

Please check my SLD parameter in my code below if that is how to do it.

Or should I stop using layer group and call the layers one by one and pass its external SLD link to it?

My Codes Below In OpenLayers

var getBasemap = function(){

    GlobalBasemapLayer = new OpenLayers.Layer.WMS(
        "Basemap",
        geoserverWMSurl,
        {   LAYERS: 'workspace:basemap',
            format: format,
            tiled: true,
            transparent:true,
            sld:'http://domain.com/getSLD/Layer1SLD,http://domain.com/getSLD/Layer2SLD,http://domain.com/getSLD/Layer3SLD'
        },
        {
            buffer: 0,
            displayOutsideMaxExtent: true,
            isBaseLayer: false
        }
    );
};


texts: ArcGIS Javascript API with Windows Authentication

texts:
I am using the ArcGIS Javascript API to consume a feature layer on an ArcGIS Server secured with Windows Authentication.  I want to satisfy the authentication programmatically so the web user isn't prompted with a login.

I've poked around the web and seen something about downloading an ASP.NET proxy page, but it is from the 9.3 section of the online ESRI resources.  I'm just wondering if there is a different way to do it in 10.x versions?


texts: Setting Graphics rendering level like ArcMap UniqueValueRenderer in ArcGIS AP for JavaScript

texts:
I have to add a Map Layer on Map which I am able render with UniqueValueRenderer. After rendering I want to have certain value matching graphics always on Top if they are over lapping. But I am not able to set Like we set symbol levels  in ArcMap.

    let green = new DirectionalLineSymbol({ style: SimpleLineSymbol.STYLE_SOLID, color: new Color([50, 200, 20,1]), width: 5, directionSymbol: 'arrowf', directionSize: 20, directionPixelBuffer: 50, showStartSymbol: false });
              let yellowf = new DirectionalLineSymbol({ style: SimpleLineSymbol.STYLE_SOLID, color: new Color([255, 210, 2, 1]), width: 5, directionSymbol: 'arrowf', directionSize: 20, directionPixelBuffer: 50, showStartSymbol: false });
              let yellowr = new DirectionalLineSymbol({ style: SimpleLineSymbol.STYLE_SOLID, color: new Color([255, 210, 2, 1]), width: 5, directionSymbol: 'arrowf', directionSize: 20, directionPixelBuffer: 50, showStartSymbol: false });
              let red = new DirectionalLineSymbol({ style: SimpleLineSymbol.STYLE_SOLID, color: new Color([255, 0, 0, 0.7]), width: 5, directionSymbol: 'arrowf', directionSize: 20, directionPixelBuffer: 50, showStartSymbol: false });
              let renderer = new UniqueValueRenderer(green, function (graphic) {
                  if (segments.find(e => e.objectid == graphic.attributes.OBJECTID && !e._deleted && !e.is_reverse)) {
                      graphic.setSymbol(yellowf);
                      return 'YELLOW_F'
                  }
                  else if (segments.find(e => e.objectid == graphic.attributes.OBJECTID && !e._deleted && e.is_reverse)) {
                      graphic.setSymbol(yellowr);
                      return 'YELLOW_R'
                  }
                  else if (segments.find(e => e._deleted == true && graphic.attributes.OBJECTID == e.objectid)) {
                      graphic.setSymbol(red);
                      return 'RED'
                  }
                  else if (graphic.attributes.MAINT_RD_TRIM && graphic.attributes.RDWAY_ID) {
                      graphic.setSymbol(green);
                      return 'GREEN'
                  }
                  else {
                      graphic.setSymbol(red);
                      return 'RED'
                  }
              })
              //get data fro m layer manually and sort and add to the map to fix road  on top
              renderer.addValue({ value: 'GREEN', symbol: green, label: 'Valid Segments' ,level :3})
              renderer.addValue({ value: 'RED', symbol: red, label: 'Invalid Segments' ,level :2});
              renderer.addValue({ value: 'YELLOW_F', symbol: yellowf, label: 'Selected Segments',level :1  })
              renderer.addValue({ value: 'YELLOW_R', symbol: yellowf, label: 'Selected Segments' ,level :0})//changed from yellowr
              layer.setRenderer(renderer);

though above code snippet I am able to render graphic in red, yellow green but I want Yellow color graphics always on TOP which is not happening. I did not find any thing in API Documentation too.


texts:
I have a shapefile (either polygon or polyline) that I need to convert to an "Area Data" .dat file to use with another program (it's already older, it says it wants "ArcInfo Format").

I've tried Googling and playing around with the different export functions, but so far I haven't had any luck.

I'm fairly new to working with ArcGIS. I have ArcGIS Desktop 10.5 and the other program I need to load the map to is a very specialized program for theodolite tracking, Pythagoras (http://www.cetaecoresearch.com/research-software-pythagoras.html).


texts: Common class range for different attributes in a same polygon shapefile

texts:
I have a polygon shapefile which has population attributes pop90,pop91 and so on for 10 years. They have different minimum and maximum value for each of the years. Now I need to do classification for each attribute(pop90,pop91) but with the same number of classes and the range for each class should be same. Same number of classes is fine but I am not sure on how can I use the same class range for each of the attribute. Any idea how can this be done.
Thank You.


texts: Seeking plugin like NetCDF Browser for QGIS 3.0 and above

texts: Is there a plugin that fulfills the same function as NetCDF Browser (i.e. to browse multi-variable and multi-dimensional NetCDF files) but is available for QGIS versions 3.0 and above?


texts:
How do you use ogrinfo with a SQL Server database?

I've tracked down documentation on it at GDAL and it mentions ...

  datasource_name:
  The data source to open. May be a filename, directory or other virtual name. See the OGR Vector Formats list for supported datasources.

and I'd hope this has an answer, but the link OGR Vector Formats gives a 404 error


texts: Identify tool and highlighting behaviour in QGIS 2.6

texts:
I have a problem with the identify tool in QGIS 2.6.

Somehow the object highlighting doesn't go off when i click on another object with the identify tool in QGIS.
In earlier versions, when you clicked on object after object, only the last object clicked was highlighted, highlighting on the previous one was turned off.
Note- it only happens when i let QGIS auto-open feature attributes form and highlight of an object stays on until i close the form.
Is there a hidden feature somewhere, that i can use to change this behaviour?


texts:
I want to draw a simple graphics rectangle\ circule in QGIS 2.12, just as i do it in ArcGIS:

I red

Digitizing rectangles/ovals in QGIS?
Drawing tools to add circles, arrows, lines for QGIS?

and installed the tool "Rectangles ovals digitizing in QGIS", but it works only when i edit a layer. I want graphics rectangle box and not editing a layer, just has in the Print Composer Layout


texts: Moving individual labels & rule-based labelling in QGIS

texts:
I have been playing around with labels and came to a problem. Hopefully someone can pinpoint whether I'm doing something wrong or there is a systemic hurdle in Qgis.

Are you allowed to move labels manually in Qgis IF the chosen type of labelling is "rule-based"?

When I change labelling type to Single labels, moving individual labels works like a dream. But under rule-based it just doesn't work, no matter what I do.

I'm using the Bonn version (Qgis 3.2.3), and have followed the BurdGis Youtube tutorial, as well as some other forum threads.


texts: How to build a query/expression to select specific points within a buffer? (in QGIS)

texts:
I have a point layer in QGIS. Each point represents a tree and each tree has information about its height as an attribute.
I would like to build a query/expression to select all those trees that have an other tree higher than 7 meters within their 20-meter-buffer (20 meter radius).

I managed to do this "manually" by creating a 20-meter-buffer around all trees higher than 7 meters (after the selection of the trees higher than 7 meters) and then selected the ones within the buffer but I was wondering if its possible with a query/expression.
(maybe with a bit more complex one)

Unfortunately I do not have PostGIS so this should be done with Select by expression or Query Builder if possible.

(It would be also helpful to know if it's simply impossible)


texts: Remove Leading Zeros with QGIS Field Calculator

texts:
I'd like to take this attribute value: "MAP_PAR_ID": "0160011000", and populate a new field with "16-11".

This works: concat(substr("MAP_PAR_ID",0,3), '-', substr("MAP_PAR_ID",5,3)) but i need to detect leading zeros and not include them

I've tried numerous variations of CASE ELSE statements but can't even detect leading zeros. Feeble attempt below:

CASE
    WHEN substr("MAP_PAR_ID", 0, 1) = '0' THEN substr("MAP_PAR_ID", 1, 2)
    ELSE  substr("MAP_PAR_ID", 0, 3)

    WHEN substr("MAP_PAR_ID", 5, 1) = '0' THEN substr("MAP_PAR_ID", 6, 2)
    ELSE substr("MAP_PAR_ID",5, 3)
END

Would the python console make things a little easier?


texts: Google Earth app KMZ overlay hidden by 3D tees and buildings

texts:
I've got a KMZ layer with an overlay layer in the <groundoverlay> tag.
I've got 3d buildings and trees in the overlay area, and they hide the overlay.

In Google Earth app it is not possible to disable 3d buildings.
Do you know if there is a tip to have the overlay always visible and attached to the ground in Google Earth app for smartphones and tablets?


texts:
I have been working with the 'lidR' package to analyze LiDAR tiles and was able to successfully complete 14 analyses on 12/1. When I started running the script on 12/2 I kept getting this error on all of my tiles, even tiles I was able to analyze the day before without issue (proj_crs_get_datum: Object is not a SingleCRS):
tile <- readLAS("w0929n1781.las")
#> proj_crs_get_datum: Object is not a singleCRS
#> proj_crs_get_datum: Object is not a singleCRS

I have tried deleting and reinstalling lidR, redowloading tiles and running the script again, updating my version of R from 3.6.1 to 3.6.3, looking for any and all suggestions.


texts:
Please look at the color legend of the image:

This image was plotted in ferret program, which has different breaks, like

0-3 with 0.5 interval,
3-9 with 3 interval,
9-10 with 1 interval,
10-20 with 2 interval and
20-70 with 50 interval.

The legend scale was uniformly graduated but labelled as per breaks given. Can we get such legends in R 'Raster - plot'? or any other packages in R?

I have tried in R (raster) and got this image:

In this, the legend is not looking good. I need the a legend like in first image.

The R codes I used are following,

library(raster)
data=raster("dfile.nc") # MODIS-AQUA satellite data
breakpoints=c(0:0.2, 0.21:2,2.01:5,15)
colpatt=rev(colorRampPalette(c("red", "orange", "blue"))(6))
plot(data,breaks=breakpoints, col=colpatt)


texts: Adding points to represent non-chosen locations?

texts:
I am looking for a way to set up a discrete choice model in ArcGIS for foraging behavior of individuals, in a grid of 15 locations. To do that I need a point for every possible choice they could have made, currently I only have data for the choice they did make. Right now my data looks a little like this:

IndivID    Date    LocationID   X    Y
    1      1/1         1       ...  ...
    1      1/2         4       ...  ...
    1      1/3         3       ...  ...
    1      1/4         2       ...  ...
    1      1/5         3       ...  ...
    1      1/6         1       ...  ...
    1      1/7         2       ...  ...

And what I am looking for would be like this (assuming I have 4 locations, and 0 is the choice to not forage for that day).

IndivID    Date    LocationID   X    Y    Choice
    1      1/1         0       ...  ...      0
    1      1/1         1       ...  ...      1
    1      1/1         2       ...  ...      0
    1      1/1         3       ...  ...      0
    1      1/1         4       ...  ...      0
    1      1/2         0       ...  ...      0
    1      1/2         1       ...  ...      0
    1      1/2         2       ...  ...      0
    1      1/2         3       ...  ...      0
    1      1/2         4       ...  ...      1

I know I would add the Choice field before adding new points and calculate it as 1, so as I add new points I can see what the original choice was, and I have X,Y coordinates for the centroids of each location and a home location I can use to spatially locate the new points, but I am not sure how to go about adding in the extra points to represent the non-chosen options. To give you an idea of the size of my dataset, I have 57 individuals and I am doing this analysis for an entire year (365 days), and actually have 15 sites (used 4 as an example).

Anyone have any ideas how I would go about setting this up, or do you think this type of process is out of the range of ArcGIS's capabilities?

I am specifically looking for an alternative to toiling away at an excel or csv file to make the points, so it can be tweaked and reapplied if I refine the location choices to more smaller plots.