texts: Limiting number of selected items in CARTO category widget?

texts:
I created a map with several widgets with CARTO but I'd like to limit the number of selected items in a category widget to one. By default you can select all of the items and in my use case that is not desirable, because there's some points overlapping making it confusing for the users.

I've already looked through the documentation and did some tests, but couldn't find a way to do it, so I'm asking for help here.


texts: Earth Engine sampleRegions: keep track of original pixels geometry/location?

texts:
I am trying to obtain (for further export), a dataset grouping my NDVI values by fields, for which I have polygons.

To do so, I use the Google earth engine function sampleRegions: this takes an image, a Feature Collection, and allows to merge the two,  giving a feature collection where each features contains the values of image's pixels contained in each polygon, associated with the values of the corresponding polygons.

The problem is that by returning a feature collection, it loses track of the geometry/location information of all pixels. That is, I can see the values of all pixels in a given field, but do no know to which location within the field these pixels belong.

Any idea on how I could retrieve these?

Or should I adopt another procedure?

// Delineate fields
var Field1 = ee.Geometry.Polygon([
  [[-122.37158060073853,38.458125582533334],
    [-122.3727822303772, 38.45752068253876],
  [-122.37192392349243,38.45646209534339],
  [-122.37063646316528,38.457033398298925],
  [-122.37158060073853,38.458125582533334]]
  ]);
var Field1_F  = ee.Feature(Field1,{id: 1});

var Field2= ee.Geometry.Polygon([
    [[-122.37115144729614,38.45536988592448],
    [-122.369863986969, 38.4560084103617],
    [-122.3706579208374,38.456966186421305],
    [-122.37196683883667,38.45636127670514],
    [-122.37115144729614,38.45536988592448]]])

var Field2_F  = ee.Feature(Field2,{id: 2});
var Fields_FC = ee.FeatureCollection([Field1_F,Field2_F]);

// get NDVI
var NDVI = ee.ImageCollection('LANDSAT/LE7_L1T_32DAY_NDVI').
mean();

// Merge NDVI and field
var pixelsID = NDVI.sampleRegions(Fields_FC, ['id'], 30);

Map.setCenter(-122.3713231086731,38.456663732197235, 17)
Map.addLayer(Fields_FC, {}, 'feature');
Map.addLayer(NDVI,
    {min:-1, max:1, opacity: 0.9, palette: ['blue', 'green', 'red']},         "NDVI")

// this is almost the perfect output, yet does not contain pixel location
print(pixelsID, "pixelsID")


texts: Cropping EUDEM with GDAL creates zig-zag terrain

texts:
I have downloaded EUDEM (20 Gb) from http://www.eea.europa.eu/data-and-maps/data/eu-dem#tab-original-data

When I try to crop it by GDAL using:

gdalwarp  -dstalpha -of EHdr -crop_to_cutline -cutline crop.shp eudem_dem_4258_europe.tif  eudem_cropped.flt

I get bad terrain. The terrain gets zig-zag in heights by 250 meters. Please check the screenshots:

What do you think is the problem?


texts: Export Address Inspector Address Points into its own layer?

texts:
I am a new ArcMap user.

I have the ArcMap Base Map, some county parcel data layers and found two important addresses from the data using the 'Address Inspector' as shown in the pictures and placed a point.

Is there a way to turn these points into its own separate layer? In the layer pallet, all I have is the Base Map, and Parcel Data layers.


texts: Is there a way to ignore the "handle bad layers" error message without removing the layers from the project?

texts:
I have multiple QGIS projects that use the same base layers. I keep the QGIS projects on dropbox so I can access them from different computers. The base layers are scattered all over the place, on multiple servers, or buried deep in the files of one computer or another. Basically, the file organization here is a mess, but fixing it is a problem above my pay grade. The problem is that sometimes I want to make changes to a project when I'm not connected to one of the servers. I can get along without these layers temporarily, but I want them to be available when I reconnect to that server. I have copies of some of the base layers on dropbox with the project, but storage space is limited so this is only practical for the layers I use most frequently.

Is there a way to ignore the "handle bad layers" error message that comes up when I open a project without removing the layers from the project?


texts: Compute vector layer area with python in qgis

texts:
I need to identify in QGIS the area of a Qgsvectorlayer .

Anyone have an idea how to do that in python ?


texts: "Delete this Add-In" event for uninitializing own ArcMap Add-In

texts:
Is there an event that gets triggered, when removing an Add-In by pressing the "Delete this Add-In" button in the Add-In Manager - using ArcMap 10.x?

I want to listen to this event to let my Add-In automatically remove any toolbars, close all dockable windows that came with my Add-In, etc.


texts: Import SHPs and DBFs to PostgreSQL in pgAdmin 4 on a Mac

texts:
I am receiving the following error when trying to import a dbf file into PostgreSQL (v 12, pgAdmin 4) using the Boundless/OpenGeo Suite for Mac pgShapeLoader GUI (which was recommended for Mac users in a Stack Exchange answer here, since there is no shapefile importer GUI included in the most recent Mac installation: Shapefile loader in pgAdmin 4):

Failed SQL begins: "CREATE INDEX ON "public"."soildef" USING GIST ("geom");
COMMIT;
"
Failed in pgui_exec(): ERROR:  column "geom" does not exist

Shapefile import failed.

The pgShapeLoader GUI worked fine for me when importing shapefiles, but it is producing the above error for the dbfs. And I know it should work because it is working for everyone using Windows for the PostGIS tutorial I'm following on Udemy (meaning, there is nothing wrong with the dbf itself).

The Windows installation of PostgreSQL/PostGIS has a shapefile loader gui shp2pgsql-gui.exe, but unfortunately this does not appear to be included in the Mac install. There is only a Unix executable shp2pgsql in the bin folder.

Are my only options for importing the dbfs to use the command line and/or to do this through QGIS, per another answer on the linked SE question/answer above? I am not finding any other solutions. Or maybe the question I should really be asking myself is simply: "How do you import a dbf into PostgreSQL on a Mac?" (Which seems a bit obvious now that I'm typing it, but I'd like to post this question anyway for anyone else who might run into this problem).


texts: How to store raster data outside postgis database

texts:
I am able to import raster data easily into the database and view it through QGIS. However,
I am unable to view the raster layer in Arcgis.
Now I want to know how to create a table with a reference to raster files outside the datbaase. However, there isn't much discuss or example with regards to how it can be done.
I want to see its benefits given that raster layers can be shared across various applications.

Please share your knowledge on this?

Thanks in advance.


texts: How to show layer attributes when mouse hovers on it with pyqgis?

texts:
My question is very clear. I want to show point features attributes when the mouse tool hovers on the feature. It can be like that:

If the mouse goes elsewhere the label should gone as well. And of course I know that the layer must be active layer for this.


texts: Install/config of QGIS-Server 3.16 on Fedora 33 & Rawhide (FC34)

texts:
Has anyone successfully installed QGIS-Server 3.16 on Fedora 33 or Rawhide (FC34)?
Thus far I have not had any success. Following the QGIS guide (which was written for Debian), the first action is to check that the install by running the executable, /usr/libexec/qgis/qgis_mapserv.fcgi, which is the default path from the Fedora install.
According to everything I have seen it should give a result as shown below:
      QFSFileEngine::open: No file name specified
Warning 1: Unable to find driver ECW to unload from GDAL_SKIP environment variable.
Warning 1: Unable to find driver ECW to unload from GDAL_SKIP environment variable.
Warning 1: Unable to find driver JP2ECW to unload from GDAL_SKIP environment variable.
Warning 1: Unable to find driver ECW to unload from GDAL_SKIP environment variable.
Warning 1: Unable to find driver JP2ECW to unload from GDAL_SKIP environment variable.
Content-Length: 206
Content-Type: text/xml; charset=utf-8

<ServiceExceptionReport version="1.3.0" xmlns="https://www.opengis.net/ogc">
 <ServiceException code="Service configuration error">Service unknown or unsupported</ServiceException>
</ServiceExceptionReport>

however I get the following:
[root@servnet /]# /usr/libexec/qgis/qgis_mapserv.fcgi
Application path not initialized
Application path not initialized
    - 'VirtualXPath'    [XML Path Language - XPath]
    - 'VirtualXPath'    [XML Path Language - XPath]
    - 'VirtualXPath'    [XML Path Language - XPath]
    - 'VirtualXPath'    [XML Path Language - XPath]
    - 'VirtualXPath'    [XML Path Language - XPath]
QFSFileEngine::open: No file name specified
    - 'VirtualXPath'    [XML Path Language - XPath]
Content-Length: 102
Content-Type: application/json
Server:  QGIS FCGI server - QGIS version 3.16.1-Hannover
Status:  400

[{"code":"Bad request error","description":"Requested URI does not match any registered API handler"}]  - 'VirtualXPath'    [XML Path Language - XPath]
    - 'VirtualXPath'    [XML Path Language - XPath]
    - 'VirtualXPath'    [XML Path Language - XPath]
Segmentation fault (core dumped)
[root@servnet /]#

The current configuration on the server is Fedora Rawhide (FC34), headless install, with xvbf installed to prevent issues with GET PRINT, but cannot even get the basics to work yet.
I am sure it is something simple like a missing package or library, but there is nothing in any log-files to give a clue.
Thus far I have checked the various questions here, and also various google searches on different keywords in the error message, but thus far, nothing seems to come close to a match
Any ideas or tips?


texts: Event change:visible not fired when map reaches layer minResolution or maxResolution

texts:
I'm listening to a layer's event change:visible. However it does not get fired when the resolution changes.

I.E. my layer as a minimum resolution of 0.5964957884324376. If I zoom past it, even though my layer is no longer displayed on the map, the event is not fired and its visible state is still true.

How can I detect my layer is no longer displayed because of its min/max resolution setup ?


texts: Is it possible to access ArcGIS's table operations "Summarize..." and "Statistics..." in python?

texts: ArcGIS has a few nifty table tools, including Summarize and Statistics. Is there a way to access them through arcpy, or will I have to hack my own methods via math and lists?


texts: Google Maps API: how to get the street name in multiple languages

texts:
I'm trying to get the street name in multiple languages from the Google Maps API V3. For example in Belgium there are multilingual street names:

If i set the language=FR i will get "Rue Brogniez" and language=NL results "Brogniezstraat".

How can i get these names together in one geocoder.geocode call?


texts: PolygonToCenterline: AttributeError: 'module' object has no attribute 'topographic'

texts:
I want to create a centerline from a river polygon using PolygonToCenterline. It works when I use the tool in ArcGIS, but when I run it from the script I get the following error:

AttributeError: 'module' object has no attribute 'topographic'

I need to be able to run it from the script. I use ArcGIS Pro 2.2.0 with an advanced license. I think I did it exactly as in the description. Any ideas what I could try? I added the code below:

# Import arcpy module
import arcpy
# Check out Production Mapping License
arcpy.CheckOutExtension ('foundation')
# Define variables
river_polygon = r'C:\zhaw\masterarbeit\data\diego\typical_rivers.gdb\birs_2011'
out_gdb = r'C:\zhaw\masterarbeit\data\diego\typical_rivers.gdb\centerline'
# Execute Polygon To Centerline
arcpy.topographic.PolygonToCenterline(river_polygon, out_folder)


texts: Ordnance Survey Tile Tracking Solutions?

texts:
We want to find a way to track how and how many tiles are drawn in the Ordnance Survey API (or through a different way if that is easier). As you may know, the royalty charges are per tile so we want to see if there is a way to track how tiles are drawn in our app without writing a software solution. It may overall end up that we'll have to write one (considering how they want to be invoiced quarterly), but for budgeting and testing purposes if there is a quick / rudimentary way to do it before we get down to finer details, that is the solution I am looking for. So far I've been able to use the inspect tool in Chrome to view how many tiles get loaded in one view of our app but a scaleable solution is what I am looking for.

For reference, each tile is 256x256 regardless of scale (same as OSM so maybe a solution people use for OSM might work in this case). I was thinking of just doing some good ol' math to figure out how many tiles might fit within a predetermined boundary we could set up in our app, however that isn't very precise given that I don't know how many tiles might be loaded outside of the boundary, whether monitor resolution has anything to do with it etc.


texts: ArcGIS JavaScript API Search widget doesn't work

texts:
I am a beginner of ArcGIS JavaScript API. I am creating a simple web map now. For some reason, the search widget of the map doesn't work. When I go to the console, it shows this error:

Here is my source code:

var navigation;
var mainMap;
require([
    "esri/symbols/SimpleFillSymbol",
    "esri/symbols/SimpleLineSymbol",
    "esri/Color",
    "esri/map",
    "dojo/dom-construct",
    "esri/dijit/Popup",
    "esri/InfoTemplate",
    "esri/layers/ArcGISDynamicMapServiceLayer",
    "esri/dijit/Search",
    "esri/layers/FeatureLayer",
    "esri/toolbars/navigation",
    "esri/dijit/Legend",
    "dojo/domReady!"
], function(SimpleFillSymbol, SimpleLineSymbol, Color, Map, domConstruct, Popup, InfoTemplate,
            ArcGISDynamicMapServiceLayer, Search, FeatureLayer, Navigation, Legend ){
    var fill = new SimpleFillSymbol("solid", null, new Color("#A4CE67"));
    var popup = new Popup({
        fillSymbol: fill
    },domConstruct.create("div"));

    mainMap = new Map("mainMap",{
        center: [-112,34.3],
        zoom: 7,
        basemap: "streets",
        infoWindow: popup
    });
    popup.startup();
    var providerInfoTemplate = new InfoTemplate();
    providerInfoTemplate.setTitle("Provider Information");
    providerInfoTemplate.setContent("Provider Name:
 ${Provider}
 Address:
 ${FullAddress}
 Type:
 ${Type_Behavior_Health_Hospital_Health_Center_}
 Intergration Model:
 ${IntergrationModel}
");
    var communityInfoTemplate = new InfoTemplate();
    communityInfoTemplate.setTitle("community Information")
    communityInfoTemplate.setContent("Community Name: ${COUMMUNITY}
")
    var mainLayer = new ArcGISDynamicMapServiceLayer("http://arcgis-cabhp-1819518128.us-west-1.elb.amazonaws.com/arcgis/rest/services/HCIC/160411_interactive_map_jzz_v01/MapServer",{});
    mainLayer.setInfoTemplates({
        0:{infoTemplate: providerInfoTemplate},
        7:{infoTemplate: communityInfoTemplate}
    });
    var mainSearch = new Search({
        enableButtonMode: false,
        enableLabel: false,
        enableInfoWindow: true,
        showInfoWindowOnSelect: false,
        map: mainMap
    }, "searchBar");
    var prividerLayer = new FeatureLayer("http://arcgis-cabhp-1819518128.us-west-1.elb.amazonaws.com/arcgis/rest/services/HCIC/160411_interactive_map_jzz_v01/MapServer/0");
    var sources = mainSearch.get("sources");
    sources.push({
        featureLayer: prividerLayer,
        searchFields:["Provider"],
        displayField: "Provider",
        exactMatch: false,
        outFields:["Provider"],
        name: "HCIC Provider Information",
        maxResults: 6,
        maxSuggestions: 6,
        infoTemplate: new InfoTemplate("Provider Information", "Provider Nmae: ${Provider}
Provider Address: ${FullAddress}
Provider Type: ${Type_Behavior_Health_Hospital_Health_Center_}"),
        enableSuggestions: true,
        minCharacters: 0
    });
    mainSearch.set("sources", sources);
    mainSearch.startup();

    navigation = new Navigation(mainMap);
    myLegend = new Legend({
        map:mainMap,
        layerInfos:[{
            title:"Legend",
            layer:mainLayer,
            hideLayers:[1,2,3,4,5,6,8,9]
        }]
    }, "legendDiv");
    myLegend.startup();
    mainMap.addLayer(mainLayer);
});

function zoomToPrev(){
    navigation.zoomToPrevExtent();
}
function zoomToNext(){
    navigation.zoomToNextExtent();
}
function zoomToFull(){
    mainMap.setExtent(new esri.geometry.Extent(-116.5,30.8,-107.2,37.9));
}


texts:
I used QGIS a lot last year for Dutch cadastral maps purposes. Since December 2020 I notice something has changed when I export maps to dxf. Exports have always been as Vectorfile. Lately they become rasterfiles.
Does anyone know how I could get exports back in vector?
Nothing has changed in the inputfiles (zip)


texts: Updating $geometry in QGIS processing modeler?

texts:
I'm creating a Qgis model with Processing modeler. I can't find how to update the $geometry field from a WKT field. If I complete the step by step process, I used the field calculator with the formula geom_from_wkt("WKT_path"), but it doesn't work within the processing modeler.

Is my approach correct?


texts: DWITHIN ECQL query of WFS is not working properly

texts:
I don't know where I am making mistake. And my URL is returning all values of the feature layer:

queryRequest: function (layer) {
    if(layer!= null || layer!= undefined){
        var center = this.draw.getGeometry().getCenter();
        var radius = this.formatRadius();
        var latLong = this.getLatLong(center);
        var latLongExtent = this.getLatLongExtent(latLong);
        var wfsUrl = layer.getWFSGetRequestUrl();
        console.log("Feature Url :"+wfsUrl);
        wfsUrl = wfsUrl+"&CQL_FILTER=DWITHIN(geometry,POINT("+latLongExtent+"),"+radius+")";
//It should return only one value ..."travel_num":4727.
//wfsUrl ="http://geoserver.saartha.net:7070/geoserver/wfs?service=WFS&version=1.0.0&request=GetFeature&typename=saartha:s03travel&outputFormat=json&srsname=EPSG:3857&CQL_FILTER=DWITHIN(geometry,POINT(76.78399222656249%2013.47409553818946),27.71,kilometers)"

                console.log("Feature Url :"+wfsUrl);
                var request = new ProxyRequest().get(wfsUrl,{
                    version:'1.0.0',
                });
                request.then(lang.hitch(this,"readFeatures"))
                console.log("Request Url :"+request);
            }
        },

## Heading ##

//converting to lat and long.
----------

 getLatLong :function(extentCircle){
    var latLong=[];
    var coord=[];
    coord=[extentCircle[0],extentCircle[1]];
    var convertedVal = ol.proj.transform(coord,'EPSG:3857','EPSG:4326' );
    latLong.push(convertedVal[0],convertedVal[1]);
    return latLong;
    },

getLatLongExtent: function (latLong) {
    var latLongExtent=latLong[0]+" "+latLong[1];
    console.log("LAT-LONG :"+latLongExtent);
    return latLongExtent;
    },

formatRadius: function () {
    var radius;
    if ( this.draw != undefined && this.draw != null)
    {
        var center = this.draw.getGeometry().getCenter();
        var endPoint = [center[0], center[1] + this.draw.getGeometry().getRadius()];
        var sourceProj = this.mapObj.getView().getProjection();
        var c1 = ol.proj.transform(center, sourceProj, 'EPSG:4326');
        var c2 = ol.proj.transform(endPoint, sourceProj,'EPSG:4326');
        radius = this.wgs84Sphere.haversineDistance(c1, c2);
    } else {
        radius = Math.round(this.draw.getGeometry().getRadius() * 100) / 100;
    }

    var output;
    if (radius > 100) {
       output = (Math.round(radius / 1000 * 100) / 100) + ' , ' + 'kilometers';
    }  else   {
       output = (Math.round(radius * 100) / 100) + ' , ' + 'meters';
    }
    return output;
    },
    //reading the features.
    readFeatures: function (response) {
        var feature = this.geoJsonformat.readFeatures(response);
            console.log("feature :"+feature);
        },


texts: QGIS - Have spreadsheet with col from 0% to 100% plus null want 100 bin gradient plus one more for null

texts:
Newbie alert. Have a spreadsheet with no geo; has a col with values in there from 0% to 100%, plus several blanks/nulls.

Ideally would like a 100 step (bin) gradient fill from 0 (100% black) to 100 (100% white) and then color the nulls pale yellow (for oops...missing data)

Have no trouble with CSV import, .csvt files, using the spreadsheet plug in, etc.  Have not had trouble with simple joins. [Other than wondering if I can join keys of different lengths using things like an expression Right(bigfield, 15) to chop down a jumbo key to a normal one in the second table?]

When the column comes in as STRING it messes up the binning because 1, 10, 100 in the land of strings are all adjacent in the sort order. (Or, do I have to tough it out and manually correct? Not bad on 1-10, worse for 1-100.)

When I make the col come in as INTEGER or DOUBLE the blanks/Nulls are assigned 0's and that masks the "missing data". Unacceptable because I need to "see" the missing data and not pretend that its there.

Ideas?

Also, having a really rough newbie time with the expression editor - doing simple things like converting string to numeric, decimal to integer, etc. I'm coming from spending way too many years with Excel and a VBA coding background (very lenient on type conversions, and very easy to get in and out of all the major databases.) I take it the "philosophy" in QGIS is quite different?  When I'm in Expression Editor am I actually in the Python world? Any thoughts on how a VBA coder should be thinking over here? Stuff like single quotes vs double, upper and lower case, fundamental data types, casting, etc...I seem to be stumbling on all. Any "summary" of this type of thing anywhere? Or should I allocate some time to learning basic Python syntax and writing simple Python programs? (If so, do I need a huge IDE like Eclipse or can I have at it with Notepad?)


texts: How can I track/graph arcserver usage statistics?

texts: I’m currently using the ArcGIS Server Properties dialog in ArcCatalog to periodically check the usage statistics for our services.  I’d like to be able to chart those to get a big picture view of what’s going on.  Does anyone know the location on an ArcServer machine where a file with these statistics is being written, ideally something that could be pulled into a spreadsheet?


texts:
I am new to QGIS.

I am trying to import a shapefile (project boundary) into QGIS but when it asks for the CRS I can't find the zone (35) for Finland (in NAD 83).

What am I doing wrong?


texts: Refreshing authentication token in Google Earth Engine?

texts:
I am trying to upload some of GeoTIFF image to GEE using batch image uploader, and receiving an error:

  Incorrect credentials. Probably. If you are sure the credentials are OK, refresh the authentication token. If it did not work report a problem. They might have changed something in the Matrix

I am sure the gmail credential that I typed is correct.

How do I refresh the authentication token in GEE?


texts: Can QGIS Measure Line tool compute correct geodetic distance?

texts:
As I understood, this tool normally shows distance in the project CRS units, not the geodetic meters (despite it actually says the distance is in meters).

However, now I read (Section III, 4 - Measure some distances) that in the Measure Tool Options you can change the "Ellipsoid for distance calculation":

Does it mean it can actually finally compute correct geodetic distance?

However, I don't see this option in the QGIS ver 1.9.0-Master.


texts: Exporting an image with a scale bar in GEE

texts:
I'm trying to export Global Surface Water imagery using earth engine code editor. I want to include a scale bar over the region of interest in the exported image... how easy is this to accomplish? I've already extracted an image as a TIFF, would it be easier to open the TIFF in ArcGIS to add the scale bar or is there some javascript that could be added to address this? Not sure where to begin!

I note that GEE shows a little scale bar already on the lower right hand corner, but this doesn't get preserved in the image I download.

Existing image port is below:

var gsw = ee.Image('JRC/GSW1_0/MonthlyRecurrence/08');
var recurrence = gsw.select('monthly_recurrence');

print(gsw);

//////////////////////////////////////////////////////////////
// Calculations
//////////////////////////////////////////////////////////////

// Create visualization layers.
var recurrence2 = recurrence.updateMask(recurrence.divide(100))
print(recurrence2);

var recurrence2RGB = recurrence2.visualize({
 min:0,
  max:100,
  palette: ['#FF7F27', '#99D9EA']
});

//////////////////////////////////////////////////////////////
// Map Layers
//////////////////////////////////////////////////////////////

print('Image to export: ', recurrence2RGB);
Map.addLayer({
  eeObject: recurrence2RGB,
  name: 'OCCURENCE IN AUGUST',
  shown: true
});

Export.image.toDrive({
  image: recurrence2RGB,
  description: '08',
  scale: 30,
  region: geometry,
  crs: 'EPSG:3857'
});


texts: How do I count the number of vector lines that pass though a given raster cell?

texts:
I am working on a project (currently in GRASS7 though I'm not against exporting to Arc or R or other software if necessary) where I am generating thousands of shortest cost paths (using the r.walk and r.drain modules) to various points around my map of Papua New Guinea. What I am interested in is figuring out how to take the generated vector lines and create a raster that counts the number of lines that pass through each cell of that raster. This will give me an estimate of which locations on the island are most likely to be located on a shortest path from one point to another.

One way to do it would be to (using python)

Create a raster initalized to 0
use a python script to run though and rasterize each line individually and
add each line to the raster using raster calculator to create a count.

The problem with this approach is that for the tens or even hundreds of thousands of lines I'm using, it takes a very very long time. I would like a way to directly count the number of lines. Something like http://grasswiki.osgeo.org/wiki/Count_points_in_raster_cells
except with lines.


texts: How to store tile maps into databases and use them for web mapping?

texts: I am making a web map using leaflet and a lot of tile maps from raster files. When using gdal2tile.py to convert these rasters into tile maps, it will automatically create a bunch of directories and thousands of PNGs. Since I will have thousands of rasters to be tiled I think it might not be a good choice to save them this way. I am thinking is it possible to save them into database and query them when needed? I am using Postgresql and PostGIS for other functions of the map. So if Postgresql can handle this, it would be the best.


texts:
I know there is a similar question, but it might be dated, and didn't completely answer my question.
I know of two plugins to export your QGis project, QPackage and QConsolidate. I'm not very familiar with either. The lather isn't available in the plug-in space in QGIS either, so that makes it a bit more complicated. But what I need:

to export all my layers, with multiple source files over multiple maps
if possible, save the styles as well at the same time
layers are grouped like this:

Project area 1

Work area 1

Coupe
Put
Vondst

Work area 2

Coupe
Put
Vondst

Project area 2

Work area 1

Coupe
Put
Vondst

etc.
So I have layers that have the same name, but fall in a different group/subgroup. I want to export all this at the same time. But I think i'll lose the layer relations (where they are placed in a group), and think the name similarities will give problems with file names.
Is what I want possible with the plugins? Or should I just export everything by itself?
Edit:
So far I've renamed all my layers so they have unique names.
With QConsolidate the layers don't gert exported, just the .qgs, so no luck there.
With QPackage everything gets exported, but i do get a warning messege, though I don't know what for. The only thing going wrong there is that the connections with the source file get broken when I move the folder. And since the .qgs is copied entirely it also shows the layers I dindn't export. But the plusside is that the styles get saved as well. So I think reaconnecting the layers with the exported files should do the trick.
This is all very time consuming though, so I still hope someone knows a faster way.


texts: OpenLayers tiles from reprojected mapdata mixes zoom levels

texts:
We have a map tile layer on OpenLayers where original projection is EPSG:3067 and it is reprojected to EPSG:3395 projection. The map is coming there quite ok, it has not the tiles as straight squares, but still ok for us.

The problem is that in the northern part of the map the tiles that are on map become from two layers. There is a horizontal line with many curves that separates the tiles in north with another zoom level and south with another.

Preview on GeoServer shows the tiles all at same zoom level, openlayers code based map does not. We had an issue about same kind when we had custom zoom levels on map, but removing them the divider appeared in north. Previously it was much souther.

What could cause the layer to take tiles from different zoom level per tile basis, having two different levels in use on the whole map, division between north and south parts.

The issue would be kinda ok, if there weren't place names on map. Now they're cut and missing partly.

Layer is initialized with a ol.source.TileWMS taking in layer name, format (image/png), transparency and a version.

The layer is a ol.layer.Tile with that tileWMS as source, title and visibility.

All that seems so basic, there is at least no hacks on there. What could cause the issue?

Edit1:

It appears the GeoServer returns tiles zoomed differently. I took the tile request urls and in the tiles from these I found tiles looking like they were zoomed to different levels.

It requests tiles 256x256 from different bounding boxes in those queries, somehow they happen have different zoom levels.

May the conversion on fly somehow mix the thing?

Edit2:

Newest finding is that the original data is a layer group, so that may cause the switch of level, if some requests give data from another group from some reason.


texts: Problems importing raster to GRASS from QGIS (error: Datum ETRS89)

texts:
When I try to import raster from QGIS to GRASS with r.in.gdal.qgis an error of projection appears.

Raster is in EPSG::25831 - ETRS897UTM zone 31N and mapset is supposed to be defined with the same.

Projection of dataset does not appear to match current location.

Location PROJ_INFO is:

name: UTM Zone 31, Northern Hemisphere
ellps: grs80
proj: utm
zone: 31
towgs84: 0,0,0,0,0,0,0
no_defs: defined

Dataset PROJ_INFO is:

name: ETRS89 / UTM zone 31N
datum: etrs89
ellps: grs80
proj: utm
zone: 31
towgs84: 0,0,0,0,0,0,0
no_defs: defined

ERROR: datum

In case of no significant differences in the projection definitions, use the -o flag to ignore them and use current location definition.

Consider generating a new location from the input dataset using the 'location' parameter.

Finished with error

I created the mapset with the same coordinate system: EPSG: 25831 - ETRS897UTM zone 31N. But looking information of my current location (with g.proj.print) no datum is printed:

-PROJ_INFO-------------------------------------------------
name       : UTM Zone 31, Northern Hemisphere
ellps      : grs80
proj       : utm
zone       : 31
towgs84    : 0,0,0,0,0,0,0
no_defs    : defined
-PROJ_UNITS------------------------------------------------
unit       : Meter
units      : Meters
meters     : 1

Furthermore, if I try to import the raster to GRASS considering the advanced option of Override projection check, the final result appears in a new coordinate system, EPSG::3043-ETRS89/UTM zone 31N (N-E) and with different values from the original one.

I am using QGIS 3.4.12 with GRASS 7.6.1 in Windows 10.

It is possible that EPSG::25831 - ETRS897UTM zone 31N is not supported as a Location (from a mapset)? Could it be solved?


texts: I am currently working with QGIS and for my project I need to start QGIS2Web automatically. Is there a way to start this plugin using the Python Console in QGIS?


texts: Extracting coordinates of segments in raster from TIFF image?

texts:
I'm working with the attached tiff file and I am trying to extract the xy coordinates of the green lines.

so far I saved the image as raster stack, and this specific image is from a layer of that raster

file_path <- "MyCells.tiff"
rs <- stack(file_path)
plot(rs[[2]])

however, I don't know how to extract specifically those green lines. I trid to follow these two questions:

Extracting value of raster with coordinates using R?
Extracting GeoTIFF data with coordinates using R?

but I couldn't find an answer to my question. google sends me to image manipulation files.


texts: automatically generate layouts based on field values in arcgis pro

texts:
I am using ArcGIS Pro 2.5.1.

I am creating maps of bird species distributions for 127 species, all over the same geographic area. I have a point feature class where each point is a bird sighting, with different symbols to represent breeding status.

I would like to create a layout for each species, using the same symbology and spatial extent for each layout. Effectively, I want to generate an identical layout for each unique value in the 'species' field of my point feature class. I thought this should be fairly straightforward to automate, but haven't been able to figure out how to do it. My first thought was creating a Map Series, but it seems my only option is to define the different maps in the series by spatial extent, rather than field values.

Is there a way to create a series of layouts defined by field values? I haven't much experience with Python so hoping that's not the only way to do this, but happy to learn if it is.


texts: Why reprojection GeoJSON with new.ol.View - not working?

texts:
I try to reproject GeoJSon file in OpenLayers 3 from :

Geojson link
With these JS in head :

<script src="js/jquery-2.2.3.min.js" type="text/javascript"></script>
     <link rel="stylesheet" href="http://openlayers.org/en/v3.15.1/css/ol.css" type="text/css">
<script src="http://openlayers.org/en/v3.15.1/build/ol.js"></script>
<script src="//cdnjs.cloudflare.com/ajax/libs/proj4js/2.2.1/proj4.js"></script>

With this code to set the new ol.View with
newProjCode :

var newProjCode = 'EPSG:27700'
proj4.defs(newProjCode, '+proj=utm +zone=31 +ellps=WGS84 +datum=WGS84 +units=m +no_defs');

var newProj = ol.proj.get(newProjCode);
var fromLonLat = ol.proj.getTransform('EPSG:4326', newProj);
var newView = new ol.View({
  projection: newProj
});
map.setView(newView);

But after setting that new ol.View i see only white screen.

Where is the problem ?


texts: arcpy: can I use SelectLayerByAttribute to compare fields from joined tables?

texts:
I can't figure out the syntax to select table records where one field's value exceeds that of another. The fields are from a pair of table views connected w/AddJoin(). One view is of a fgdb table and the other is of a point feature class in a separate fgdb. My ArcGIS10.0 syntax looks like:

arcpy.MakeTableView_management(r'C:\Workspace\wk_NHDPlus\NHDPlus\NHDPlus_national.gdb\POIs\GagesII', 'gagesTV', '"REGION" = \'02\' AND "COMID" = 0')
arcpy.AddJoin_management('msrTV', 'ComID_Measures_REACHCODE', 'gagesTV', 'REACHCODE', 'KEEP_ALL')
arcpy.SelectLayerByAttribute_management('msrTV', 'NEW_SELECTION', '"GagesII.MEASURE" >= "ComID_Measures02.ComID_Measures_BOTMEAS" ')

I've also tried:

newName = arcpy.AddFieldDelimiters(r'C:\Workspace\wk_NHDPlus\NHDPlus\NHDPlus_national.gdb\ComID_Measures02', "ComID_Measures02.ComID_Measures_BOTMEAS")
arcpy.SelectLayerByAttribute_management('msrTV', 'NEW_SELECTION', '"GagesII.MEASURE" > '+newName)

The result in both cases is the entirety of the original 'msrTV', which I know is incorrect.

Is what I'm trying to do possible? Hopefully I'm just fouling up the syntax.

* update: also tried:

newName = arcpy.AddFieldDelimiters(r'C:\Workspace\wk_NHDPlus\NHDPlus\NHDPlus_national.gdb\ComID_Measures02', "ComID_Measures02.ComID_Measures_BOTMEAS")
oldName = arcpy.AddFieldDelimiters(r'C:\Workspace\wk_NHDPlus\NHDPlus\NHDPlus_national.gdb\ComID_Measures02', "GagesII.MEASURE")
sqlquery=oldName+'>='+newName
arcpy.SelectLayerByAttribute_management('msrTV', 'NEW_SELECTION', sqlquery)


texts: Display attributes for created features disappear upon exiting "Edit"

texts:
I have a very simple map with a land grid. After adding a new polygon, the symbology shows up appropriately, and remains through Save Edits, but disappears as soon as I hit Stop Editing. The polygon remains, both visibly (if I select it) and in the Attribute Table, but the symbology changes to blank, and will only show up if I check "All Other Values" in the Symbology tab.

Things I have tried, to no avail:

Copy/pasting the value within the field from a polygon that shows up correctly into the appropriate field of the new polygon. Polygons with correct symbology have blank fields. I've also tried deleting anything that might have gotten into the field (spaces..?).
Many of the suggestions on this post: Created features disappear when saving edits?, including checking Feature Exclusions, and there are no Definition Queries.
Copy/pasting the Attributes from one "good" polygon to the "bad" polygon. It changes the symbology to the correct color, keeps it after Save Edit, and then reverts after Stopping Editing.
Delete feature templates, and then create one using a "good" polygon that shows up as expected. Use that template to draw a new polygon. Same result, stays through Save Edit, disappears after Stop Editing.
Copy/paste a "good" polygon, then use Editor to Move Vertices over to where I want the new polygon to be.

Also, the color/symbology doesn't seem to disappear until I zoom out. But when I zoom back in, it doesn't reappear. I don't have any extents set on this layer.

I am using ArcGIS 10.4 for Desktop, and I am very much a beginner (I know how to navigate, but my vocabulary is terrible).


texts: Exporting specified layers in ArcMap using ArcPy?

texts:
I am trying to export specific layers from an ArcMap .mxd. From the TOC below I would like to export "irrig_2, citylyr, boundlyr" as a GIF. Once this is done the script would then export "irrig_1, citylyr, boundlyr" as another GIF and so on. If possible, the script would also export a legend that corresponds to the layers being exported. I'm not sure if it's a good place to start, but I'm trying to work off the code below.

import arcpy

mxd = arcpy.mapping.MapDocument(r"E:\IrrigatedLands\FC_test\irrig2.mxd")
df = arcpy.mapping.ListDataFrames(mxd, '')[0]

GIFPath = r"E:\IrrigatedLands\FC_test"

for lyr in arcpy.mapping.ListLayers(mxd, '', df):
        if lyr.name == "citylyr":
            lyr.visible = True
        if lyr.name == "boundlyr":
            lyr.visible = True
        if lyr.name == lyr:
            lyr.visible = True
            arcpy.mapping.ExportToGIF(mxd,GIFPath+"\\" + lyr.name + ".gif")
            lyr.visible = False
arcpy.RefreshActiveView()
del mxd


texts: Config.json shown instead of the setting UI

texts:
I copied the setting page from the Demo widget and I try to adjust it to suit my widget. My config.json contains the configText variable and several others as well. When I open the config window for my widget in app1, I see only the config.json contents, not UI as for other widgets. I can edit the JSON code freely.

So far I found no way to use any more complex debugging options or even how to debug the settings.js through FireBug (I can debug Widget.js, which works), so I can't confirm that the setting code is not even loaded, but it doesn't seem to respond to changes in my setting.js or setting.html files. The files seem to be on their places (config is read from WABRoot/server/apps/6/configs/MyWidget/config_MyWidget.json and the widget from WABRoot/server/widgets/MyWidget/), though I might have forgotten to check something. The setting code is exactly like in the Demo widget, except I renamed baseClass from jimu-widget-demo-setting to jimu-widget-my-widget-setting.

What else can I do, and especially, why do I get the json instead of the UI?

I reposted this on the ESRI support forum.


texts: Unable to generate token 'referer' must be specified

texts:
I have created a HTML page which contains a feature layer that is published in the ArcGIS Online. The feature layer is shared with organization. When I run the page it is asking for the username and password, and whenever I enter the credentials it shows the username or password you entered is incorrect, but I know my username and password is absolutely correct.

When I checked in network, in generateToken file it gives error like:
message: "Unable to generate token."
details: ["'referer' must be specified."]

And when I tried the same code in ArcGIS Online sandbox it worked fine.


texts: advice on browsing multiple points with leaflet - query only the points in the viewport

texts:
This a very basic question, but I don´t quite know where to start

I have a postgis table with aproximately 15.000 points hosted on CartoDB

And I'm developing an application that should let the user navigate a map, and show the points that are inside that area, with several filter applied

I was thinking about building a query like this to pass to the web service

select * from my_table
where
  ST_Contains(
    the_geom,
    ST_GeomFromText(
      'POLYGON((-71.177 42.390,-71.1776 42.3903, -71.1776 42.3903,-71.17758 42.3903))'
    )
  ) and
  [more filter conditions]

I just need to know the four points of the view port from leaflet, and update it everytime the user moves or zooms...

I guess this it just what the map view on cartoDB does...

Can anybody give me a hint on how to achieve it, or at least tell me if my approach is completely wrong.

Thanks a lot


texts: Unable to establish a persistent ODBC connection to an access mdb database

texts:
I’m an MS Access diehard and I’m running a 64bit version of Win 10 on my PC. I’ve installed both the 32bit (2.14.7) and 64bit (2.16.3) versions of QGIS on it, and MS Office 2016.  I would like to establish a persistent ODBC connection to an access mdb database but can’t get this to work on either on versions of QGIS.

Using the EVis plugin I have no issues in terms of linking to mdb, and identically accdb, databases on both the 32bit and 64bit versions of QGIS.  However, try as I may I can’t form a successful link via QGIS’s add vector layer ‘database’, ‘ODBC’ setting – the only thing I get is a database table header, its field structure, but no data whatsoever and the name assigned to the layer when I attempt this ends with ‘Unknown’.

It has been suggested that a modification to QGIS's bat file should overcome the issue I'm experiencing, as per adding the following two lines:

set OGR_SKIP=ODBC
set PGEO_DRIVER_TEMPLATE=DRIVER=Microsoft Access Driver (*.mdb, *.accdb);DBQ=%%s

I've done this, but in my two installations of QGIS the second line of that code doesn't appear to have any effect, since when I drag an mdb file into QGIS's workspace I get an invalid file format error.

How I can resolve this issue and establish in QGIS a persistent ODBC connection to an access mdb database?


texts:
I am writing a plugin to import some files into qgis. The files initially have coordinates in lat lon, but I want to reproject them to a local projection during the import.

So, I can set the CRS on the layer I am importing in and and I know that my imported points are in WGS 84, so what I want to do is similiar to what I would have done in postgis with

ST_Tranform(ST_SetSRID(ST_MakePoint(15, 60),4326),32633);

to import the point into a layer WGS84 UTM zone33N

But in qgis, how can I make a PointXY() with a CRS and then reproject it to another CRS?

(I know I can read everything into a layer with CRS 4326 and then reproject the entire layer into a new layer, but I would prefer to do it without having to make an extra layer)


texts:
In China address is usually like: No.few homes, building number, which unit

or:

No.23,Unit 2,Bldg.1 47th Yard, Fengchan Road Jinshui District

Does OSM support this operation?


texts: ArcGIS/Adobe Illustrator/Adobe Photoshop: Create blend between two shapes

texts:
I have a landuse layer as an ArcGIS shapefile which looks like this (here only a small sample part)

The round shape of the lightgreen area results from a buffer and as the borders of this feature are not exactly clear I would like to create a blend/gradient change for the final cartographic layout (only between the round light green 'buffer' shape and the darkgreen area). This landuse layer is in vector format (ArcGIS shapefile).

I was wondering if there is an easy way to achieve this either in ArcGIS 10, Adobe Illustrator or Photoshop (would include a Polygon to Raster conversion first and therefore lead to information loss as it is rasterized)?

Later I would like to drape the landuse layer over a hillshade with burned in rivers using a transparency of about 20% like this.


texts: fixing polygons with "ring not closed" errors

texts:
I have a shapefile with a single Polygon feature. This feature has a "ring not closed" geometry error (seen here in QGIS):

I tried using buffering and the GRASS v.clean command as suggested in other answers, but no luck. I also can't find the gap by inspecting the feature visually (it looks perfect).

I need a valid geometry here. Ideally, I'd like to have a Python/command line solution. Any way of doing this manually in QGIS/ArcGIS would also be great.

UPDATE: this shapefile can be downloaded here.


texts: Checking if polygon exceed in at least one direction another polygon in ArcMap?

texts:
I have a shapefile with some polygons,I need to compare these polygons with another set of polygons in another shapefile and see if the second are larger in AT LEAST one direction.

How can I do it?

I'm using Arcmap in the trial version for the moment

As you can see, I have a set of yellow polygons that have to be compared with the blue ones (that are several and of really different dimensions): I need to know for each yellow polygon if the blue polygon that contain his centroid is TOTALLY included in the yellow one or is larger for at least one direction

What I mean is that if the blue polygon that contain the centroid of the yellow one is larger or not. Maybe it's easier with a graphical support:

Case 1: the blue is greater than the yellow
Case 2: the blue is smaller than the yellow
Case 3: the blue is greater than the yelloe in one direction (so for my purpose it's greater

I need to refer to the polygons in which the centroid fall, because the blue are consecutive (as you can see from the previous image) and the yellow can exceed in another one.


texts: Are there templates for adding a search filter on a CartoDB embed?

texts:
I'm looking for any examples of an embedded CartoDB map with a search filter built in. So, a map of schools, say, with a pulldown that will let you select one school and then zoom to it.

Or a map of schools and a pull down that will let you toggle between Elementary, Middle and High schools.

I know I can use filters or a custom query but I'm looking for something viewers of the final map can play with themselves.

Here's the deal: I'm trying to figure out what it would take to make a map that a colleague wants. If I have to do a lot of JS work, that's a bigger project. If there's a cut-and-paste solution that makes it easier. And if there's a GUI, they can make it themselves.


texts: Create shapefile from millions of geocoded addresses

texts:
I am trying to make a map of administrative units in a state. I have a list of ~10 million addresses with a variable indicating which unit they fall into.  My plan is to geocode the addresses using an Amazon EC2 instance or similar geocoding solution, map the points and then draw polygons around the sets of points.

Given my relative inexperience with GIS, I am wondering if there is a more sensible way to approach this problem.  The addresses are quite clean, with zip codes attached (sometimes zip+4), and 99% or more of the addresses should be clustered together correctly in a single contiguous polygon.  But I anticipate that after geocoding some addresses will fall outside the logical bounds of the polygon and need to be corrected post hoc.  Each unit will have ~1200 addresses geocoded on average, with significant variance.  My address list represents approximately 40% of all addresses within the units, so the points I will have to work with are quite dense.

Is this the only way to approach the problem?  Note I am a grad student, so any solution needs to be cheap or free (my budget is under $300).

Edited to add: All addresses are in a single US state. I am working in ArcMap 10.4 and even figuring out how to use an API has been tough. My addresses are just in a .csv file and I eventually want a .shp file. I am a statistician, not a GIS professional, unfortunately.


texts: Modify coordinate string in parsed KML file

texts:
I am new to Python. I found this code on the internet to parse names and coordinate strings from KML files and write them into a .csv file. Works excellent! However, for my purpose, I need to convert the coordinate string to GeoJSON as output. For my GeoJSON format I need to remove the height/accuracy from each coordinate. So I only need the lat,lon coordinates and want to remove the height/accuracy. Besides that I need to enclose each individual lat,lon coordinate with square brackets []. So, for example, I now have this output format for my coordinate string:
-89.42462055297086,16.35478756523923,507.7306687727497
-89.41330216185102,16.33874177258091,521.3689105782825

And need to convert this coordinate string to this format:
[-89.42462055297086,16.35478756523923],[-89.41330216185102,16.33874177258091]

Below is the code I am using now. I want to add some code to create the GeoJSON coordinate string format I need. How can I do this?

import xml.sax, xml.sax.handler

inputname = 'Test4.kml'
kml = open(inputname, 'r')

class PlacemarkHandler(xml.sax.handler.ContentHandler):
    def __init__(self):
        self.inName = False # handle XML parser events
        self.inPlacemark = False
        self.mapping = {}
        self.buffer = ""
        self.name_tag = ""

    def startElement(self, name, attributes):
        if name == "Placemark": # on start Placemark tag
            self.inPlacemark = True
            self.buffer = ""
        if self.inPlacemark:
            if name == "name": # on start title tag
                self.inName = True # save name text to follow

    def characters(self, data):
        if self.inPlacemark: # on text within tag
            self.buffer += data # save text if in title

    def endElement(self, name):
        self.buffer = self.buffer.strip('\n\t')

        if name == "Placemark":
            self.inPlacemark = False
            self.name_tag = "" #clear current name

        elif name == "name" and self.inPlacemark:
            self.inName = False # on end title tag
            self.name_tag = self.buffer.strip()
            self.mapping[self.name_tag] = {}

        elif self.inPlacemark:
            if name in self.mapping[self.name_tag]:
                self.mapping[self.name_tag][name] += self.buffer
            else:
                self.mapping[self.name_tag][name] = self.buffer
        self.buffer = ""

parser = xml.sax.make_parser()
handler = PlacemarkHandler()
parser.setContentHandler(handler)
parser.parse(kml)

def build_table(mapping):
   sep = ';'

   output = 'Name' + sep + 'Coordinates' + sep + 'test column\n'
   points = ''
   lines = ''
   shapes = ''
   for key in mapping:
       coord_str = '{"type":"FeatureCollection","features":[{"type":"Feature","geometry":{"type":"Polygon","coordinates":[[[' + mapping[key]['coordinates'] + ']]]}}]}' + sep

       if 'LookAt' in mapping[key]: #points
           points += key + sep + coord_str + "\n"
       elif 'LineString' in mapping[key]: #lines
           lines += key + sep + coord_str + "\n"
       else: #shapes
           shapes += key + sep + 'coord_str\n'
   output += points + lines + shapes
   return output

outstr = build_table(handler.mapping)
out_filename = 'output-result.csv'#output filename same as input plus .csv
f = open(out_filename, "w")
f.write(outstr)
f.close()
print(outstr)
°°°


texts: How to get a list of all data sources used in a ArcGIS Pro project?

texts:
Is there a way to generate (or access) a list of all data sources used in a ArcGIS Pro project?

I would like to avoid going through each layer's property dialog.


texts: Merge uniquely named point shapefiles in 50,000 different folders into a single layer

texts:
I have point shapefiles in 50,000 different folders (each shapefile file name is unique) - about 33,400 shapefiles (not all folders will have a point shapefile). I like to open the point shapefile in each of the folder into QGIS Desktop 2.2.0 map frame and then hopefully run a merge on them to get a single point shapefile layer. I tried drag-dropping from Windows Explorer to QGIS, and QGIS froze initially and continue to add the point shapefiles for several hours and didn't seem to stop. I found the Load Them All plugin (http://goo.gl/HbJhhS) that seem to specify the feature type and also specify a string filter, this plugin also didn't seem to work...

I was thinking that had I planned well, I could've automated and copied the point shapefiles into another single folder as soon as they were created and I could 've dragged-dropped into mapframe from this single folder. I don't know if this would've had worked. My analysis run took 18 days to run, so I don't want to re-run just to do that. I tried copy-pasting with Windows Explorer and it hung, may be because it was 100,000+ files...

Would you know any trick to do this?


texts: Using Proportional Symbol Map vs Graduated Symbol Map?

texts:
I am having trouble understanding the difference between a Proportional Symbol Map and a Graduated Symbol Map, or why it matters.

How do these two differ from each other?


texts:
Is there a simple way to convert an AVL file into an SLD or QML file?
deegree does not work an my laptop (Java problems)

Or is there a way to open the original AVL file in QGIS?


texts: Do any of the SQL Server 2008 R2 SP 2 cumulative updates contain the MAXDOP spatial index fix?

texts:
I am troubleshooting ArcGIS 10.1 performance issues with SQL Server 2008 R2 SP 2 (10.50.4000) and the issue appears to be this MAXDOP issue:
https://connect.microsoft.com/SQLServer/feedback/details/625896/spatial-index-not-used-on-multiproc-machine-unless-maxdop-set-downward

The issue is that the query optimizer is not using the spatial index when it should. Setting a OPTION (MAXDOP 1) query hint on a SQL statement will cause the spatial index to be used as will an index hint, but neither are an option for the Arc-clients.

I've read that SQL Server 2008 R2 SP 1 - Cumulative Update 3 and higher has the fix but I'm seeing it at SP2.  And I don't see the fix listed in the documentation for any of the CUs, but I'm hoping it is undocumented and present in one of the cumulative updates.  Has anyone had this issue at SP2 and fixed it with one of the SP2 cumulative updates?


texts: How to create an elevation profile in google earth from a kml?

texts:
I have a KML file of a stream that I originally created as a line shapefile in ArcMap (one stream extracted from a larger watershed data set). I need to create an elevation profile of this stream. While I do have DEMs for this area, unfortunately I do not have spatial analyst or 3D analyst in ArcMap, and my employer is unwilling to purchase these tools.

Since I do not need super detailed elevations, my idea was to create a simple elevation profile in Google Earth, which can be done from a path. However, that option is "grayed out" for my KML file in google earth. Is there a way to create this elevation profile for a KML, or potentially convert the KML to a path?

This stream is very windy and while I have tried tracing in in google earth to create the path a few times, my traced path comes out very wonky, and would be fairly inaccurate.


texts: Simple way to plot netCDF variables on Mollweide projection in R

texts:
Is there a very simple way to plot up variables from a netCDF file that are stored in a regular lat-lon grid (-90to90˚N, 0-360˚E) onto the now more popular Mollweide projection?

For example say I have something like the following:

lon <- seq(-90,90,2.5) # really these data would come from the netCDF file...
lat <- seq(0,356.25,3.75)
myvar <- runif(length(lon)*length(lat))
mymat <- array(dim=c(length(lon), length(lat)), myvar)

Which I can plot up easily with something like image:

image(lon, lat, myvar)

But what I'd like to do is convert the myvar array and reproject it so that I can plot it in the Mollweide projection:

... convert to Mollweide projection ...
image(myvar_converted)

I'm looking for a very generic approach and ideally something that would allow me to add on points, lines etc on to the plots using base functions if possible.


texts: Split many multi-band files into single-band files using PyQGIS?

texts:
I wanted to know how to take many multi-band (RGB) files in a folder and extract a single band from each one, putting everything in an output folder with the single band transformed files.

How should I proceed without extracting every single file?


texts:
I've built a map in Leaflet which has a really large GeoJSON file. Currently the only way I can get it to work is if I severely simplify the GeoJSON file down from 100mb to 25mb. I'm looking at using the mapbox gl leaflet binding as a way of solving the problem of dealing with such a large file. However, I can't figure out how to interact (on hover and click) with the layer I'm pulling in using Mapbox GL.

This is my code initialising my map and pulling my layer in from Mapbox. The style includes a basemap and a vector tile polygon layer (originally uploaded to Mapbox from a GeoJSON file):

var map = L.map('map', {maxZoom: 17}).setView([37.8, -96], 4);

var gl = L.mapboxGL({
    accessToken: 'myAccessToken',
    style: 'mapbox://styles/myStyle',
    onEachFeature: onEachFeature
}).addTo(map);

Here's the onEachFeature function:

    function onEachFeature(feature, layer) {
        var popupText = "";
        var labelText = "";
        if (feature.properties && feature.properties.GEOID10) {
            labelText = "Zip Code: " + feature.properties.GEOID10.toString();
        }
        layer.bindLabel(labelText).addTo(map);
        layer.on({
            click: highlightFeature,
            mouseout: resetHighlight
            //click: zoomToFeature
        });
    }

    function highlightFeature(e) {
        var layer = e.target;

        layer.setStyle({
            weight: 5,
            color: '#666',
            dashArray: '',
            fillOpacity: 0.7
        });

        if (!L.Browser.ie && !L.Browser.opera) {
            layer.bringToFront();
        }

        info.update(layer.feature.properties);
        map.fitBounds(e.target.getBounds());
    }

    function resetHighlight(e) {
        geojson.resetStyle(e.target);
        //info.update();
    }

I don't know how to implement the on click functionality to interact with the vector layer being brought in by Mapbox GL.


texts: How to convert MDB to DXF with labels using ogr2ogr?

texts:
I have to convert some 700 MDB files (parcel layer). I want to generate DXF from from them. I tried
ogr2ogr -f "DXF" d:\01Chitapol_1Ka.dxf d:\01Chitapol_1Ka.mdb
This gives me the lines. I want the parcel numbers as labels in the DXF file. This post mentions using HEADER and LABEL. I couldn't find examples of using this.


texts: I have many point features in close distance to each other and when I label them, the labels overlap. Is there a way that QGIS automatically seperates the lapels from each other and resolves the overlap?


texts: Definition of "most northwesterly point" on a landmass?

texts:
If we are using a 2D representation of a landmass, with north represented as up and east as right, then we can conceptualise the most northwesterly point as being the last point we reach when we place a ruler on the map from bottom left to top right and move it diagonally upwards perpendicularly to its length. So far, so simple.

But how do we define the most northwesterly point of a landmass if our model of the earth is a sphere or a spheroid?

Let us assume that the point is in the northern hemisphere. One definition I have thought of is that the most northwesterly point on the landmass is the one from which the journey to the equator along a geodesic leaving in a southeasterly direction is longest.

Is there a better definition?


texts:
Recently I have been trying to use the TIGER WOIS package which looks really useful. I'm trying to use the Flood mapping tool with Sentinel-I data, but I encountered a problem right away.

In the first step of the pre-processing (WOIS Workflow 1c - Sentinel I) I get the following error: "Failed to initialize velocity engine". In the History and Log I see that I have errors for all SAGA and Workflow algorithms (error in line...) and this is probably why I run into this mistake (I guess?). Did I do something wrong with the installation of Tiger WOIS or should I make some changes in my settings somewhere?


texts: Update EVLR header info for .LAS file using laspy

texts:
I'm trying to change the header info for a LAS file using the python 2.7 (64-bit) laspy module (version 1.2.5). I'm attempting to follow the example shown here (in the section "Writing Data + EVLRS"): http://laspy.readthedocs.io/en/latest/tut_part_3.html

However, when I update the Extended Variable Length Record (EVLR) info with the new content, I get this error message:

  Error text: ValueError: offset must be non-negative and no greater than >buffer length (1094)

This error shows up on line 151 of laspy/base.py:

    self._pmap = np.frombuffer(self._mmap, self.pointfmt,
                offset = self.manager.header.data_offset,
                count = self.manager.header.point_records_count)

The "buffer length" value of 1094 seems to be the same as the length of the EVLR text record, while the "data_offset" value above is a much higher value of 4151.

import numpy as np
import laspy
import shutil
import datetime

print('Update started: {:%Y-%m-%d %H:%M:%S}'.format(datetime.datetime.now()))

 # Set input / output filenames
evlrContentFile = "./LAS_DATA/evlrContentFile.txt"
inFileName = './LAS_DATA/infile.las'
outFileName = './LAS_DATA/outfile.las'

# Get EVLR string content
with open(evlrContentFile, "r") as evlrFile:
    evlrString=evlrFile.read();

# make a copy
shutil.copy(inFileName,outFileName);

# open the outfile for modification
outFile_v14 = laspy.file.File(outFileName, mode = "rw")

# create a new EVLR header entry
new_evlr = laspy.header.EVLR(user_id = 10, record_id = 2,
                VLR_body = evlrString)

# update the outfile EVLR header entries
old_evlrs = outFile_v14.header.evlrs
old_evlrs.append(new_evlr)
outFile_v14.header.evlrs = old_evlrs
outFile_v14.close()
print('Update Completed: {:%Y-%m-%d %H:%M:%S}'.format(datetime.datetime.now()))

I want to modify the existing files with the new header info with a minimum of processing, so I want to avoid using laspy to copy the contents of the LAS points to the new output file if possible. Hence why I am using the "rw" modify option instead of "w" to create the output file.


texts:
New to using GeoPandas instead of Pandas. I am trying to export a GeoPandas dataframe to a PostgreSQL database (with PostGIS extension).

The import is from the same database as I try to export to using:

gpd.GeoDataFrame.from_postgis(sql, con, geom_col='geom')

The CRS seems to be imported as well as srid=25832. When I print gdf.crs it yields as expected:

{'init': 'epsg:25832'}

A sample of the geodataframe looks like this, where the last column geom (the geometry) seems to be the problem child:

row_num     model  ...  mean_res                            geom
0          1  644862.0  ...  37.40605  POINT (573187.125 6131561.456)
1          2  644862.0  ...  36.36506  POINT (573189.374 6131560.369)
2          3  644862.0  ...  35.42172  POINT (573191.623 6131559.282)
3          4  644862.0  ...  34.90574  POINT (573193.873 6131558.194)
4          5  644862.0  ...  34.36104  POINT (573196.122 6131557.107)
..       ...       ...  ...       ...                             ...
116      117  644862.0  ...  41.61498  POINT (573448.099 6131435.305)
117      118  644862.0  ...  41.26093  POINT (573450.350 6131434.217)
118      119  644862.0  ...  41.39089  POINT (573452.600 6131433.129)
119      120  644862.0  ...  41.49015  POINT (573454.851 6131432.041)
120      121  644862.0  ...  41.56080  POINT (573457.101 6131430.954)

My problem is when I try to export the dataframe using:

df.to_sql(tablename, con=engine, schema=schemaname, if_exists='replace', index=False)

The engine is from the SQLAlchemy lib using psycopg2, which is working for a "normal" Pandas dataframe (i.e. without geometry)

I get the error:

"Table could not be edited: databasename.schemaname.tablename
'GeometryDtype' object has no attribute 'base'"


texts:
Is it possible to hide some icons from QGIS toolbar(s) which are not used often?
I tried with Settings > Interface Customization but I didn't find an option for removing a single icon, however here is possible to remove icon groups but I would like to remove individual icons.
Is there such an option(s)?


texts:
I've created marker cluster in this example and now I want to show the same thing with data from GeoJSON, but it doesn't work. I'm trying to reproduce same situation in this example, but marker cluster doesn't show, there are only markers.

Does anybody knows why Marker Cluster in the second example does not show? Did I miss some attribute in GeoJSON?


texts: Adding columns with x/y data to SpatiaLite Layer?

texts:
I am quite new at working with QGIS.

My Spatialite Database looks like this (arch_locations.sqlite)

Fields:
pkuid
short_ID
name
geometry

I´d also like to have two additional columns in that table in QGIS (in the Spatialite Database) where it displays the X and Y coordinates. In the "geometry" field I only see a very long combination of numbers and letters.

While looking for a way to adding those columns I found this short guide:

http://pvanb.wordpress.com/2012/02/21/extract-x-and-y-coordinates-from-a-spatial-table-in-spatialite/

Using the Plugin "QSpatialite" and modifying the following code passage I managed to create another Table in QGIS with two columns containg the x and y value.
I adapted this code:

CREATE TABLE 'Test' AS

SELECT PKUID, X(Geometry) as X, Y(Geometry) as Y

FROM Points

I did not even know where to enter the code until I installed the plugin.

Though, what I really want is to add two additional columns in the original table and to display the x/y values in those.

How do I do this?


texts: PostGIS compilation error on Windows - missing libpq

texts:
I have built PostgreSQL source using MSVC2015 compiler for 64bit support and installed on Windows. Now, I want to build PostGIS also from source (I am referencing from this documentation). I am using MinGW. First, I had problems with a missing PGXS makefile which I have solved by installing the PostgreSQL server development package (from Debian). Now, the configure file is complaining about missing libpq:

checking for PQserverVersion in -lpq... no
configure:error: could not find libpq

Upon closer inspection at the config.log file, this stands out:

error C:/PostgreSQL/9.6/lib/libpq.dll: file not recognized: File format not recognized

I have found this which suggests either an incompatibility between compilers or trying to combine a 64bit version with a 32bit one. Well, I made sure to compile in 64bit throughout - maybe I can somehow verify this?
Otherwise, if it is a problem with compilers, the PostgreSQL compilation was only working with MSVC and I cannot use this for PostGIS to my knowledge.
I am trying to make a contribution to PostGIS so there is no point using the binary installation. I understand it is not ideal to use Windows, however I would very much like the solution to work on Windows at the end of the day.


texts: Relation does not exist when configuring Mapnik with PostGIS table

texts:
I am trying to create a Mapnik layer based on a PostGIS table, and I keep getting an error message when trying to load the layer

Postgis Plugin: ERROR: relation "public.myTable" does not exist ...

. This is my first time working with PostGIS and I've gotten Mapnik working with other layer types, but I was wanting to experiment with the possibility of using PostGIS.

I created the table I'm trying to use with QGIS, and I can see the table inside of QGIS in the database manager and through psql.

This is what my layer definition looks like in Mapnik

  <Layer name="myLayer" srs="+init=epsg:4326">
    <StyleName>My Style</StyleName>
    <Datasource>
      <Parameter name="type">postgis</Parameter>
      <Parameter name="host">localhost</Parameter>
      <Parameter name="dbname">postgis</Parameter>
      <Parameter name="user">postgres</Parameter>
      <Parameter name="password">[password]</Parameter>
      <Parameter name="table">public.myTable</Parameter>
      <Parameter name="estimate_extent">true</Parameter>
    </Datasource>
  </Layer>

I've tried "public.myTable", "public"."myTable", and myTable. I've tried uppercase, lowercase, and actual case.

Why is Mapnik failing to load this table?


texts: Avoiding SQL injection when updating tables through CartoDB API

texts:
I'm trying to write a routine which synchronises data in our database with the CartoDB database, so we can display it on a map. It's simple things like names and contact details, which our users have the ability to change online. It's unlikely to be a problem but nevertheless it's prudent to worry about SQL injection.

Is there a recommended way to do this through the SQL API? I'm writing queries such as:

https://ouraccount.cartodb.com/api/v2/sql?q=update our_table set surname = 'Smith' where cartodb_id = 5&api_key=123456789

I'd rather not try to sanitize the 'Smith' part myself, and I can't find any documentation saying that the API supports parameterisation. What do people do here?


texts: Automated Shapefile/Feature Class source credit from metadata in ArcGIS layout

texts:
In ArcGIS (10.1), is there a way to automatically place the attribution (= source credits) of my data onto the map layout? The source credits are stored within the metadata (Item Description) and the shapefiles/feature classes are located within the folder structure.

The data within the folder structur is accessed by many users and not all of them know which attribution has to be used when they add the layer to their mxd.

Im thinking of the case of web mapping solutions, where one can specify that the source credits of visible layers are placed in the bottom left or right corner.

Any hints or ideas are welcome.


texts:
I created some layers and upload it to PostGIS and then GeoServer and then styled it using GeoExplorer. I have done all this in localhost. I used the iframe generated in GeoExplorer and used it in my webpage. So I can open it in my computer. But I don't understand how to keep this in remote server so that everyone can access it.
I used OpenGeoSuite for doing all this.
Is the data I uploaded in PostGIS and GeoServer just located in my computer or in remote server? Please help!!

Question Extension

This is the code that I am trying to run. But it's not working even in the same machine:

var map;

function init() {
   map = new OpenLayers.Map('map_element', {});

   var road_name = new OpenLayers.Layer.WMS(
            "IIRS:river_canal_project - Tiled",
"http://localhost:8080/geoserver/IIRS/wms",
            {
                LAYERS: 'IIRS:river_canal_project',
                STYLES: '',
                format: format,
                tiled: true,
                tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom
            },
            {
                buffer: 0,
                displayOutsideMaxExtent: true,
                isBaseLayer: false,
                yx : {'EPSG:4326' : true}
            }
        );

   Var industrial_park = new OpenLayers.Layer.WMS(
            "IIRS:future_industrial_project - Tiled",
"http://localhost:8080/geoserver/IIRS/wms",
            {
                LAYERS: 'IIRS:future_industrial_project',
                STYLES: '',
                format: format,
                tiled: true,
                tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom
            },
            {
                buffer: 0,
                displayOutsideMaxExtent: true,
                isBaseLayer: false,
                yx : {'EPSG:4326' : true}
            }
        );

   var land_use = new OpenLayers.Layer.WMS(
            "http://localhost:8080/geoserver/IIRS/wms",
            {
                LAYERS: 'IIRS:current_landuse_project',
                STYLES: '',
                format: format
            },
            {
               singleTile: true,
               ratio: 1,
               isBaseLayer: true,
               yx : {'EPSG:4326' : true}
            }
        );
   map.addLayers([land_use, industrial_park, road_name]);
   if(!map.getCenter()){
       map.zoomToMaxExtent();
   }
}

Along with the code, I inserted openlayers.js in script and only that.
I found myself some mistakes and now I that can load the map but..the images are not loading.
This is my edited code :

map = new OpenLayers.Map('map_element', {});
            // setup tiled layer
            var tiled = new OpenLayers.Layer.WMS(
                "earth:ne_10m_railroads - Tiled",
"http://localhost:8080/geoserver/earth/wms",
                {
                    LAYERS: 'earth:ne_10m_railroads',
                    STYLES: '',
                    format: OpenLayers.Format.WMSGetFeatureInfo,
                    tiled: true,
                    tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom
                },
                {
                    buffer: 0,
                    displayOutsideMaxExtent: true,
                    isBaseLayer: true,
                    yx : {'EPSG:900913' : false}
                }
            );

            map.addLayer(tiled);
var bounds = new OpenLayers.Bounds(
                -135.326858520508, -30.5615100860596,
                179.357788085938, 69.604377746582
);

    map.zoomToExtent(bounds);
}


texts: Joining over 255 fields to shapefile using ArcGIS Desktop?

texts:
I am trying to create a time series visualization of unemployment rates for every county in the continental US using  ArcMap 10.3.1.

I have data for every month going back to January 1990, but when I try to join it to my county shapefile in ArcMap, I get an error message because the spreadsheet has too many columns.

How can I join this spreadsheet to a shapefile such that I can symbolize the same shapefile by each field for which I have data?


texts: QGIS integer raster stats calculation error

texts:
I am comparing two raster temperature files generated from two different sources that use different algorithms. One of them is provided as a floating point (temp(C)) and the other is integer (temp(C)*10). I generated statistics for both using QGIS Raster Layer Statstics. Both files have extensive areas of no data values (-3.402...e+38 and -9999).

When I calculate statistics, QGIS provides an apparently erroneous mean value for the integer-scaled file:

  Valid cells: 153660
  No-data cells: 200134
  Minimum value: 141
  Maximum value: 270
  Sum: 36204112
  Mean value: 147
  Standard deviation: 84.0595027347

In particular note that 36204112/153660 = 235.61, not the reported 147.
This problem appears to affect the Standard deviation as well.

The same operation appears to work for corresponding the floating point file:

  Valid cells: 153773
  No-data cells: 200021
  Minimum value: 14.6999998093
  Maximum value: 26.0
  Sum: 3524999.0977
  Mean value: 22.923394209
  Standard deviation: 1.56633585615

In this case, the operation 3524999/153773 produces the same result as the output Mean (22.92...)

Indeed, when I convert the integer file to float (raster calculator: * 0.1), the stats calculation seems to work as expected. And when I convert the other file that was originally in float to integer, the stats calculation becomes erroneous.

In other words, there appears to be an error in stats calculation of integer raster files.

Is this a bug or am I missing something?

Perhaps the best way to avoid is to always use float.

I am running QGIS 2.18.2 on MacOSX10.10.5.


texts: QGIS upgrade from 1.8 to 2.12: Vector layers are not displayed

texts:
I just updated my very outdated version of QGIS. I was running 1.8 and thought the new version would work out a couple kinks that had developed. Well, now I can't even figure out how to make my vector data files appear after I add the layers.

The vector files are in the layers panel but nothing on the screen. Everything is projected and I try to zoom to layer and nothing! I'm just trying to add .csv files of camera trap locations and other files with gps points. Is there some simple fix like a setting I have to activate?


texts: How to get "inner" buffer from closed multiline?

texts:
I have a vector layer containing a closed boundary multiline, and I'd like to extract a certain buffer which would exclude the part that fall outside of the closed boundary.

How would I accomplish this?
I'm using QGIS 2.18.19.


texts: Can I get a list of public tables from CartoDB without authentication?

texts:
I've been running
"SELECT cdb_usertables FROM CDB_UserTables('public')"
over the CartoDB SQL API, but I'm getting 0 rows in return (although, notably, no authentication error).

It was my impression that callers no longer needed an API key to query public tables:
https://github.com/CartoDB/cartodb-postgresql/issues/36

Am I doing something wrong?

Related: Get list of all tables in CartoDB account with SQL statement?


texts:
ers,
At my prior position, we "often" had to capture Geology data from geo-referenced images.  Yes ... sometimes a necessary evil, but it has been a while since my hardcore digitizing days.  Regardless, we used this handy add-on for ArcMap called ET GeoTools.  We were able to simply digitize the polylines for the adjacent Geological categories, then run an analysis of over/undershoots.  Once we had that issue fixed, you could simply convert the polylines to polygons and build your polygon attribution in to the Database from there.  It was a beautiful thing!!!!  I didn't have to concern myself with snapping to every node on adjacent polygons, and this little add-on had a number of other neat little capabilities as well.
Now, I looked in to it, and it turns out they want $195.00 for the software.  I don't have any knowledge of programming languages - I should make that clear too before I move on.  Anyway,  I wondered 2 things, and perhaps a few of you may have a suggestion or two.  Does ArcMap 10 have this functionality (I am in 9.3 but soon to be using 10).  OR .... is there a free add-on that can do this?  Even if there's something out there that will export to .shp format, that will suffice.  I'll just take care of the rest of the work in good ol' ESRI following the conversion.
Anyone?  Thanks in advance


texts: GRASS GIS v.net.iso prompt not calculating correct travel time distances

texts:
When I run the v.net.iso command then export it into QGIS, the distance isochrones are not calculating my prompt distances (e.g. 1km, 2km, 5km).

1km distances (shown in green) are next to the 5km distances (purple), is close proximity of my starting location. Is this something to do with my road shapefile? I've tried another road file in another country and GRASS is calculating it correctly. Any ideas on how I can fix the problem?


texts: getImageData() returns all 0 values in data for an OpenLayers3 layer set as Transparent (Safari 8 on OSX Yosemite)

texts:
Please pardon my ignorance if this seems trivial, but I am very new to using openlayers3 for map rendering.

I am rendering a Tiled layer over the world map layer for the KML data from CGI mapserver 7.0.0. The KML layer tiles come up fine on display but when I click on the layer and call the forEachLayerAtPixel() [ from the onMouseClick()] it seems not to fire the callback associated with it as the data element in getImageData() returns all zeroes for R G B Aplha values.

When I make this layer as transparent=false, the data array gives me non zero values and the layer information comes up.

I am using Safari Browser (8) on Mac OSX for viewing the map and OL layers.

Here is how I declare the KML layer:

$rootScope.kmlLayerSource = new ol.source.TileWMS({
                                        url: someURL,
                                        params: {'LAYERS': 'kmltiles', format:'image/png', TRANSPARENT:true},
                                        crossOrigin:''
                                    });

$rootScope.kmlLayer = new ol.layer.Tile({ source: $rootScope.kmlLayerSource });
$rootScope.map.addLayer($rootScope.kmlLayer);

But when I change the above to layer def to TRANSPARENT: false, the data array in the getImageData() call returns positive integers.

Am I missing any settings here while creating my layer?


texts: Using NASA's GIBS WMTS servers with Mapbox mgl-vector-source

texts: I try to render https://gibs.earthdata.nasa.gov/wmts/epsg4326/best/wmts.cgi?TIME=${date}&layer=VIIRS_NOAA20_Thermal_Anomalies_375m_All&tilematrixset=500m&Service=WMTS&Request=GetTile&Version=1.0.0&FORMAT=application%2Fvnd.mapbox-vector-tile&TileMatrix={z}&TileCol={x}&TileRow={y} by Mapbox vector source mgl-vector-source. But when I implement it I get some shift in thermal anomalies data visualisation. I cant get along with it.


texts:
I just got a new MacBook Pro and installed QGIS 3.8 yesterday. I had several GDBs working on Ubuntu and QGIS 3.4 but that same GDB will not import in 3.8 on the MacBook Pro. They are vector layer databases and not raster files.

I've been converting the GDB to GeoPackages (on the Unbuntu/QGIS 3.4) which are much faster and flexible. But I can't get much to work on the MacBook in the way of the GDB files.

I used the new package installer from the QGIS download page and noticed it streamlined the installation from previous MacOS installs but I'm wondering if maybe the 3.8 package was ready for distribution yet.


texts: Exporting set of layers into multiple PDF (or Adobe Illustrator) files, one at a time?

texts:
I am completely NEW to PYTHON(only took a 3-hr entry course in ESRI tutorial. I have a dataset contains files named as "1_A, 1_B, 1_C....12_A,12_B, 12_c,....200_A, 200_B, 200_C..etc." I need to export each feature layer, one at a time , to a PDF or AI file through a COMPLETELY AUTOMATED process.

I tried to use Addlayer and Removelayer tool in Model Builder but they all run for one single feature layer once at a time. I coudl not connect the add and remove process to each other, nor could I connect the process through Iterator to the script tool anyway.

This means although I use model builder or write python script, I have to change the layers (or the name of layers) hundreds of times and save every layer into a new MXD file and export then to pdfs.

I was thinking about turning on one of the layers at a time, export a .ai file with that layer visible, then turn that layer off and move on to the next layer in the list. I'm wondering if you have any suggestion about building the model or the scripts?


texts: Difference between files of extension .asc, .asc.aux.xml, .aux and .rrd?

texts:
I have received ASCII data from a client. These ASCII files have been created from LiDAR data.

For each tile there are 4 different file types:

Filename.asc
Filename.asc.aux.xml
Filename.aux
Filename.rrd

What is the purpose of these files (asc.aux.xml, aux, .rrd)?


texts: Understanding Openlayers MapServer configuration?

texts:
I need to implement functionality like this http://www.mpbhuabhilekh.nic.in/bhunakshaweb/.
I searched on net and found it can be done using Openlayers. I need simple and clear steps on how to to configure Openlayers, Mapserver or WMS.

I am working on windows 7 with xampp. I am using Openlayers2 with Mapserver in which we need map file to show the map. My map file name is sample.map and its location is C:/ms4w/apps/mapserv-demo/sample.map

Content of my sample.map file:

#Start of map file
MAP
    NAME "sample"
    STATUS ON
    SIZE 600 600
    EXTENT 388107.634400379 5203120.88405952 500896.339019834 5310243.30613897
    IMAGECOLOR 255 255 255

    #Start of web interface definition
    WEB
        #CONFIG "MS_ERRORFILE" "/ms4w/tmp/ms_error.txt"

        MINSCALE 1000
        MAXSCALE 1550000
        IMAGEPATH "C:/ms4w/tmp/ms_tmp/"
        IMAGEURL "/ms_tmp/"

        METADATA
            WMS_TITLE "Sample Demo"
            WMS_ABSTRACT "This is a sample demo"
            WMS_ACCESSCONSTRAINTS "none"

            #change this value to match your setup
            WMS_ONLINERESOURCE "http://localhost/cgi-bin/mapserv.exe?"

            WMS_SRS "EPSG:30134"
        END
    END #WEB

    #Start of layer definitions
    LAYER
        NAME 'global'
        TYPE RASTER
        STATUS DEFAULT

        CLASS
            NAME 'Digital Raster Graphic'
            KEYIMAGE graphics/drgs_keyimage.png
        END

        METADATA
            WMS_TITLE "Sample Demo"
            WMS_ABSTRACT "This is sample demo"
            WMS_ACCESSCONSTRAINTS "none"

            WMS_SRS "EPSG:30134"
        END
    END #LAYER

END #MAP

I am running only my map file to test if its working or not like this:

http://localhost/cgi-bin/mapserv.exe?mode=map&map=C:/ms4w/apps/mapserv-demo/sample.map

But it's not working, it's showing the blank page.


texts: Is there anyone who has worked with FARSITE?  I can build some of the files needed, such as slope and aspect, but not sure where to start when it comes to building files for 'canopy', 'fuels', etc.  I am currently using QGIS 1.8.0.  The tutorial provided with FARSITE is excellent, and I worked through that easily.  It's in creating the files for my area of interest that I'm encountering problems.  These would have to be created for South African conditions.


texts: GeoExt Popup always loading after first click

texts:
I'm trying to put popups on my map to show the feature's info. I have done that, however after the first click the mouse pointer starts loading and never stops.

var control = new OpenLayers.Control.WMSGetFeatureInfo({
              url: "/geoserver/prod/wms",
              layers: [restaurantes,ruas,linhas_agua]
});

control.events.on({
        getfeatureinfo: function(event) {
        if (popup)  { popup.destroy(); }
        popup = new GeoExt.Popup({
        title: "Informação",
        map: mapPanel,
        layers: [restaurantes,ruas,linhas_agua],
        width: 250,
        autoScroll: true,
        collapsible: true,
        bodyStyle: {padding: 5},
        html: event.text,
        });
        popup.show();
        }

});

mapPanel.map.addControl(control);
control.activate();


texts: Deriving Z-score values for given field using ArcGIS for Desktop?

texts:
I have a shapefile with an attribute field with a range of numerical values.

How do I derive z-score values for this field using ArcGIS for Desktop?


texts:
I would like to create a map of point density of vector features /archealogical sites/ per hectar in QGIS /something like in an attached figure/. I doubt if the kernel  density is a proper tool for this objective. Are there any other more appropriate procedures?
Or in another words what do numbers in KDE output actually represent. Is it so that:

"Raw values" in KDE = number of points per cell ?,
"Scaled by kernal size" = number of points per map unit ?
Is the term "kernal size" a synonym for the chosen radius or does it refer to the cell area " ?


texts: Is clip and ship possible with ArcGIS Server map service?

texts:
Is clip and ship possible with a map service?

Following the tutorial, it appears you can publish a clip and ship GP service, but it will only work with Arc Map layers (SDE, FGDB layer, etc...) and not a Map/Feature Service. Is that correct?

We'd like to be able  to publish a clip and ship that allowes users to extract data from a web map.

If this is possible, are there any configuration examples?

Not finding any on Esri web site.


texts: Placing label designs in legend of ArcGIS Desktop?

texts:
I'm working in ArcGIS 9.3.1 and apart from "regular" slayers with symbology, I have a few layers which I do not show their symbol on the map, but only their label (for example: cities layer in which I've made the point symbols transparent and their labels printed in blue).

Is there any way I can show the labels of the layers on the legend instead of their symbols (which is the default)??


texts: ArcGIS SDE next rowid value doesn't increase at a constant rate

texts:
Using ArcGIS 10.5.1 and Postgresql 9.5 database, I created a point feature class named Locations. After I inserted 8 points in my Locations table, the objectid value of last row is 8.

An hour later, when I inserted a new point to the Locations table the object id value created as 408. But it should be 9.

Several hours later when I inserted a new Location, the objectid created 1220.

Why don't objectid values increase at a constant rate?


texts: Extract the elements not in common between two different point sources

texts: I've downloaded accommodations from two different sources, now I have the list of elements from them in two separated shape files and also the list of elements in common between these two sources. How can I extract the elements from the attribute table that are only in the first source and not in the second and viceversa? I need something like a difference between two point vectors.


texts: Execute GRASS command within QGIS python console

texts:
We try to calculate catchments with the r.water.outlet function from python console. To do so we use the following code:

import processing
rlayer = QgsRasterLayer("/my/path/DrainageDirection.tif", "DrainageDirection")
processing.runalg("grass:r.water.outlet", rlayer , 2621378.611131 ,1262320.888212, "2589700, 2640500, 1242200, 1270000", 73, "/my/path/zzz.tif")

The result in the zzz.tif file is only one pixel with no value. This result differs when running the r.water.outlet is called from the GRASS Tools GUI.

So, three questions:

1) Do we handle the input raster (rlayer) correct?
2) Do we handle the output raster correct?
3) Would it be possible to read a GRASS Mapset raster directly from python console?


texts: Opening'set location' box of geopoint within repeats in Survey123?

texts: In my survey, each page contains a repeat for a set of questions and geopoint. The geopoint set location box is open for the repeat on the first page, but it is closed for the rest of the pages. Is there a way I could make sure all the set location boxes are automatically open instead of having to click on them?


texts: How to set the order of layers in mapservers legend graphic?

texts: Is it possible to set the order of the displayed elements (Layers) in the legend graphic recieved by WMS GetLegendGraphicREQUEST or will MAPSERVER always return in the order of layers defined in the map-file?


texts: What is correct entry to Input Location (GPFeatureRecordSetLayer)*?

texts:
I use:

http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Network/ESRI_DriveTime_US/GPServer/CreateDriveTimePolygons/execute

What is correct entry to Input Location (GPFeatureRecordSetLayer) ?

I try:

{"type": "Point","coordinates": [-105.01621,39.57422]}

but I get : "Invalid value for parameter 'Input_Location'"


texts:
After loading a point layer out of a text file (plug in add delimited text layer) saved in ISO 8859-15 (Latin 9) out of Access, I did save the layer under Qgis 1.8 with the ISO 8859-15 encoding.

When going into the properties of the layer, it stays in UTF-8 (layer properties - general - Provider specific options - Encoding). No matter what encoding I set there, after closing and reopening the properties window it still stays in UTF-8.

Encoding seems to be still correct within the desktop Qgis, but when opening the project in the Qgis-web-client, not all special characters are correctly displayed.

Does anybody know a remedy?


texts: How To Convert A Set of ARC/Info ASCII .e00, .e01, .e02, etc Files Into .shp Files

texts:
I have a set files that I downloaded with extensions (.e00, .e01, .e02, etc).  I have read that these are files used by ARC Info mainly, but I do not have ARC and nor am I really that interested in downloading it.  The funnier thing is that I am able to convert a .e00 file in QGIS to a .shp set, but if the file has been broken into many .e0n files, QGIS does not recognize the others and leaves me with only have my shapefile.  I have tried taking the other files and copying and pasting their contents to the end of the .e00 file, but this resulted in a very erroneous shape in QGIS.  I have been searching for some while now if there is a way to do this with ogr2ogr, but haven't found anything for a set of files like this.  I also downloaded e002shp, but had no idea on how to go about using it and didn't find much online for that one.  If anyone can tell me any command-line tools I could use (preferably ogr2ogr) or a way to do this in QGIS, I'd gladly appreciate it :)

Actually I'd appreciate any tips, but that is the software which I am more familiar with ;)

Thanks,

elshae


texts: How to create several header files to open AVHRR NDVI3g (GIMMS) files?

texts:
I have downloaded GIMMS NDVI3g data from http://ecocast.arc.nasa.gov/data/pub/gimms/3g/

I now need to open the files in arcgis. I have the follwing GRID PARAMETERS:

grid-name: Geographic Lat/Lon
    pixel-size: 1/12=0.0833 degrees

size-x: 4320
    size-y: 2160

upper-left-lat: 90.0-1/24
    upper-left-lon: -180.0+1/24
    lower-right-lat: -90.0+1/24
    lower-right-lon: 180.0-1/24

I want to ask if someone has worked with these files yet since the release is recent and how you opened the files? And if a header file is needed how do I create multiple header files in one step? (there are 700 files covering 3 decades)...

Best regards,
Johan


texts: Add character to end of layer names of rasterStack in R


texts: Decrease resolution of vector shapes (polygon) i.e. reduce size and number of points in outer ring

texts:
I'm trying to get a kmz export down to less than 3 megabytes and am wanting to do so by reducing the number of points in my polygons' outer rings.

I've seen this concept in action when it comes to buffer resolutions. If I were to create a buffer around a point with a resolution of 1 it would be a square, if it had a res of 2 it would be an octagon, and so on.

Is there something similar I could do to simple polygons whose outside rings are too detailed for my purposes?

Thanks


texts: What's wrong in using latitude and longitude as X and Y on a Cartesian space?

texts: What can go wrong and what kind of error can build up if I use geodetic coordinates as if they were cartesian and place my objects into my world with no conversion formulas?


texts: Python code to determine distance to specific polygon boundary

texts:
I am trying to automate a process to assign location codes to parcels.  The method is to measure from the centroid of the parcel polygon to the eastern edge of the section polygon, then measure from the same centroid to the northern edge.  The order is important.

I was looking at the PointDistance_analysis tool/command, but since the points in question lie within the polygon edges that I am measuring to, it won't return a viable number.  In addition, I have been unable to determine how to tell it to measure to a specific polygon side.

Doing this manually is simple and takes about 2 minutes per parcel.  However, I have almost 100,000 parcels to do, so an automated process would be preferred.

Any suggestions?


texts: How do I draw the shortest imaginary boundary around a point dataset?

texts:
I would like to produce a range map for my species (smooth snakes in the south of England) using the IUCN-defined 'extent of occurrence' (EOO) as my metric. EOO is defined as 'the area contained within the shortest continuous imaginary boundary which can be drawn to encompass all the known, inferred or projected sites of present occurrence of a taxon, excluding cases of vagrancy'- see below where A is the location of all records and B is their extent of occurrence:

I would like to do this in ArcMap 10.3.1 if possible. I first tried it using the 'Minimum Bounding Geometry' tool (Data Management) but this just produced a square polygon encompassing all records rather than the 'shortest continuous imaginary boundary' that I need:

After some research I concluded what I probably needed to do was to create a concave hull, for which I found Richard Fairhurst's 'ConcaveHullByCase' tool (I was not able to post a link as I had already included two images - not sure why this restriction exists, surely adding links and screenshots helps to illustrate/clarify a point and it for everyone's benefit?)

The result I got from this is not quite what I'm after either, however. Firstly, it does not encompass all points. Secondly, it draws a line between all points, rather than the shortest continuous line as in example B above. It also creates a 'polygon within a polygon' (white rectangle within turquoise rectangle, see below) that I don't fully understand - it does not appear in the output's attribute table but seems to be a separate feature:

Is there a solution to my problem besides manually drawing these polygons myself?? Apologies for the rather descriptive/untechnical way I'm phrasing this, I'm still a novice and had to look up the difference between convex and concave before posting this, which should give you an idea of how much I still have to learn.