from:"Even Rouault"

Re: [Qgis-user] QGIS w/ Oracle 19c

2024-06-24 Per discussione Even Rouault via QGIS-User


Will,

I presume Oracle must have changed their WKTEXT definition.  The issue 
here was the naming "Albers_Conical_Equal_Area" that must be a Oraclism 
and a not-yet supported alias for this projection name. I've queued 
support for it in PROJ (https://github.com/OSGeo/PROJ/pull/4181), for 
the next PROJ release


Even

Le 24/06/2024 à 23:15, Burt, William WLRS:EX via QGIS-User a écrit :


Hello QGIS users! I am wondering if anyone has any experience with 
this issue I am having working with Oracle DB and QGIS.


When new layers are added from Oracle via the QGIS Browser panel I get 
an error message in QGIS:


/No transform is available between NAD83 / BC Albers and EPSG:4326 - 
WGS 84./


/Unknown error (code 4096) /

Some details:

 1. I am certain the root of the issue is not QGIS as I get the error
message on both current and older versions of QGIS that were
working fine a few months ago.
 2. Our enterprise Oracle db was recently upgraded to 19c. This may or
may not be the cause but seems re
 3. The error is only thrown for data with BC Albers coordinate
systems. WGS 84 datasets work fine
 4. Datasets where the URI is defined by .qqz or data added via our
SLYR plugin via .lyr .lyrx  all work fine.
 5. Dataset can be made to work if the assigned CRS is changed to
EPSG:3005 – NAD83/BC Albers from the auto detected crs: NAD83 / BC
Albers

This leads me to think that the CRS metadata is being delivered to 
QGIS differently than before. Below is the WKTEXT for SRID 3005 from 
the Oracle medadata table MDSYS.CS_SRS (I can’t tell if this has 
changed from the older Oracle and I don’t see why it would):


PROJCS["NAD83 / BC 
Albers",GEOGCS["NAD83",DATUM["North_American_Datum_1983",SPHEROID["GRS 
1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6269"]],PRIMEM["Greenwich",0],UNIT["Decimal 
Degree",0.0174532925199433]],PROJECTION["Albers_Conical_Equal_Area"],PARAMETER["Latitude_Of_Origin",45],PARAMETER["Central_Meridian",-126],PARAMETER["Standard_Parallel_1",50],PARAMETER["Standard_Parallel_2",58.5],PARAMETER["False_Easting",100],PARAMETER["False_Northing",0],UNIT["Meter",1],AXIS["Easting",EAST],AXIS["Northing",NORTH]]


Thank you,

Will


___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info:https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.
___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Vermont Coordinate System

2024-05-15 Per discussione Even Rouault via QGIS-User

I couldn't spot anything that directly maps to "Vermont Coordinate 
System", but https://spatialreference.org/ref/?search=vermont lists a 
few serious candiates. Choose your disease :-)



Le 16/05/2024 à 00:43, Chris Olson via QGIS-User a écrit :

Is it possible to have the Vermont Coordinate System (VCS) used by QGIS?

___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.

___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Search through folders and combine SHP files

2024-05-12 Per discussione Even Rouault via QGIS-User


Hi,

The GDAL ogrtindex might be of good help to create an index of various 
vector sources: https://gdal.org/programs/ogrtindex.html


And you can possibly "gdalmanage identify -r some_path" to find all GDAL 
recognizable datasets: https://gdal.org/programs/gdalmanage.html


Even

Le 12/05/2024 à 23:41, Tony Shepherd (FarmMaps NZ) via QGIS-User a écrit :


Hi Nigel

Thank you, that works a treat, although it does fail if I try to merge 
a line feature at the same time as a polygon.  That is not a big deal 
and can get around that.


Your comment, “Then parse these through to the merge vector layer”.  I 
can’t see a way of parsing the filelist.txt to the tool. I am not 
experienced in using python.


Cheers

TIby

*From:*Nigel Berjak - S3 Technologies 
*Sent:* Friday, May 10, 2024 6:03 PM
*To:* Nyall Dawson 
*Cc:* Tony Shepherd (FarmMaps NZ) ; 
qgis-user@lists.osgeo.org

*Subject:* Re: [Qgis-user] Search through folders and combine SHP files

Hi

For your listing, although it would not provide the file type i.e. 
point, line, polygon etc, you could simply create a DOS batch command 
(dir/s *.shp > shapefilelist.txt) to locate all the file and generate 
a list of their locations. Then parse these through to the merge 
vector layer. If you are merging them all into a Geopackage, then file 
type would not be relevant and you could then sort them once they are 
in a single location. I don't think the Geopackage would have a 
limitation on the number of files and records in it.


Perhaps this can assist.

---

Regards,

Nigel Berjak
S3 Technologies
Geographic Information Systems & Large Format Printing specialists
T: +27 33 3423681
F: +27 86 6721242
E: ni...@s3.co.za 
Website: http://www.S3.co.za 

Please consider the environment before printing this email.

On 2024-05-10 04:06, Nyall Dawson via QGIS-User wrote:

On Fri, 10 May 2024 at 09:06, Tony Shepherd (FarmMaps NZ) via
QGIS-User mailto:qgis-user@lists.osgeo.org>> wrote:

Hi

I have loads, literally thousands, of SHP, GPKG, TAB files
located in loads of different folders and sub-folders.

Is there a plugin that will trawl through a set of folders,
copy say just the polygons from each file into a new file, and
populate a field with say the filename and folder the
polygon(s) came from?

Essentially combining a heap of files into one with some
metadata about where the data came from?

Why?  I have files from clients dating back about 20 years.
Clients sell properties.  Properties get subdivided.  New
clients come and see us, and often we don't realise we have
GIS info on file for that property from many years ago,
usually under a different name.  It would be handy to see a
polygon on the screen and easily go back to some old GIS data.

Maybe I am missing an easy tool to do this, but I can't see
anything obvious.

You could use the "Merge Vector Layers" tool from the Processing
toolbox, that does everything you want. The trickiest bit would be
getting a master list of ALL your shapefiles/etc you want to
combine and passing this to the tool. It's not going to be easy
from the GUI itself, as that only allows you to add files from a
single folder at once.

I'd run the tool on the files from a single folder, then checkout
the processing history dialog and copy the python command which
corresponds to what you just did. You could then adapt this
command to pass a complete list of all the files you want to
combine. (It's quite straightforward in Python to build a list of
files which recursively match a file pattern -- see eg
https://stackoverflow.com/a/2186565
)

Hope that helps!

Nyall

Cheers

Tony


--

*Tony Shepherd*

/  GeoSpatial Manager @ FarmMaps NZ & Photographer @ Shepherd
Photos/

/*Phone*/*–* 027 435 6193  | /*Website*/*–
*shepherdphotos.co.nz 

/*Email*/*– Maps*t...@farmmaps.nz |
/*Email*/*– Photos* t...@shepherdphotos.co.nz


/*Facebook*/TonyShepherdPhotos


FM_logo1_small (Custom)

___
QGIS-User mailing list
QGIS-User@lists.osgeo.org 
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user

Unsubscribe:
https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] building gdal for mrsid

2024-01-21 Per discussione Even Rouault via QGIS-User

No, this won't be a .deb, and using the past .deb for mrsid won't work
since the switch to CMake

General instructions to build GDAL from source are at
https://gdal.org/development/building_from_source.html

To install it in a location where it won't mess with with your packaged
GDAL, do for example

cmake .. -DCMAKE_INSTALL_PREFIX=/opt/gdal-3.8.3 (others options here
such as -DMRSID_ROOT=... -DECW_ROOT=... )

(or any other installation directory of your choice)

Even

Le 21/01/2024 à 16:08, Don Harter a écrit :

OK. Thanks. Those scripts don't use .deb packages so nothing would be
registered with dpkg. I have had problems before with unregistered
packages lying around. I just fixed such an issue with GCC. I don't
even see a package that is available for libecwj2. In what directory
is the build script run in? I already have gdal installed. Will it
overwrite the existing ubuntu locations for the gdal packages? If I
were to purge those packages then other packages that rely on them
would have unfulfilled dependencies. I see that script is for an
earlier version of ubuntu. I am running : Ubuntu 22.04.3 LTS. And my
gcc (gcc (GCC) 12.3.1 20230626) is newer so I have downloaded the
other/newer mrsid file.

So if I download the source for the .deb packages and replace the
scripts in them with those 2 that you located for me, will that work?

On 1/21/24 8:51 AM, Even Rouault wrote:

Hi,

cf
https://github.com/OSGeo/gdal/blob/master/.github/workflows/ubuntu_20.04/Dockerfile.ci#L144
and
https://github.com/OSGeo/gdal/blob/master/.github/workflows/ubuntu_20.04/build.sh#L15
for a build recipee to build GDAL >= 3.5 with MrSID support

Even

Le 21/01/2024 à 03:10, Don Harter via QGIS-User a écrit :

I have been following this procedure to build gdal with mrsid
capabilities. I had done this a while back previously but the
different dependencies between packages cause a mess.

https://trac.osgeo.org/ubuntugis/wiki/TutorialMrSid

I tried to run this command:

sudo apt-get install libgdal-mrsid-src

But that package is not available in the unbuntugis repository.

___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info:https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user

--
http://www.spatialys.com
My software is free, but my time generally not.

--
http://www.spatialys.com
My software is free, but my time generally not.
___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] building gdal for mrsid

2024-01-21 Per discussione Even Rouault via QGIS-User


Hi,

cf 
https://github.com/OSGeo/gdal/blob/master/.github/workflows/ubuntu_20.04/Dockerfile.ci#L144 
and 
https://github.com/OSGeo/gdal/blob/master/.github/workflows/ubuntu_20.04/build.sh#L15 
for a build recipee to build GDAL >= 3.5 with MrSID support


Even

Le 21/01/2024 à 03:10, Don Harter via QGIS-User a écrit :


I have been following this procedure to build gdal with mrsid 
capabilities.  I had done this a while back previously but the 
different dependencies between packages cause a mess.


https://trac.osgeo.org/ubuntugis/wiki/TutorialMrSid

I tried to run this command:

sudo apt-get install libgdal-mrsid-src

But that package is not available in the unbuntugis repository.

___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info:https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.
___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] NTv2-transformation: different output values in QGIS and GDAL

2024-01-16 Per discussione Even Rouault via QGIS-User


Le 16/01/2024 à 11:59, andy via QGIS-User a écrit :

Hi,
I'm doing NTv2-transformation, using this gsb file:
https://web.archive.org/web/20201018143532/https://provincia.agrigento.sitr.it/scaricabili/NTv2/ItalyED50toWGS84_NTV2_GN.gsb

I have copied it in "C:\OSGeo4W\share\proj" and I have defined this 
custom projection:
+proj=utm +zone=33 +ellps=intl +nadgrids=ItalyED50toWGS84_NTV2_GN.gsb 
+units=m +no_defs +type=crs


If I test it in QGIS, in the custom CRS window, using lon=13 and 
lat=38, I have in output:

324388.7743
4207772.1972


Did you try with the full path 
+nadgrids=C:\OSGeo4W\share\proj\ItalyED50toWGS84_NTV2_GN.gsb ?


If I do that (well on Linux: "+proj=utm +zone=33 +ellps=intl 
+nadgrids=/home/even/qgis/qgis/build_ubuntu2204/ItalyED50toWGS84_NTV2_GN.gsb 
+units=m +no_defs"), I do get the correct result 324454.6136 
4207894.9261 in QGIS (master)


But, what is the history of this "ghost" grid that requires fetching it 
from web.archive.org ? Has it had some official status, and could it be 
provided under an open data license ? If so, it would be great if 
Italian authorities could register to IOGP/EPSG this grid for ED50 -> 
WGS 84 for Italy... Then we could ship it in PROJ-data.


But I remember complicated discussion at FOSS4G Firenze with the guys 
from the geodetic authority (https://www.igmi.org) Not sure if that 
has lead to anywhere.


Even

--
http://www.spatialys.com
My software is free, but my time generally not.
___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Does QGIS require units "m" for the projection coordinate variables in a netCDF/CF file?

2023-09-16 Per discussione Even Rouault via QGIS-User

Thomas,

I am not from the GIS world and not even a QGIS user, so bear with me
if my answers do not use the expected vocabulary.
It takes years to get fully familar with GIS oddities, and netCDF is an
area where even a life long of experience will not be enough to fully
overcome the "creativity" of the netCDF community...

From my point of view, my netCDF files describe their geolocation in a
CF-compliant manner, with a grid_mapping CRS variable, and the x and y
projection coordinate variables. But QGIS does not recognize the CRS,
seemingly because I have my x and y as km (and I do specify
:units="km", I am not hiding this).

Since I am also a programmer, I tried to find where in the QGIS code
the netCDF/CF geolocation is read and decoded, to see if there was a
strict test on :units="m". But I failed to find this in the code.

This is actually not really done in QGIS, but in PROJ (more exactly
QGIS's QgsProjUtils::identifyCrs() delegates to PROJ's proj_identify()
in
https://github.com/qgis/QGIS/blob/8723b82ec19a3bec7f39b46128c798c7a2ee4230/src/core/proj/qgsprojutils.cpp#L271
to try to correlate the CRS definition from the netCDF file to one in
the database).

As you likely open it as a raster (and not a mesh as Richard
mentionned), the GDAL netCDF driver is also involved since it is it that
reconstructs a CRS object from the attributes of the netCDF CF conventions.

From what I can tell, both PROJ and GDAL do the job "as expected". It
is just that the way this CRS is encoded is not going to find any match
with a known CRS. To get what you want, the GDAL netCDF driver should
both modify the CRS definition to change the km unit to m, and alter the
geotransform matrix (which gives the coordinate of the top left corner
and pixel size) to multiply their value by 1000. This could be done, but
should it be done...? I'm not so sure. The issue is more than the
"philosophy" of netCDF data producers is sometimes at odds with the
usual CRS definitions.

That said, for that very particular case, there's a proj4_string
attribute in Lambert_Azimuthal_Equal_Area variable, with value
"+proj=laea +lat_0=90 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m
+no_defs +type=crs", which indicates the intent to have a metre unit.
There's also a global attribute geospatial_bounds_crs = "EPSG:6931"
which further correlate this. Hence this enhancement to the GDAL netCDF
driver to renormalize to metre: https://github.com/OSGeo/gdal/pull/8407

Even

--
http://www.spatialys.com
My software is free, but my time generally not.

___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] OAPIF does not honour the view extent

2023-09-15 Per discussione Even Rouault via QGIS-User


Jukka,

Fixed per https://github.com/qgis/QGIS/pull/54617 . The infinite items 
request came from when requesting the feature count.


The server returns sub-optimal content for QGIS usage:

- the collection metadata does not advertize a extent.spatial.bbox

- the /items answers don't include a "|numberMatched|" property that 
would help to get the feature count (QGIS emits a /items?limit=1 hoping 
that it returns numberMatched)


Even

Le 12/09/2023 à 17:01, Rahkonen Jukka via QGIS-User a écrit :


Hi,

I hoped that the OAPIF driver would honor the check box “Only request 
features overlapping the view extent” but at least in QGIS 3.32 that 
does not happen.


How to test:
- Set project into EPSG:4326

- Open the Debugging/development tools

- Zoom for example to lat-lon 63.189111,23.538514 and set scale into 
1:5000


- Add collection “building” from an open OAPIF service 
https://beta-paikkatieto.maanmittauslaitos.fi/inspire-buildings/features/v1/. 
Use feature paging with page size 1000.


- Remove the layer from the project because otherwise QGIS continues 
downloading a few million building polygons


- Check the http requests

The http history shows that after a bunch of probing OAPIF requests 
QGIS sends plain = request without = and thus asks the 
server to send the whole collection. What QGIS does right is that it 
follows the next-links and continues downloading.


https://beta-paikkatieto.maanmittauslaitos.fi/inspire-buildings/features/v1/collections/building/items?limit=1000=http://www.opengis.net/def/crs/EPSG/0/4326 



I wonder if the view extent limit should work also for OAPIF or is a 
not-yet implemented feature?


-Jukka Rahkonen-


___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info:https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.
___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] OGC API Features Part 4 - Editing with QGIS Desktop

2023-07-11 Per discussione Even Rouault via QGIS-User


Jakob,

yes this has been added in 3.32 per https://github.com/qgis/QGIS/pull/52504

Even

Le 11/07/2023 à 15:51, Jakob Miksch via QGIS-User a écrit :


Hello everyone,

I'm writing to ask about editing vector data with OGC API Features 
Part 4 (OGC API - Features ). I know 
that QGIS Desktop can read OGC API Features (see 22. Working with OGC 
/ ISO protocols — QGIS Documentation documentation 
).


But can it also edit them, like WFS-T does?

I couldn't find any information on this topic.

Any help would be appreciated.

Thanks, and best regards,

Jakob





Jakob Miksch
Consultant

+49 170 9036956
siticom.online



	siticom GmbH • Werner-von-Siemens-Straße 2-6 • 76646 Bruchsal • 
Deutschland

Geschäftsführer: Jürgen Hatzipantelis, Stefan Höltken
Sitz der Gesellschaft: Weiterstadt • Amtsgericht Darmstadt • HRB 9276
Umsatzsteuer-ID Nr.: DE 197 925 589



___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info:https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.
___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Export to GML in 3.28?

2023-03-22 Per discussione Even Rouault via QGIS-User

Le 23/03/2023 à 01:25, Eric Ortmayr a écrit :

Hi Evan,

Thank you for the information, the document was very helpful.

Since the default changed to GML3.2 do you think the FORMAT dropdown
should include a GML2 option? I am not sure if this was missed or if
this was done by design.

Good point ==> https://github.com/qgis/QGIS/pull/52350

Regards,

Eric

*From:*Even Rouault
*Sent:* March 21, 2023 5:24 PM
*To:* Eric Ortmayr ; qgis-user@lists.osgeo.org
*Subject:* Re: [Qgis-user] Export to GML in 3.28?

Le 22/03/2023 à 01:16, Eric Ortmayr via QGIS-User a écrit :

Hello,

Exporting to GML using version 3.16 results in a different
coordinate format than in version 3.28. In 3.16 the X and Y values
are separated by a comma (GML2?), but in 3.28 the X and Y values
are separated by a space (GML3?) . Both exports were created using
the FORMAT under the Datasource Options section in the
Save Vector Layer as… dialog.

Are the 3.16 and 3.28 FORMATS different by design? I
checked the changelogs and could not find anything obvious.

Yes, most file format export (and import) is done by GDAL. GDAL 3.4.0
changed to GML 3.2 as the default variant of GML. Cf
https://gdal.org/drivers/vector/gml.html#dataset-creation-options
<https://can01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgdal.org%2Fdrivers%2Fvector%2Fgml.html%23dataset-creation-options=05%7C01%7Ceortmayr%40crgl.ca%7Cc5533e1316c64ed7365c08db2a6baeba%7C2aa156d4b33c41688f39045b19f30a21%7C0%7C0%7C638150414206590656%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C=iLfZK7gmSVxOrNQDEwpRJqwiiDIiYL5GQkMINhJ6FY0%3D=0>

--
http://www.spatialys.com
<https://can01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.spatialys.com%2F=05%7C01%7Ceortmayr%40crgl.ca%7Cc5533e1316c64ed7365c08db2a6baeba%7C2aa156d4b33c41688f39045b19f30a21%7C0%7C0%7C638150414206590656%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C=CEO3ReB4CP1RSXlkBHPmECVoCbFwFN2Lo8M7ML6IxOk%3D=0>
My software is free, but my time generally not.

Re: [Qgis-user] Export to GML in 3.28?

2023-03-21 Per discussione Even Rouault via QGIS-User



Le 22/03/2023 à 01:16, Eric Ortmayr via QGIS-User a écrit :


Hello,

Exporting to GML using version 3.16 results in a different coordinate 
format than in version 3.28. In 3.16 the X and Y values are separated 
by a comma (GML2?), but in 3.28 the X and Y values are separated by a 
space (GML3?) . Both exports were created using the  FORMAT 
under the Datasource Options section in the Save Vector Layer as… dialog.


Are the 3.16 and 3.28  FORMATS different by design? I checked 
the changelogs and could not find anything obvious.


Yes, most file format export (and import) is done by GDAL.  GDAL 3.4.0 
changed to GML 3.2 as the default variant of GML. Cf 
https://gdal.org/drivers/vector/gml.html#dataset-creation-options



--
http://www.spatialys.com
My software is free, but my time generally not.
___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] File Geodatabase raster datasets readable in GDAL 3.7.0 (was Re: Unable to open some LISTdata open data ESRI File Geodatabases)

2023-03-09 Per discussione Even Rouault via QGIS-User

Hi,

Just to inform you that GDAL 3.7.0 will include a driver to read
directly File Geodatabase raster datasets:
https://github.com/OSGeo/gdal/pull/7389 . I've successfully tested it
against the below datasets, among others.

They will consequently be directly be usable by QGIS once QGIS is built
against this GDAL version (whose release will be in May)

Even

Le 16/12/2022 à 00:36, Nyall Dawson a écrit :

On Fri, 16 Dec 2022 at 09:18, Even Rouault via QGIS-User
wrote:

Hi Peter,

yes I confirm this dataset is a raster gdb. I can sort-of open it with
my QGIS + GDAL master build using the OpenFileGDB driver, but you'll see
the raster tables in a raw unusable way. For example layers like
fras_blk_NCH_ES_WATER_LOGGING_HAZARD which contains the tiles as
compressed data in a BLOB column.

You could try the https://github.com/r-barnes/ArcRasterRescue utility to
convert into GeoTIFF

FWIW, I've had more luck with the fork at
https://github.com/justinfisk/ArcRasterRescue

Nyall

You could also kindly or not kindly suggest your administration to not
use proprietary formats for open data whereas thereare so many open
alternatives, like GeoTIFF...

Small "advertizing": for those interesting, there's a funding
opportunity to develop a GDAL open source raster driver to be able to
read filegdb rasters that QGIS could use to open directly such files.
Contact me off-list if interested.

Even

Le 15/12/2022 à 23:49, JOLLY, Peter via QGIS-User a écrit :

Hello

I've downloaded some datasets from the LISTdata open data repository
(https://listdata.thelist.tas.gov.au/opendata/index.html) - which is a great
service.
However, I can't open any datasets that are solely an ESRI File Geodatabase,
e.g.
https://listdata.thelist.tas.gov.au/opendata/data/NCH_ES_WATER_LOGGING_HAZARD_STATEWIDE.zip
(279 MB)
Maybe these are raster gdb datasets? How can I tell if that is the case?
Is there any way the open a raster gdb other than with ESRI software?

I can open gdb files in QGIS from LISTdata that are included as part of a set with
shapefile and MapInfo tab files (e.g. in the "Tasmanian Irrigation Irrigation
District" dataset). I'm guessing that is because these are vector gdb files?

I'm running QGIS 3.10.6 ltr on Windows 10 64-bit (yes I know it's not the latest ltr, but
it was the most recent "portable" build I have been able to find and run
successfully).

Thanks
Peter
Tasmania, Australia

___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

--
http://www.spatialys.com
My software is free, but my time generally not.

___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

--
http://www.spatialys.com
My software is free, but my time generally not.

___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS (OGR) updating ESRI file geodatabases

2023-03-01 Per discussione Even Rouault via QGIS-User



Le 28/02/2023 à 18:26, chris hermansen via QGIS-User a écrit :

Thanks Randy,

If I understand Even Roualt's earlier reply it's only the FGDB version 
10 and newer that can be edited.  I was really hoping for the version 
9 ones...


If there was a strong need for version 9 editing, that would require 
further reverse engineering, which is complicated by the fact that the 
FileGDB SDK doesn't handle version 9. If my memories are correct, the 
main missing point was the .gdbindexes files that collect indices. Not 
sure about the content of the indexes themselves.


Even

--
http://www.spatialys.com
My software is free, but my time generally not.

___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS (OGR) updating ESRI file geodatabases

2023-02-28 Per discussione Even Rouault via QGIS-User




However today I stumbled upon this:

https://gdal.org/drivers/vector/openfilegdb.html

which clearly states:

The OpenFileGDB driver provides read, write and update access to
vector layers of File Geodatabases (.gdb directories) created by
ArcGIS 9 and above.

That was misleading. I've just corrected it. Write support is only for 
ArcGIS 10.x geodatabases



--
http://www.spatialys.com
My software is free, but my time generally not.
___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Two LAZ file possible issues iwth some new data we received to processs

2023-02-20 Per discussione Even Rouault via QGIS-User


To report to https://github.com/qgis/QGIS/issues

Le 20/02/2023 à 18:48, Mike Stanton via QGIS-User a écrit :


Hi Nicolas, thank you for the well researched and detailed answer.  I 
was running on the latest build and it had a the new feature for Cloud 
Optimization COPC.  That feature doubled the number of open files for 
the process and I dropped back to 3.22 to get my processing done. It 
appears that during the indexing process  of 3.28 the COPC files are 
left open even after moving on to another LAZ file. The 3.22 level did 
not create the COPC files and everything worked out for the images.


*From:* Nicolas Cadieux 
*Sent:* Friday, February 17, 2023 7:20 AM
*To:* Mike Stanton 
*Cc:* qgis-user@lists.osgeo.org
*Subject:* Re: [Qgis-user] Two LAZ file possible issues iwth some new 
data we received to processs


Hi,

I did run into this layer number limite in the past.  If you check the 
user list or if you Google “QGIS, maximum amount of layers… “ you will 
find multiple post. The issue, at the time was caused by limites in 
the OS.  (Windows 7 at the time), OSX had and easy work around.  Same 
thing for the various flavours of  Linux.


At the time, I either combined files when using QGIS or opened them in 
Saga that did not have this problem. I believe this was  because Saga, 
at least back then, loaded all the file into memory and did not keep a 
handle on the files once is was open.


This is the type of thing that got me into learning Python so that I 
could create loops and speed up operations on LiDAR data.


It’s normal (or at least expected)  for large batches of file to have 
bugs in them like missing crs in headers.  File transfers can be 
difficult to handle.  Talk with the provider to see if they have the 
same issue.  If you know a bit of Python, I could send you a script 
that reads the LiDAR header directly.


LasPy is a good little library for this 
https://laspy.readthedocs.io/en/latest/intro.html#header


You could find tools in WhiteboxTools also.

Nicolas Cadieux

https://gitlab.com/njacadieux



Le 16 févr. 2023 à 17:50, Mike Stanton via QGIS-User
 a écrit :



I have LAZ files  for a square mile of the local area where my
office is and when it is loaded under QGIS I hit two problems. It
is a bit over 2000 files and I seem to have hit a capacity limit
where it stops processing the images at about 200 files.  The
message being output for the files is “Failed to Calculate
Statistics of the Point Cloud of .LAZ”.   Breaking it
into smaller pieces gets the project around the scale limits of
the product because the doing only the rejected files gets me more
files processed, but it is not convenient considering I expect to
be processing a lot more of these files.

The other problem is that QGIS tells me I don’t have a CRS for the
images and the image has a field that says “No CRS Found” in the
log file and what I have is data in the header that says “COMPD_CS
[“Projected”, PROJCS “UTM 15N”, GEOGCS [WGS84 …. “  but the
strings have unpaired square brackets round the relevant text. Is
this what may be causing the error to occur?

___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user


___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info:https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.
___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Two LAZ file possible issues iwth some new data we received to processs

2023-02-16 Per discussione Even Rouault via QGIS-User



Le 16/02/2023 à 23:41, Mike Stanton via QGIS-User a écrit :


I have LAZ files  for a square mile of the local area where my office 
is and when it is loaded under QGIS I hit two problems. It is a bit 
over 2000 files and I seem to have hit a capacity limit where it stops 
processing the images at about 200 files.  The message being output 
for the files is “Failed to Calculate Statistics of the Point Cloud of 
.LAZ”.   Breaking it into smaller pieces gets the project 
around the scale limits of the product because the doing only the 
rejected files gets me more files processed, but it is not convenient 
considering I expect to be processing a lot more of these files.


The other problem is that QGIS tells me I don’t have a CRS for the 
images and the image has a field that says “No CRS Found” in the log 
file and what I have is data in the header that says “COMPD_CS 
[“Projected”, PROJCS “UTM 15N”, GEOGCS [WGS84 …. “  but the strings 
have unpaired square brackets round the relevant text. Is this what 
may be causing the error to occur?


Can you paste the exact CRS WKT string ? Though I doubt it is related 
with the 200 file limit issue



___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info:https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.
___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Export to ESRI File Geodatabase

2023-01-17 Per discussione Even Rouault via QGIS-User


Randal,

Support for the OpenFileGDB driver in the Save As dialog only dates back 
to a few days ago per https://github.com/qgis/QGIS/pull/51397 , so 
you'll have to wait for 3.28.3


What you must see with the driver being interminently available is that 
some of your installations might have the FileGDB driver (not 
OpenFileGDB) that relies on the ESRI closed source SDK.


Even

Le 17/01/2023 à 21:56, Randal Hale via QGIS-User a écrit :
I don't think this is just me - BUT - I wanted to check before filing 
a bug report. This is all with QGIS 3.28.2 On Windows 10 and using 
GDAL 3.6.1


I was checking on a problem and had used Network Installer to install 
QGIS - I had File Geodatabase as a choice for Export (Right Click on 
Layer -> Export As ->pick  ESRI File Geodatabase).


I installed QGIS using the Stand Alone installer (after removing QGIS 
Network Install) and didn't have ESRI File Geodatabase as a choice on 
export Features. This is confirmed by another user using the 
standalone installer.


I removed QGIS and Reinstalled using the Network Installer and didn't 
have ESRI File Geodatabase as a choice (but I did before).


So it seems like I'm missing something - I've been digging through 
options and creating new profiles but I can't make the ERI File 
Geodatabase consistently appear or disappear. Is there something else 
I should try?


Randy



--
http://www.spatialys.com
My software is free, but my time generally not.

___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] wrong interpretation of .prj-file

2023-01-02 Per discussione Even Rouault via QGIS-User


Uwe,

do you know how / by which software this .prj is generated ? The 
',AUTHORITY["ESRI",102328]' in it is quite suspicious as it is malformed 
WKT (102328 should be quoted in double quotes), and Esri software 
doesn't no generate such authority node.


Looking in the current official projected CRS list of Esri at 
https://raw.githubusercontent.com/Esri/projection-engine-db-doc/master/csv/pe_list_projcs.csv, 
which is incorporated in PROJ, ESRI:102328 is 7 Stellen and there's no 
trace of a 6 Stellen CRS definition:


102328,102328,"PE_PCS_ETRS_1989_UTM_32N_7STELLEN","ETRS_1989_UTM_Zone_32N_7stellen","PROJCS[""ETRS_1989_UTM_Zone_32N_7stellen"",GEOGCS[""GCS_ETRS_1989"",DATUM[""D_ETRS_1989"",SPHEROID[""GRS_1980"",6378137.0,298.257222101]],PRIMEM[""Greenwich"",0.0],UNIT[""Degree"",0.0174532925199433]],PROJECTION[""Transverse_Mercator""],PARAMETER[""False_Easting"",250.0],PARAMETER[""False_Northing"",0.0],PARAMETER[""Central_Meridian"",9.0],PARAMETER[""Scale_Factor"",0.9996],PARAMETER[""Latitude_Of_Origin"",0.0],UNIT[""Meter"",1.0]]","PROJCRS[""ETRS_1989_UTM_Zone_32N_7stellen"",BASEGEOGCRS[""GCS_ETRS_1989"",DATUM[""D_ETRS_1989"",ELLIPSOID[""GRS_1980"",6378137.0,298.257222101,LENGTHUNIT[""Meter"",1.0]]],PRIMEM[""Greenwich"",0.0,ANGLEUNIT[""Degree"",0.0174532925199433]],CS[ellipsoidal,2],AXIS[""Latitude 
(lat)"",north,ORDER[1]],AXIS[""Longitude 
(lon)"",east,ORDER[2]],ANGLEUNIT[""Degree"",0.0174532925199433]],CONVERSION[""Transverse_Mercator"",METHOD[""Transverse_Mercator""],PARAMETER[""False_Easting"",250.0,LENGTHUNIT[""Meter"",1.0]],PARAMETER[""False_Northing"",0.0,LENGTHUNIT[""Meter"",1.0]],PARAMETER[""Central_Meridian"",9.0,ANGLEUNIT[""Degree"",0.0174532925199433]],PARAMETER[""Scale_Factor"",0.9996,SCALEUNIT[""Unity"",1.0]],PARAMETER[""Latitude_Of_Origin"",0.0,ANGLEUNIT[""Degree"",0.0174532925199433]]],CS[Cartesian,2],AXIS[""Easting 
(X)"",east,ORDER[1]],AXIS[""Northing 
(Y)"",north,ORDER[2]],LENGTHUNIT[""Meter"",1.0]]","ETRS 1989 UTM Zone 
32N 7Stellen","Esri","10.0.1","no","Germany - west of 
12~E",47.27,55.9,5.5,12.0


However googling for ETRS_1989_UTM_Zone_32N_6Stellen, I can find a few 
references to it using the same definition as the one you point to, but 
none of them seem to reference it as ESRI:102328.


If you replace the content of the .prj file with the following one 
(which uses equivalent EPSG:25832 definition), that should work:


PROJCS["ETRS_1989_UTM_Zone_32N",GEOGCS["GCS_ETRS_1989",DATUM["D_ETRS_1989",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",50.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",9.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]

I've also submitted a fix/improvement in PROJ per 
https://github.com/OSGeo/PROJ/pull/3532, so that it trusts less 
authority codes. With that change, QGIS will use EPSG:25832 for that file.


Even

Le 02/01/2023 à 12:00, Uwe via QGIS-User a écrit :


Hello list,

first: best wishes for a peaceful 2023 to everyone!

I often have to process shapefiles that come with a prj file including 
the „Authority: ESRI 102328“ statement (complete content of that file 
see below).


This .prj is processed wrong by QGIS 3.28, because QGIS sets the false 
easting to 2.500.000 (which has to be 500.000) and it interpretes the 
projection as „7 Stellen = decimal places“ which is „6 Stellen“ in fact.


Can I change the way QGIS interpretes the .prj-File? It is annoying to 
change the projection setting for every single shapefile I load in the 
legend for each new shape delivery. In ESRI-Software, the .prj is 
being read correct.


Here is the complete .prj which is being interpreted wrong:

PROJCS["ETRS_1989_UTM_Zone_32N_6Stellen",GEOGCS["GCS_ETRS_1989",DATUM["D_ETRS_1989",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",50.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",9.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0],AUTHORITY["ESRI",102328]]

Mit freundlichen Grüßen,

Uwe Fischer

Ingenieurbüro Fischer

Esbecker Str. 8

*31036 Eime*

/T. 05182/8325/

/T. 0172/8876934/


___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info:https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.
___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Unable to open some LISTdata open data ESRI File Geodatabases

2022-12-15 Per discussione Even Rouault via QGIS-User

Hi Peter,

You could try the https://github.com/r-barnes/ArcRasterRescue utility to
convert into GeoTIFF

You could also kindly or not kindly suggest your administration to not
use proprietary formats for open data whereas thereare so many open
alternatives, like GeoTIFF...

Even

Le 15/12/2022 à 23:49, JOLLY, Peter via QGIS-User a écrit :

Hello

I'm running QGIS 3.10.6 ltr on Windows 10 64-bit (yes I know it's not the latest ltr, but
it was the most recent "portable" build I have been able to find and run
successfully).

Thanks
Peter
Tasmania, Australia

___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

--
http://www.spatialys.com
My software is free, but my time generally not.

___
QGIS-User mailing list
QGIS-User@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] [QGIS-Developer] Offering QGIS installers through Windows store

2022-11-23 Per discussione Even Rouault via Qgis-user



- There is a requirement for every library in the package to be "UWP 
compatible".
  At that time, libpq (postgresql) didn't have that, so I stopped. 
It's hard to say how many other packages will also need special 
attention. [1]


Interesting. I know that people have contributed in the past fixes to 
make PROJ *compilable* with UWP, but as we don't have a CI target for 
that, it is unknown how it behaves at runtime.


GDAL likely fails to build with UWP according to 
https://github.com/OSGeo/gdal/issues/1201


Without a CI target that can be hooked into github action (maybe that 
exists. didn't check), UWP compatibility on the long term is going to be 
hard to achieve.


--
http://www.spatialys.com
My software is free, but my time generally not.

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] multiple csv file import

2022-11-15 Per discussione Even Rouault via Qgis-user

See the X_POSSIBLE_NAMES and Y_POSSIBLE_NAMES open options of the CSV 
driver mentioned at https://gdal.org/drivers/vector/csv.html#open-options


and the last example of https://gdal.org/drivers/vector/csv.html#examples

Even

Le 15/11/2022 à 20:58, Hugh Kelley via Qgis-user a écrit :

David, this was my first thought when i saw this question as well.

however, I didn't look for very long but I haven't seen a way to tell 
ogr2ogr to read columns in a csv as the lat/lon and write those as 
points to the shapefile.  I generally write a csv to postgres as a 
non-spatial table and then process the lat lon columns with postgis.


 Are there arguments for ogr2ogr that can do this?

On Tue, Nov 15, 2022 at 12:30 PM David Strip via Qgis-user 
 wrote:


You might consider ogr2ogr as an alternative approach. You can run
this from the command line allowing you to use shell scripts to
iterate through all your .csv files. There are also python
bindings for ogr2ogr if you're more comfortable with python than
shell scripts.

On 11/15/2022 9:59 AM, Salvatore Mellino via Qgis-user wrote:

Hello,

thank you for your answer. I have many csv (about 100), so I need
an automatic procedure. Maybe a python script...

Il 15/11/2022 16:55, Nicolas Cadieux ha scritto:

Hi,

Yes, you can do that very easily using QGIS.  Layer/add
layer/add delimited text layer.  Then just export the layer in
the format of your choice. You may need to convert the
coordinates in decimal degrees (ex 75 05 30.4 ->
75.0917800 ).

You can do this in Excel using =(A6)+(B6/60)+(C6/3600)+(D6/3600)
A= Deg, B= Min, C=Sec, D= Decimal Sec.  Then export to csv.

Nicolas


On 2022-11-15 10:08 a.m., Salvatore Mellino via Qgis-user wrote:

Hello,

I would like to know if it is possible to import multiple csv
files contained in a folder and to convert them in shapefiles
(1 for each csv). All csv files are structured as "lat long
value" separate by space and without any header line.

Thank you for your help! Regards,

Salvatore




___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user



--
Hugh Kelley


___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info:https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Readonly gpkg

2022-01-31 Per discussione Even Rouault


Resolution attempt in https://github.com/qgis/QGIS/pull/47098

Even

Le 11/01/2022 à 15:10, perg...@gmail.com a écrit :


I’ve packaged a number of read-only layers (e.g. contours, 
watercourses) for my area of interest into one gpkg. However, QGIS 
writes changes to this gpkg every time it is opened, even if the 
layers themselves are never edited. I think this is related to SQLite 
wal journaling.


It’s an annoyance since it’s a large file that gets unnecessarily 
reuploaded in every incremental backup, and when I’m in the field, 
that’s sometimes over a bad 3G cell connection.


I’ve tried making the gpkg read-only, which hasn’t crashed anything in 
brief testing and generates no warnings. However, it feels risky, and 
clearly QGIS or OGR aren’t fully happy, since the temporary shm and 
wal files don’t get deleted on exit. So I’m hesitant to continue doing 
this for fear of a crash when I’m actually doing something useful.


Is there a better way to tell QGIS to not write to a gpkg in which all 
layers are only being accessed read-only?


Thanks!

---

Martin Pergler


___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info:https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Readonly gpkg

2022-01-13 Per discussione Even Rouault





I don't immediately grasp that, but it seems wrong to say that a
particular journal mode needs to be used to get the right behavior.  As
I read the sqlite3 docs, the journal mode should only affect efficiency,


We turn WAL on only because of the second bullet of 
https://www.sqlite.org/wal.html:


"WAL provides more concurrency as readers do not block writers and a 
writer does not block readers. Reading and writing can proceed 
concurrently."



--
http://www.spatialys.com
My software is free, but my time generally not.

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Readonly gpkg

2022-01-12 Per discussione Even Rouault





I think the simplest and safest change is to remove all code that
changes the journal mode as a side effect of access.  People that want
WAL can pragma it on, with DBA hat on.  In my view sqlite's WAL option
breaks the previous concept that you can have a database without
understanding databases, and as long as it's a DBA choice, that's ok.
WAL is really a programmer level concept. I'm not sure a QGIS advanced 
user/admin should be aware of that.


Separately from the "DBA chooses journal mode" view, I feel that
constantly flipping the journal mode is asking for trouble.   I suppose
one could write a test with N processes that each connect, set a random
choice of WAL or DELETE, wait a bit, do a transaction to increment a
value, perhaps repeat that, and then if they set WAL set DELETE and
exit, as a way to look for races.   Maybe sqlite3's own tests already do
that.
If you do no turn WAL on when editing, the following QGIS tests will 
fail (and some real world situations like where you have a huge layer 
being refreshed in the background while trying to edit it):


https://github.com/qgis/QGIS/commit/b6b8759efbeb833d0d3dbf6df008087701361ad3#diff-56354e2446fe2cb6d1ee92d4e984091172e964e90f3be4d3d42276e033c4986eR92

--
http://www.spatialys.com
My software is free, but my time generally not.

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Readonly gpkg

2022-01-12 Per discussione Even Rouault




So it looks like qgis is also resetting the mode to delete, or perhaps
GDAL is doing that.   I haven't found that in the sources so far.


Both attempt at doing this:

- QGIS at 
https://github.com/qgis/QGIS/blob/master/src/core/providers/ogr/qgsogrproviderutils.cpp#L1172 
to #L1230 when it closes the last connection instance of a GPKG file


- and OGR too : 
https://github.com/OSGeo/gdal/blob/master/ogr/ogrsf_frmts/sqlite/ogrsqlitedatasource.cpp#L367



Probably the most reasonable & simple change for QGIS would be turn WAL 
only when QGIS update mode is turned on. Would require testing that 
deadlocks don't occur with existing readers prior to turning WAL on (or 
accept that risk...)





___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Readonly gpkg

2022-01-11 Per discussione Even Rouault

Le 11/01/2022 à 18:58, Greg Troxel a écrit :

Jürgen E. Fischer writes:

On Tue, 11. Jan 2022 at 09:10:57 -0500, perg...@gmail.com wrote:

I've tried making the gpkg read-only, which hasn't crashed anything in brief
testing and generates no warnings. However, it feels risky, and clearly QGIS
or OGR aren't fully happy, since the temporary shm and wal files don't get
deleted on exit. So I'm hesitant to continue doing this for fear of a crash
when I'm actually doing something useful.

We do the same with world_map.gpkg - see
https://github.com/qgis/QGIS/issues/35383.

I had seen that (NetBSD 9, pkgsrc, qgis 3.16.x) and wondered if it was a
packaging error I should look into. Now I understand. I have long
found it problematic that the wal/shm files exist and persist for long
periods, and that gpkd files I don't intend to change are modified at
the OS level.

It strikes me that the basic bug needs one of two things:

1) The geopackage driver should have a read-only access mode, and qgis
should start out that way, trying to change to read/write only when
the user hits the edit button on one of the layers. And, if the file
is read only, there should be a nice error like "can't make this layer
writable" and the toggle to write should fail.

The main issue I see is that it is unknown (and non documented AFAICS in
SQLite documentation) what happens if :

- reader A has opened the file with default journaling mode (for read-only)

- reader B (other QGIS process for example, or possibly the same QGIS
process like a renderer thread) has opened the file with default
journaling mode (for read-only)

- let's say A wants to edit and we thus enable WAL to avoid
reader/writer blocks

- is B which has opened without WAL still in a safe state ? Does it
react properly to WAL having been enabled behind its back ?

The default enabling of WAL even in initial read-only access mode is to
avoid that potential issue.

There are currently way of disabling WAL (for advanced users) :

- setting the OGR_SQLITE_JOURNALenvironment variable to JOURNAL

- or setting the QGIS setting "qgis/walForSqlite" to false

See
https://github.com/qgis/QGIS/blob/master/src/core/providers/ogr/qgsogrproviderutils.cpp#L1005

With this option, there perhaps needs to be some kind of lock outside
of sqlite, which can be multireader and single writer.

2) somehow, the geopackage layer could hide all this and logically open
the layer read/write, but refrain from writing to the filesystem.
Then, when there is an actual write, the db open-for-write can
happen. After the database is written, if no layers [are
editable]or[have unwrittend changes] for some period of time (5m?) the
database should be changed back to read only. This is sort of like 1
but at the gdal layer vs the qgis layer.

Because of this I'm starting to think about moving all my data into
postgis, and not use geopackages any more.

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

--
http://www.spatialys.com
My software is free, but my time generally not.

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] How make reading from protocol to use /vsicurl_streaming/ instead of /vsicurl/?

2021-09-21 Per discussione Even Rouault


Jukka,

Somewhat counter-intuitively, you can use 
"https://api.digitransit.fi/routing/v1/routers/hsl/isochrone?fromPlace=60.168,24.938=WALK,TRANSIT=2019-11-01=08:00:00=500=1800=3600; 
as a filename when using the File source type. This is due to GeoJSON 
driver having a particular logic when passed a http[s]:// filename to 
ingest it fully in memory.


Even

Le 21/09/2021 à 12:46, Rahkonen Jukka (MML) a écrit :


Hi,

This command gives a result with ogrinfo:

ogrinfo 
/vsicurl_streaming/"https://api.digitransit.fi/routing/v1/routers/hsl/isochrone?fromPlace=60.168,24.938=WALK,TRANSIT=2019-11-01=08:00:00=500=1800=3600; 
-al


If I try to use the URI part with QGIS with Add layer – Vector – 
Protocol – HTTP/HTTPS/FTP then QGIS seems to add automatically 
/vsicurl/ into the service request. With this service /vsicurl/ does 
not work but /vsicurl_streaming/ works. It there any way to force QGIS 
to use the streaming mode?


The debugging tools do not seem to log the network traffic from the 
“read from protocol” requests. Is that a right observation?


-Jukka Rahkonen-


___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Grid shift FIle for MSI version

2021-08-31 Per discussione Even Rouault



Le 31/08/2021 à 19:56, Dario a écrit :

First of all all, thank you for your reply.

Regarding the grid issue, you know why whit the exe release I don’t 
have this kind of problem?
I'm not familiar with the Windows packaging side. Perhaps an older PROJ 
version where the existence of the grid wouldn't even be present in the 
database. But do you actually get different results when using both 
versions ?




_

Il giorno 31 ago 2021, alle ore 19:48, Even Rouault 
 ha scritto:




If I recall well, this grid is produced by Istituto Geografico 
Militare (IGM), www.igmi.org but they don't propose it for free. 
You'll have to pay to get it.


Even

Le 31/08/2021 à 19:40, Dario C a écrit :

Hi,

I'm working on a 3003 project with 7791 layers. Once I updated to 
3.16.10 msi release I received this message:
/"The preferred transfer between EPSG: 7791 - RDN2008 / UTM zone 32N 
and EPSG: 3003 - Monte Mario / Italy zone 1 is not available for use 
in the system./
/This transformation requires the lattice file 
"35160622_47161840_R40_F00.gsb", which is not available for use on 
the system.
The current transformation "Inverse of UTM zone 32N + Ballpark 
Geographic offset from RDN2008 to Monte Mario + Italy zone 1" has an 
unknown accuracy, while the current transformation "Inverse of UTM 
zone 32N + Inverse of Monte Mario to RDN2008 (5) + Italy zone 1 "has 
an accuracy of 0.1 meters."/

/
/
and actually I have translated projection between 7791 and 3003 
(around 100 meter)


I tried to install 3.16.10 exe (instead of msi) and everything works 
well.


Is it a bug or I just need a .gsb file? In this case, where can I 
find .gsb file? Thank you


cheers

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info:https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe:https://lists.osgeo.org/mailman/listinfo/qgis-user

--
http://www.spatialys.com
My software is free, but my time generally not.


--
http://www.spatialys.com
My software is free, but my time generally not.

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Grid shift FIle for MSI version

2021-08-31 Per discussione Even Rouault

If I recall well, this grid is produced by Istituto Geografico Militare 
(IGM), www.igmi.org but they don't propose it for free. You'll have to 
pay to get it.


Even

Le 31/08/2021 à 19:40, Dario C a écrit :

Hi,

I'm working on a 3003 project with 7791 layers. Once I updated to 
3.16.10 msi release I received this message:
/"The preferred transfer between EPSG: 7791 - RDN2008 / UTM zone 32N 
and EPSG: 3003 - Monte Mario / Italy zone 1 is not available for use 
in the system./
/This transformation requires the lattice file 
"35160622_47161840_R40_F00.gsb", which is not available for use on the 
system.
The current transformation "Inverse of UTM zone 32N + Ballpark 
Geographic offset from RDN2008 to Monte Mario + Italy zone 1" has an 
unknown accuracy, while the current transformation "Inverse of UTM 
zone 32N + Inverse of Monte Mario to RDN2008 (5) + Italy zone 1 "has 
an accuracy of 0.1 meters."/

/
/
and actually I have translated projection between 7791 and 3003 
(around 100 meter)


I tried to install 3.16.10 exe (instead of msi) and everything works well.

Is it a bug or I just need a .gsb file? In this case, where can I find 
.gsb file? Thank you


cheers

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] ERS Files

2021-08-15 Per discussione Even Rouault

If your ERS file is a ERMappers .ers file (see 
https://gdal.org/drivers/raster/ers.html), you should be able to open it 
with the Add Raster Layer command


Even

Le 15/08/2021 à 14:39, kelvin fox a écrit :

Hi Everyone

How do I open a ERS file containing gridded airborne magnetic data in 
QGIS 3.18.


Thank You

Kelvin

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] [QGIS-Developer] Bug with GeoPackage and Expressions?

2021-07-28 Per discussione Even Rouault

SQLite, which is the underneath database engine for GeoPackage, doesn't
support ILIKE. The SQLite LIKE operator is actually case insensitive for
ASCII characters. See paragraph 5 "The LIKE, GLOB, REGEXP, and MATCH
operators" of https://www.sqlite.org/lang_expr.html

Le 28/07/2021 à 21:43, C Hamilton a écrit :
I think I found a bug and was wondering if someone could verify this.
I am on Windows 10 and tried this with both QGIS 3.16.9 and 3.20.1.

I tried the following filter expression on the Natural Earth data set
ne_10m_populated_places
(https://www.naturalearthdata.com/http//www.naturalearthdata.com/download/10m/cultural/ne_10m_populated_places.zip
).
Once loaded into QGIS right mouse click on the layer and select
"Properties...", then click on the "Source" tab, click on "Query
Builder" and then use this query

"NAME" ILIKE '%Bom%'

This works fine. If I export ne_10m_populated_places as a GeoPackage
and try the same thing with the GeoPackage version of the layer, I get
the error:

"An error occurred when executing the query, please check the
expression syntax."

If I use LIKE instead of ILIKE then it works. It is like ILIKE was not
implemented for GeoPackage. Can anyone verify this? I think you can
use any vector layer for this test with one that is a shapefile and
the other a GeoPackage.

Thanks,

Calvin

___
QGIS-Developer mailing list
qgis-develo...@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-developer
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-developer

--
http://www.spatialys.com
My software is free, but my time generally not.

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] No transform available between XXX and EPSG:4326 when using non-Earth projections

2021-03-18 Per discussione Even Rouault



Le 18/03/2021 à 23:25, Nyall Dawson a écrit :

On Fri, 19 Mar 2021 at 07:04, Even Rouault  wrote:

Sebastian,

I can reproduce the 2 error messages you got. There's indeed something in QGIS 
that tries a terrestrial reprojection and should be fixed, but I don't think 
this is really critical apart of the annoyance of being spammed about that. 
Please file a QGIS issue to track this.

There's two parts in particular which cause this warning to be shown,
and indeed make an assumption that lat/lng === EPSG:4326
- The status bar coordinate widget
- The CRS selector "area of use" preview widget

We probably need api for testing "is non earth crs" and "given a
non-earth crs, give me a lat/long crs for the same planetary body" to
resolve these. But as Even has pointed out neither are critical
functionality, and the warning can be ignored.


At the PROJ API level, you can extract the underlying geodetic CRS 
(which in 99% of the cases will be a geographic one) of a CRS (be it a 
geographic, projected or compound one) with proj_crs_get_geodetic_crs(). 
proj_normalize_for_visualization() might be needed to get a consistent 
lon, lat axis order.


The "advantage" of this is that if you just want to see the lat, long of 
any projected CRS without requiring the datum to be WGS 84, it will work 
and without the complications of datum adjustments. And you don't need 
to distinguish Earth vs non-Earth use cases.


Pedantic note: if the CRS is a Geocentric CRS, then 
proj_crs_get_geodetic_crs() will return it, and not a geographic CRS as 
one might perhaps expect. A geodetic CRS may be a geocentric (X,Y,Z) or 
a geographic (lat,lon lon,lat) one). To get a geographic CRS from that, 
you'd need to extract its datum ( proj_crs_get_datum() and then build a 
geographic CRS with proj_create_geographic_crs_from_datum() ), but 
that's quite unlikely to be needed in QGIS use cases.




Nyall



That said, I managed to do a non-Earth reprojection by adding a custom CRS with 
the following WKT definition:

GEOGCRS["MARS Geographic",

DATUM["D_MARS",

ELLIPSOID["MARS",3396000,0,

LENGTHUNIT["metre",1,

ID["EPSG",9001,

PRIMEM["Reference meridian",0,

ANGLEUNIT["degree",0.0174532925199433,

ID["EPSG",9122]]],

CS[ellipsoidal,2],

AXIS["longitude",east,

ORDER[1],

ANGLEUNIT["degree",0.0174532925199433,

ID["EPSG",9122]]],

AXIS["latitude",north,

ORDER[2],

ANGLEUNIT["degree",0.0174532925199433,

ID["EPSG",9122

and selecting it as the project CRS. Using ESRI:104905 (GCS_Mars_2000) seems to 
also work (I haven't deeply checked but I do get a reprojected raster and 
geographic coordinates), but you'll get each time the warnings about the 2 
terrestrial coordinate operations. Hopefully in coming months, I'll be able to 
ingest the IAU catalog into the PROJ database as well.

Even

Le 18/03/2021 à 21:07, Sebastian Walter a écrit :

Hi,

Since Version 3.X I'm not able to perform transformations on planetary 
(non-Earth) data (from e.g. Mars or Mercury, example dataset: 
https://maps.planet.fu-berlin.de/mex4/h1104_.ihs.53.tif). The error message 
is:

No transform is available between SINUSOIDAL MARS and Custom CRS: 
GEOGCRS["unknown",DATUM["unknown",ELLIPSOID["WGS 8….
proj_create_operations: Source and target ellipsoid do not belong to the same 
celestial body
No transform available between SINUSOIDAL MARS and EPSG:4326

The consequence is that transformations to lat/lon are not possible any more. With the 
older QGis there was no problem like this. My suspect is the "new" Proj6 
database layout, but the above error messages also point to the fact there there is some 
EPSG:4326 projection hard-coded into Qgis.

Is anybody successfully using Qgis with non-Earth data and proper 
transformations?

Thanks for any advise...

Sebastian



___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

--
http://www.spatialys.com
My software is free, but my time generally not.

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] No transform available between XXX and EPSG:4326 when using non-Earth projections

2021-03-18 Per discussione Even Rouault


Sebastian,

I can reproduce the 2 error messages you got. There's indeed something 
in QGIS that tries a terrestrial reprojection and should be fixed, but I 
don't think this is really critical apart of the annoyance of being 
spammed about that. Please file a QGIS issue to track this. That said, I 
managed to do a non-Earth reprojection by adding a custom CRS with the 
following WKT definition:


GEOGCRS["MARS Geographic",

DATUM["D_MARS",

ELLIPSOID["MARS",3396000,0,

LENGTHUNIT["metre",1,

ID["EPSG",9001,

PRIMEM["Reference meridian",0,

ANGLEUNIT["degree",0.0174532925199433,

ID["EPSG",9122]]],

CS[ellipsoidal,2],

AXIS["longitude",east,

ORDER[1],

ANGLEUNIT["degree",0.0174532925199433,

ID["EPSG",9122]]],

AXIS["latitude",north,

ORDER[2],

ANGLEUNIT["degree",0.0174532925199433,

ID["EPSG",9122

and selecting it as the project CRS. Using ESRI:104905 (GCS_Mars_2000) 
seems to also work (I haven't deeply checked but I do get a reprojected 
raster and geographic coordinates), but you'll get each time the 
warnings about the 2 terrestrial coordinate operations. Hopefully in 
coming months, I'll be able to ingest the IAU catalog into the PROJ 
database as well.


Even

Le 18/03/2021 à 21:07, Sebastian Walter a écrit :


Hi,

Since Version 3.X I'm not able to perform transformations on planetary 
(non-Earth) data (from e.g. Mars or Mercury, example dataset: 
https://maps.planet.fu-berlin.de/mex4/h1104_.ihs.53.tif). The 
error message is:


No transform is available between SINUSOIDAL MARS and Custom CRS: 
GEOGCRS["unknown",DATUM["unknown",ELLIPSOID["WGS 8….
proj_create_operations: Source and target ellipsoid do not belong to 
the same celestial body

No transform available between SINUSOIDAL MARS and EPSG:4326

The consequence is that transformations to lat/lon are not possible 
any more. With the older QGis there was no problem like this. My 
suspect is the "new" Proj6 database layout, but the above error 
messages also point to the fact there there is some EPSG:4326 
projection hard-coded into Qgis.


Is anybody successfully using Qgis with non-Earth data and proper 
transformations?


Thanks for any advise...

Sebastian



___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user


--
http://www.spatialys.com
My software is free, but my time generally not.

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Fwd: GDAL2tiles in local crs in QGIS

2021-01-14 Per discussione Even Rouault

On jeudi 14 janvier 2021 10:07:02 CET Marco Bernasocchi wrote:
> Hi all,
> just a short/probably stupid question that I'm loosing way too much time on
> already :)
> 
> is there a way to have gdal2tiles create tiles in a local crs (for example
> EPSG:2056) and load them into QGIS?
> I'm doing
> gdal2tiles.py --profile=raster --zoom=0-5 --s_srs=EPSG:2056
> CH_over_2000_lv95.tif test_tiles
> 
> and it looks like the resulting tilemapresource.xml and the openlayers
> example are correct, but loading them in QGIS using
> 
> file:///home/marco/gis/test_tiles/{z}/{x}/{y}.png
> 
> I get a wrong placement (crs of the tiles is set to 3857, which is
> "workaroundable") and wrong tiles alignment (looks like 0,0 top left vs 0,0
> bottom left).
> 
> I know the TMS "standard" is old but it would still be nice to be able to
> quickly take a tiff, chop it up in tiles and use it in QGIS. And using
> --xyz parameter results in correct tiles but in 3857
> 
> Am I doing something stupidly wrong :) ? or is WMTS really the only way to
> have tiles in local crs?

I believe --s_srs=EPSG:2056 will be ignored in --profile=raster mode. It is 
only used for reprojecting to the SRS implied by other profiles.

In --profile=raster mode, the (0,0) tiles is the origin of your raster. Is 
that what you want ? Or do you want some tiling scheme based on EPSG:2056 
where the origin of the tiling scheme is some point generally far from your 
raster origin ? If the later, you should create a tms_X.json file and 
place it in GDAL data directory (see the existing profiles such as 
tms_MapML_APSTILE.json). Then gdal2tiles (GDAL >= 3.2) will be able to use it.

But re-reading your email, the issue might be more assumptions on QGIS side 
regarding the CRS (EPSG:3857) and the tile origin

Even


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Visualizing multi-page GeoTIFF in QGIS

2021-01-11 Per discussione Even Rouault

On lundi 11 janvier 2021 13:24:18 CET Eric Younkin - NOAA Federal wrote:
> Hello,
> 
> I have a GeoTIFF that I have created in the following manner:
> 
> gdal_translate first_grid supergrid.tif
> 
> gdal_translate second_grid supergrid.tif -co APPEND_SUBDATASET=YES
> 
> Dragging this into QGIS appears to only allow me to visualize the first
> dataset used. Is there some option to allow me to view all added pages?
> This does seem to work if I use a GeoPackage instead of a GeoTIFF.  I can
> also pull in a VRT for these TIFFs and that works.  I'm going the
> APPEND_SUBDATASET route to construct one TIFF file that will work with PROJ
> vgridshift, which appears to support multipage TIFFs.

Erik,

Looks like this is a limitation of the QGIS UI. If the GDAL driver exposes a 
valid raster when opening the file (which it does for a multi-page GeoTIFF 
file, by presenting the first page), then QGIS doesn't propose to open the 
other subdatasets. The GPKG driver is different in that when there are several 
raster tables, it doesn't expose anyone by default as a valid raster, hence 
QGIS UI proposes to select the one to open.
You may file a ticket in QGIS tracker about that, if there's not an existing 
one.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Proposal to remove DB2 provider

2021-01-10 Per discussione Even Rouault

Hi,

see https://github.com/qgis/QGIS-Enhancement-Proposals/issues/204

If you use the DB2 provider, please provider feedback in the above ticket with 
how critical this function is for you

Best regards,

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] DEM error with OpenDataNRW

2020-11-16 Per discussione Even Rouault

On lundi 16 novembre 2020 16:53:16 CET Tobias Wendorff wrote:
> Am 16.11.2020 um 16:16 schrieb Even Rouault:
>  > Please provide a link where to download this DEM1 dataset
> 
> Single tiles:
> https://www.opengeodata.nrw.de/produkte/geobasis/hm/dgm1_xyz/dgm1_xyz

Improvement to handle such datasets queued in
https://github.com/OSGeo/gdal/pull/3185

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] DEM error with OpenDataNRW

2020-11-16 Per discussione Even Rouault

Martin,

> could somebody please explain me why QGIS doesn't process the DEM1 data from
> OpenDataNRW (.xyz, 1 meter, Nordrhein-Westfalen) while the DEM25 data from
> OpenDataRLP (.xyz, 25m, Rheinland-Pfalz, but congruent structure to NRW)
> works out without any problem???

Please provide a link where to download this DEM1 dataset

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] GDA94 and WGS84 transformations

2020-10-01 Per discussione Even Rouault

Hi,
> 
> I'm having trouble understanding how QGIS handles the conversion from the
> Australian datum to WGS84 in the UTM grid.
> 
> GDA94 (EPSG:4283) /  MGA94 Zone 5X (EPSG:2835x)
> WGS84 (EPSG:4326) /  UTM South (EPSG:3265x)
> 
> As far as I know the only difference between those two is the ellipsoid
> flattening parameters:
> GDA94: 298.257222101
> WGS84: 298.257223563
> 
> When converting from WGS84 <--> GDA94 in degrees, I get exact same
> coordinates
> 
> GDA94 115.8231834665 -31.94828325747
> WGS84 115.8231834665 -31.94828325747
> 
> This is expected, since the CRSs transformation assumes the two ellipsoids
> are identical:
> 
> EPSG:1150 GDA94 to WGS 84 (1)
> Approximation at the +/- 3m level using inappropriate assumption that GDA94
> is equivalent to WGS 84. Accuracy changed from 1m to 3m due to tectonic
> plate motion over more than 15 years.
> +proj=noop

This is a much stronger assumption than the 2 ellipsoids are identical. This is 
an assumption 
that the 2 *datums* are considered identical because you don't care about 
positional 
accuracy below 3 meters.
This approximation is only used during the step where geodetic datum changes 
are used, not 
for forward or reverse map projections

> 
> However when going from degrees to UTM or MGA, I get different values:
> 
> GDA94 115.8231834665 -31.94828325747
> To
> MGA50 388777.573657*112836372* 6464692.233*647795394063*
> UTM50 388777.573657*624074258* 6464692.233
> *548169955611*
> WGS84 115.8231834665 -31.94828325747
> To
> MGA50 388777.573657112836372 6464692.233647795394063
> UTM50 388777.573657624074258 6464692.233548169955611
> 
> MGA50 388777.573657112836372 6464692.233647795394063
> To
> UTM50 388777.573657624074258 6464692.233548169955611
> 
> I don't know why UTM and MGA coordinates differ, since the transformation
> used still assumes WGS and GDA have the same ellipsoid:
> 
> INVERSE(EPSG):16150, INVERSE(EPSG):1150
> Inverse of UTM zone 50S + Inverse of GDA94 to WGS 84 (1)
> +proj=pipeline +step +inv +proj=utm +zone=50 +south +ellps=WGS84 +step
> +proj=unitconvert +xy_in=rad +xy_out=deg

If you go from GDA94 to MGA50, GRS80 will be used as the ellipsoid for map 
projection 
purposes:

$ projinfo -s GDA94 -t "GDA94 / MGA zone 50" -o PROJ
[...]

EPSG:17350, Map Grid of Australia zone 50, 0 m, Australia - 114°E to 120°E

PROJ string:
+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert 
+xy_in=deg 
+xy_out=rad +step +proj=utm +zone=50 +south +ellps=GRS80


If you go from GDA94 to UTM50S, WGS84 will be used as the ellipsoid for map 
projection 
purposes:

$ projinfo -s GDA94 -t EPSG:32750 -o PROJ --spatial-test intersects
[...]

unknown id, GDA94 to WGS 84 (1) + UTM zone 50S, 3 m, Australia - GDA

PROJ string:
+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert 
+xy_in=deg 
+xy_out=rad +step +proj=utm +zone=50 +south +ellps=WGS84


Which accounts for those submillemtric differences in projected coordinates.

But remember that if you went from GDA94 to WGS84 or a projected CRS based on 
WGS84, 
you have an inaccuracy of 1 - 2 meters... So if you care about submillemetric 
accuracy, you 
must remain in the same datum (WGS84, UTM50 / WGS84) or (GDA94, MGA50 / GDA94)

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] More user-information?

2020-06-06 Per discussione Even Rouault

> After some Q, it appeared that he had opened the shapefile 
from a ZIP...

Not addressing your general concern, but upgrading to GDAL 3.1.0 
will fix that particular point: update of .shp.zip is now supported. 
(performance will probably not be rocketting if the .shp.zip is too 
big)

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] [QGIS-Developer] (how to) register gsb transformation in proj.db?

2020-06-05 Per discussione Even Rouault

> Yes, I already thought of that, too. Can the extend be derrived 
from a
> gsb-file? Otherwise I have to estimate.

Command line with "gdalinfo Minden.gsb" and look at Corner 
Coordinates
Or just open the gsb file in QGIS as a raster and look at the 
properties.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] [QGIS-Developer] (how to) register gsb transformation in proj.db?

2020-06-05 Per discussione Even Rouault

On vendredi 5 juin 2020 11:13:09 CEST Thomas Schüttenberg wrote:
> Hi all!
> 
> Because I better should not roll out QGIS with the afore mentioned premanent
> ballpark-warning in place, I investigatetd further and stumbled across an
> "allowFallback"-setting:
> 
> By placing a default datum transformation in the settings/options.../CRS you
> get new entries in your QGIS3.ini [find it on win at
> %APPDATA%\QGIS\QGIS3\], which are the operation parameters themselve and a
> setting to allow the approximate “ballpark” transformation [...] between a
> source and destination CRS pair, in the case that the preferred coordinate
> operation fails".
> 
> [Projections]
> EPSG%3A31467\EPSG%3A25832_allowFallback=true
> EPSG%3A31467\EPSG%3A25832_coordinateOp="+proj=pipeline +step +inv
> +proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=350 +y_0=0 +ellps=bessel +step
> +proj=hgridshift +grids=Minden.gsb +step +proj=utm +zone=32 +ellps=GRS80"
> 
> This is also writen to any project file containing layers with the crs',
> covered by the default datum transformation. It appears within in the tag
>  and causes the warning
> to show up every time this project will be loaded, although the datum
> transformation is already working an no fallback seems necessary.
> 
> If the _allowFallback=-setting in QGIS3.ini is (manually) edited to 'false',
> newly saved project files also contain 
> and the warning will not be displayed again on project load.
> 
> Reference to my quite hacky findings can be found in the QGIS Python API
> Doc, search for 'allowFallback':
> https://qgis.org/pyqgis/3.10/search.html?q=allowFallback Some of the
> allowFallback stuff is labeled "New in version 3.12", however I currently
> investigate on 3.10.5-LTR. Might this indicate a backporting issue?
> (Otherwise, the implementation of my datum transformation very well might
> still be incomplete, don't know.)
> 
> Just wanted to inform you about my (current) workaround, any further hints
> are certainly appreciated. Greetings!

I'm not super familiar with the QGIS details about the selection and preference 
settings of 
transformations. I now that Nyall made a lot of efforts to make this as smooth 
as possible 
(difficult challenge to have user friendliness and accurate behaviour in that 
field!).
Looking at your record for minden.gsb, it uses a very large area of use. Once 
integrated into 
EPSG, it will probably use a much smaller one fitted to the grid extent, and 
I'm not sure how 
QGIS will react when zooming out sufficiently so that the grid doesn't cover 
the displayed 
extent.
You could simulate that by creating a new record into the area table which uses 
the real 
extent of your grid

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] (how to) register gsb transformation in proj.db?

2020-06-03 Per discussione Even Rouault

> Is it correct, that the tranformations in this dialog are read from
> \share\proj\proj.db and therefore this is the place to add my
> configuration? 

Yes

> I found that there is a record for the (not shipped)
> BWTA2017.gsb. in table grid_transformation of proj.db. So, do I just have
> to add a record to grid_transformation for my gridfile in the same way?

Yes

> What would be the correct parameters? Is there another table that has to be
> adapted to make it work?

The critical parameters to set are the EPSG code of the source and target CRS. 
You must use 
the geographic CRS ones that correspond to the projected one you mentions. So 
EPSG:4314 
(DHDN) as the source and EPSG:4258 (ETRS89) as the target. Starting from the 
record for 
BWTA2017 should be a good idea.
For the area of use, you may want to look at the "area" table for a code that 
fits the bounding 
box of the grid, or create a new one if there is no match
For the auth_name of the transformation use 'PROJ', and for the code give some 
name like 
'EPSG_4314_TO_EPSG_4258_FOR_AREA_'

Utimately, if this grid is open data matching the criteria mentioned at 
https://github.com/
OSGeo/PROJ-data, you may want to submit for inclusion here *and* also make sure 
the grid 
is registered in the EPSG database by contacting IOGP: 
http://www.epsg.org/EPSGDataset/Makechangerequest.aspx
Once donce, that would avoid you in the future to have to hack the PROJ 
database.

Note: there are a number of other grids for Germany listed at
https://github.com/OSGeo/proj-datumgrid/issues that would be candidates for 
such a 
process.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] S57 creation/export

2020-05-24 Per discussione Even Rouault

On dimanche 24 mai 2020 20:23:16 CEST Jonathan Moules wrote:
> Hi Matt,
> 
> I may be wrong, but I don't believe so. GDAL (which sits behind much of
> the format abilities of QGIS) is read-only for S-57:
> 
> https://gdal.org/drivers/vector/s57.html

Not completely true, if you go to the bottom of the doc:
https://gdal.org/drivers/vector/s57.html#s-57-export

Basically you need to provide input layers with the same structure as the ones 
of the read-
side of the driver

Some tests of S57 to S57 translation at:
https://github.com/OSGeo/gdal/blob/master/autotest/ogr/ogr_s57.py#L205

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] WFS error in 3.4.15 but not 3.4.9

2020-01-28 Per discussione Even Rouault

Lars,

Please file an issue with the URL to the server and the layer you load (check 
before if there are not similar issues opened)

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] ITRF00

2019-12-13 Per discussione Even Rouault

On vendredi 13 décembre 2019 13:01:36 CET malcolm wrote:
> The proj definition below is what I am using to convert data using PostGIS
> from NAD83 to WGS84 in a separate application. This creates data that
> matches what is transformed from ArcMap using the ITRF00 transformation
> (which is default in the latest version). The NAD83 datum in PostGIS is
> simply towgs84=0,0,0,0,0,0,0. I am using the 3735 projection for NAD83
> (Ohio SP South).
> 
> +proj=lcc +lat_1=40.03 +lat_2=38.73 +lat_0=38
> +lon_0=-82.5 +x_0=60 +y_0=0 +ellps=GRS80
> +towgs84=-0.9956,1.9013,0.5215,0.025915,0.009246,0.011599,-0.00062
> +units=us-ft +no_defs
> 

(adding back the list)

OK, so from the above parameters and digging into the EPSG database, it seems 
to be the NAD83(CORS96) to ITRF2000 transformation, with the omission of the 
time dependant terms:

$ projinfo -s "NAD83(CORS96)" -t ITRF2000 --spatial-test intersects -o PROJ
Candidate operations found: 2
-
Operation n°1:

unknown id, Conversion from NAD83(CORS96) (geog2D) to NAD83(CORS96) 
(geocentric) + Inverse of ITRF2000 to NAD83(CORS96) (1) + Conversion from 
ITRF2000 (geocentric) to ITRF2000 (geog2D), 0 m, USA - CONUS and Alaska; PRVI

PROJ string:
+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert 
+xy_in=deg +xy_out=rad +step +proj=cart +ellps=GRS80 +step +inv +proj=helmert 
+x=0.9956 +y=-1.9013 +z=-0.5215 +rx=0.025915 +ry=0.009426 +rz=0.011599 
+s=0.00062 +dx=0.0007 +dy=-0.0007 +dz=0.0005 +drx=6.7e-05 +dry=-0.000757 
+drz=-5.1e-05 +ds=-0.00018 +t_epoch=1997 +convention=coordinate_frame +step 
+inv +proj=cart +ellps=GRS80 +step +proj=unitconvert +xy_in=rad +xy_out=deg 
+step +proj=axisswap +order=2,1


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] ITRF00

2019-12-13 Per discussione Even Rouault

On vendredi 13 décembre 2019 12:00:23 CET malcolm wrote:
> Question on how QGIS chooses its transformations. I am assuming it is from
> some proj defaults?? I want to know because I do not see the current
> NAD83<>WGS84 transformation (ITRF00) which is used by HAZUS and Esri. Is
> there any way to add this upstream so it is used if not by default at least
> as an option in QGIS?

Which transformations do they use ?

I can get non-null transformations with PROJ if using NAD83(2011) as source 
CRS and ITRF2000 as target:

They go through ITRF2008 or ITRF2014 as intermediate. However as those are 
time-dependent transformations, you'd need to specify a coordinate epoch to 
get very accurate results, which is not (yet?) possible through QGIS AFAIK.


$ projinfo -s "NAD83(2011)" -t ITRF2000 --spatial-test intersects -o PROJ
Candidate operations found: 2
-
Operation n°1:

unknown id, Conversion from NAD83(2011) (geog2D) to NAD83(2011) (geocentric) + 
Inverse of ITRF2008 to NAD83(2011) (1) + Inverse of ITRF2000 to ITRF2008 (1) + 
Conversion from ITRF2000 (geocentric) to ITRF2000 (geog2D), 0.01 m, USA - 
CONUS and Alaska; PRVI

PROJ string:
+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert 
+xy_in=deg +xy_out=rad +step +proj=cart +ellps=GRS80 +step +inv +proj=helmert 
+x=0.99343 +y=-1.90331 +z=-0.52655 +rx=0.02591467 +ry=0.00942645 
+rz=0.01159935 +s=0.00171504 +dx=0.00079 +dy=-0.0006 +dz=-0.00134 
+drx=6.667e-05 +dry=-0.00075744 +drz=-5.133e-05 +ds=-0.00010201 +t_epoch=1997 
+convention=coordinate_frame +step +inv +proj=helmert +x=0.0019 +y=0.0017 
+z=0.0105 +rx=0 +ry=0 +rz=0 +s=-0.00134 +dx=-0.0001 +dy=-0.0001 +dz=0.0018 
+drx=0 +dry=0 +drz=0 +ds=-8e-05 +t_epoch=2000 +convention=position_vector 
+step +inv +proj=cart +ellps=GRS80 +step +proj=unitconvert +xy_in=rad 
+xy_out=deg +step +proj=axisswap +order=2,1

-
Operation n°2:

unknown id, Conversion from NAD83(2011) (geog2D) to NAD83(2011) (geocentric) + 
Inverse of ITRF2014 to NAD83(2011) (1) + Inverse of ITRF2000 to ITRF2014 (1) + 
Conversion from ITRF2000 (geocentric) to ITRF2000 (geog2D), 0.01 m, USA - 
CONUS and Alaska; PRVI

PROJ string:
+proj=pipeline +step +proj=axisswap +order=2,1 +step +proj=unitconvert 
+xy_in=deg +xy_out=rad +step +proj=cart +ellps=GRS80 +step +inv +proj=helmert 
+x=1.0053 +y=-1.9092 +z=-0.5416 +rx=0.0267814 +ry=-0.0004203 +rz=0.0109321 
+s=0.00037 +dx=0.0008 +dy=-0.0006 +dz=-0.0014 +drx=6.67e-05 +dry=-0.0007574 
+drz=-5.13e-05 +ds=-7.01e-05 +t_epoch=2010 
+convention=coordinate_frame +step +inv +proj=helmert +x=-0.0007 +y=-0.0012 
+z=0.0261 +rx=0 +ry=0 +rz=0 +s=-0.00212 +dx=-0.0001 +dy=-0.0001 +dz=0.0019 
+drx=0 +dry=0 +drz=0 +ds=-0.00011 +t_epoch=2010 +convention=position_vector 
+step +inv +proj=cart +ellps=GRS80 +step +proj=unitconvert +xy_in=rad 
+xy_out=deg +step +proj=axisswap +order=2,1


Even


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS 3.10 - wgs84 (epsg4326) file disappears in pseudo-mercator-project (epsg:3857)

2019-11-24 Per discussione Even Rouault

On dimanche 24 novembre 2019 16:40:02 CET Grün, Michael wrote:
> Hi all,
> since QGIS 3.10 I have a problem with in wgs84 projected files. When I
> change the whole qgis-project to pseudo-mercator (epsg:3857) it was up to
> QGIS 3.8. no problem to draw the wgs84-projected files correctly within
> epsg:3857-projects. Now in QGIS3.10. the wgs84-projected files
> disappear...the layer is turned on, but the data are not visiblethey
> only get drawn when I change the whole project to wgs84 (epsg:4326), which
> is not comfortable for cartographic reasons... The settings options for crs
> in 3.10 are the same as in 3.8...what changed there? Thanx for any
> hintsbest regards, Michael

Michael,

It sounds pretty much like https://github.com/qgis/QGIS/issues/32973 that has 
been fixed 2 days ago in master, and has been backported to the 3.10 branch as 
well for the next point release.
You might try the latest nightly build of qgis-dev of OSGeo4W to see if that 
fixes things.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Crash when opening project with QGIS 3.10

2019-10-30 Per discussione Even Rouault

> What's this about? A social experiment?

Yes, I'm "Ben Hutcher". Too bad you destroyed it so quickly :-)

I just discovered this functionnality of QGIS and this scared me. I believe I 
would have accepted the warning because I ignored that QGIS projects could 
contain Python code, and I presume a lot of users would. IMHO clicking on the 
Enable macros link should raise a dialog box with a more explicit message 
about the potential risks to double confirm.

"Python macros cannot currently be run." sounds to me as "you should enable 
macros so things work as expected". It is not even clear that the macros come 
from the project itself. They could be some trusted code in QGIS itself. Once 
you know that functionnality is, then yes the current behaviour is probably 
fine. But if you don't know it, there's a high chance you could run untrusted 
code without realizing it.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS and WFS limitations

2019-10-29 Per discussione Even Rouault

Lars,

https://issues.qgis.org/issues/6535 is about QGIS server side, and here if I 
understand well, you are considering QGIS client side, so completely different 
parts of the code

> Did the developers forget to implement WFS POST in 3.x ?

It seems https://github.com/qgis/QGIS/issues/32348 was created recently about 
that.
If you can add a way to reproduce your issue in that ticket (server + QGIS 
filter), that could be useful

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] GeoPackage deadlocks

2019-10-25 Per discussione Even Rouault

On vendredi 25 octobre 2019 22:51:39 CEST Jésahel Benoist wrote:
> I could provide some background information as we're working for about a
> year in a multi-users (5 users) configuration mainly based on network
> shared geopackages in a Windows environment. I've tried and tested a lot of
> configurations and still looking for a near PostGIS solution only in terms
> of easiness and security.
> 
> At first, we were using SQLite and I have to say that it was good,
> especially in terms of speed. But it lacked the possibility to alter
> database schema, which was a pain.
> So I switched to Geopackage and very bad problems started.
> First as you all know was the speed problem. Sometimes, it was just
> unusable with only two users, and sometimes even with only one user,
> especially with big database, with some display locks. I discovered that
> the drawing by itself was like a concurrent user... Disabling the drawing
> restored usage, even if sometimes it needed to restart qgis.
> Then I discovered the wal problem and switched to delete mode. Speed
> problem was gone but... I got a file lock. One of my db file was then
> unusable. I had to restore it. Happily it was with the most experienced
> user and at the end of the day, and we had only one session data lost.
> So I decided to divide the biggest db to several little files and to make
> apply stricter usage rules to users to be sure that there isn't much
> concurrent access.
> For now it's usable but we still have to restart qgis sometimes. Sometimes
> the display disappears and reappears with a delay or by creating an entity.
> Note that it's the same with SQLite.
> 
> As a coder, I wonder why there isn't always a file locking mechanism with a
> threaded waiting loop, even for read access and for all file formats,
> because if you change the file content, nothing can say a read wasn't done
> at the same time, and so the read could obtain bad values. Even multiple
> threads on a single machine could access the same file anytime, so it would
> be safer to avoid this kind of problem.

To disable all locking, you can try setting the SQLITE_USE_OGR_VFS environment 
variable to YES that will switch to a custom SQLite virtual file system 
implementation (designed for GDAL /vsicurl/ mostly, but that can work with 
regular files as well), that doesn't take any lock at all. But beware: it will 
eat your children if you try to do concurrent edits, and you'll get reading 
errors if reading while editing.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] GeoPackage deadlocks

2019-10-25 Per discussione Even Rouault

On vendredi 25 octobre 2019 17:10:40 CEST Matthias Kuhn wrote:
> Hi all,
> 
> We had the same experience in a project of ours. Multiple users with
> read-only access to a gpkg on a (Windows) network share. The results
> were freezes on unloading the layer (e.g. closing QGIS or opening
> another project) on system A until system B would have closed the QGIS
> layer there. The behavior seemed surprising to me and on first sight I'd
> have put it more into the category "bug" rather than "feature".
> 
> Good news, this behavior could be changed by switching the sqlite
> journal mode to delete (environment variable OGR_SQLITE_JOURNAL=DELETE).

Caution: WAL was enabled because DELETE (which is the default journal mode of 
SQLite) was causing locking issues, even in local file, single user, single 
QGIS application context. There might have been later changes in the OGR 
provider that makes it no longer needed, but this should be carefully checked. 
In particular if the test added in https://github.com/qgis/QGIS/commit/
b6b8759efbeb833d0d3dbf6df008087701361ad3#diff-17fd782357fcfed115b06d24b4416bbdR92

still pass in DELETE

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Version 3.4.7 (LTR) is locking MapInfo tab files

2019-06-25 Per discussione Even Rouault

On mardi 25 juin 2019 13:34:47 CEST Paul Wittle wrote:
> Hi,
> 
> I believe there may already be a bug ticket for this but I've found my
> deployment of QGIS is causing MapInfo errors when opening MapInfo tab
> files. It happens when you load the TAB file and also when you open the
> project and it even happens when a TAB file is open but not rendering.

Was reported at https://github.com/qgis/QGIS/issues/29885

What is annoying is that I fixed that issue 2 years ago per
https://github.com/qgis/QGIS/pull/3034.
I can't see what has changed on GDAL or QGIS side regarding how and when files 
are opened in read-only vs read-write mode since then. So perhaps this would 
be related to a change on MapInfo side how it opens files in newer versions.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Geolocating .flt files

2019-05-10 Per discussione Even Rouault

On jeudi 9 mai 2019 19:31:55 CEST Quinn, Gary wrote:
> I am new to Qgis.
> 
> I have an .flt (floating point raster) file and an associated .hdr file.
> I can import the .flt  image file into Qgis (using Add raster layer), and
> the image displays correctly.
> But how do I geolocate this image using the associated .hdr file?
> In other words, how do I make Qgis use the .hdr file to let me know the
> lat/lon of each pixel on the display?
> 
> Note: The .hdr file is like this:
> geo points  = { 1, 1,  -40.212464,  146.355212, 1,  2590,
> -44.006326,  145.100364,  1024,  2590,  -44.252673,  146.901021,  1024,
>  1,  -40.442148,  148.049722}
> 
> So I want Qgis to use this .hdr information to geolocate the image.
> I have searched various info about using .flt images and geolocation but
> can't find anything.

Gary,

Those are ground control points in the ENVI format:
https://www.harrisgeospatial.com/docs/ENVIHeaderFiles.html

GDAL can write those, but not read them, hence the georeferencing is ignored. 
I've filed https://github.com/OSGeo/gdal/issues/1528 about that

In the meantime you could use the georeferencer to add manually the ground 
control points.

As documented, the values in geo points are a succession of
pixel_x_coordinate, pixel_y_coordinate, latitude, longitude
Note: you must substract 1 to pixel_x_coordinate and pixel_y_coordinate since 
the GDAL/QGIS convention for origin of a raster is (0,0) whereas ENVI uses 
(1,)

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Import of Large CSV File

2019-04-16 Per discussione Even Rouault

On mardi 16 avril 2019 16:06:36 CEST Grant Boxer wrote:
> I have a 2.5 GB csv file which contains 9 fields and 17,167,595 records.
> QGIS does not seem to be able to import this using the "Delimited Text"
> option. Any suggestions?

Try opening it with the "Add vector layer" functionality. This will go through 
the OGR CSV driver.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] [gdal-dev] Convert KML to Shapefile PyQGIS

2019-03-19 Per discussione Even Rouault

Julierme,

the issue here is that your KML file is seen as a multi-layer dataset, with
one feature per layer, so you can't use the logic you put below. You need to
merge all the layers in a single one.

You can do that with the ogrmerge.py script that comes with GDAL:
https://www.gdal.org/ogrmerge.html

$ ogrmerge.py -overwrite_ds -single -src_layer_field_name src_layer -o out.shp 
A_Cycle06.kml

With released versions of the script, this will be very slow. I've just pushed 
a fix to speed that up, so grab:
https://raw.githubusercontent.com/OSGeo/gdal/master/gdal/swig/python/scripts/ogrmerge.py

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Supported versions of Geopackage

2018-10-21 Per discussione Even Rouault

On dimanche 21 octobre 2018 18:44:36 CEST Stéphane Henriod wrote:
> Dear all
> 
> I have a quick question regarding the support for Geopackage.
> 
> I want to prepare a few geopackages in Qgis and then serve their data as
> WMS from GeoServer. But apparently GeoServer only supports Geopackage
> version 1.2, while Geopackages created with Qgis are an older version.

Doh. What are they so strict on this ? GeoPackage 1.0, 1.1 or 1.2 are pretty 
similar, just bringing up more capabilities. So a v1.2 reader should be able 
to handle previous versions without effort.

> Consequence: my layers are not recognized by GeoServer.
> 
> Is this a known problem or is there a trick that I am missing?

If you upgrade to latest OSGeo4W with GDAL 2.3.2, it will create Geopackage 
v1.2 by default
Otherwise you can use ogr2ogr to force the GPKG version

ogr2ogr out_v1.2.gpkg in.gpkg -f GPKG -dsco VERSION=1.2

Even


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Potential Bug in writing gpkg from database manager?

2018-10-06 Per discussione Even Rouault

On vendredi 5 octobre 2018 17:00:12 CEST Michael Treglia wrote:
> I'll file a bug when I can deal with an issue logging into the issue
> tracker system, but if somebody else can do it quickly and beats it to me,
> just report back. Thanks all!

Done as https://issues.qgis.org/issues/20040
I'm working on a fix

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] WFS failure

2018-06-15 Per discussione Even Rouault

Jeroen,

It looks like this server doesn't properly implement WFS 2.0, and wrongly 
expects a TYPENAME parameter to be provided for the GetFeature request, 
whereas QGIS provides the TYPENAMES one as per the standard.

The solution is in the window defining the connexion parameters to force 
version 1.0 or 1.1, instead of letting the automatic version selection

Even


> Hi all,
> 
> Using Qgis 2.18.20 (64bit, windows) I can't download/open a specific wfs
> layer. Using 2.14.18 (64bit, windows) the same WFS link works without a
> problem.
> 
> What is going wrong in 2.18.20?
> The error is the following (in Dutch)
> Laag se:OGC_Gasnetbeheerdervlak: Downloaden van objecten voor laag
> se:OGC_Gasnetbeheerdervlak mislukt of gedeeltelijk mislukt: Server
> genereerde een uitzondering in het antwoord GetFeature: Internal error. U
> kunt met F5 proberen om de laag opnieuw te laden
> 
> The wfs link is the following (all public data)
> https://ogc.polygonentool.nl/wfs?
> Username: MijnAansluitingOGC
> Password: MApoly2017
> 
> Thanks in advance.
> 
> Jeroen
> 
> ___
> Qgis-user mailing list
> Qgis-user@lists.osgeo.org
> List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
> Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] How to Prevent GeoPackage file from updating modified date when viewed in a Qgis project when it is not being edited

2018-04-10 Per discussione Even Rouault

On mardi 10 avril 2018 08:11:30 CEST Richard Duivenvoorde wrote:
> On 09-04-18 22:22, Even Rouault wrote:
> > Or perhaps have some option in a dialog box somewhere! "I promise that I
> > will not edit this file in any concurrent session, or in this session, so
> > you can safely open it in DELETE mode" ? Or spend more time digging again
> > into the details of sqlite locking to see if we could make a clever use
> > of it, and/or possibly improving it in sqlite itself (rather involved!)
> 
> I would be ok with an option in the dialog, though it is hard to promise
> that somebody *else* will not open it too in a non readonly mode, in a
> networked work environment ...

You'll have issues with sqlite in a networking environment, as WAL cannot work 
in that context (QGIS does some effort to detect if the file location is a 
network driver or not, before deciding to turn WAL on). So you'll get the 
default journaling behaviour, which is lock prone.

> 
> Is it a possibility (at least on linux) to look at the file permissions?
> That would open the possibility to put a read-only (permissioned) file,
> and in that way you would promise to not edit/delete?

Just tested: QGIS can open file a read-only .gpkg on Linux (it will fallback 
to the default journaling mode), and presumably on Windows as well.

Or were you talking about some side car files that would be used to 
application level locking ?

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] How to Prevent GeoPackage file from updating modified date when viewed in a Qgis project when it is not being edited

2018-04-09 Per discussione Even Rouault

On lundi 9 avril 2018 15:24:41 CEST rjwill...@gmail.com wrote:
> I would really like to use and promote GeoPackage as the go to data format
> in my organization but I am having the following problem that I haven't been
> able to solve.
> 
> 
> 
> When ever someone opens a Q template project that has a layer with a
> GeoPackage data source, the GeoPackage file (not its contents) will be
> altered and it will show that it has been updated in the Windows file
> manager. I have been told that this happens because the write ahead logging
> has been triggered (creation of *.wal and *.shm) files by simply viewing the
> contents of the GeoPackage in QGIS and/or in Intramaps Roam (works on top
> of QGIS for tablets).
> 
> 
> 
> This is problematic because most users in our organization will then look at
> the windows modified date time stamp and conclude that the file's contents
> have been modified. Additionally, all of the GIS layers are stored in
> Dropbox and subsequently Dropbox concludes that the file has been modified
> and it re-uploads a new copy of the GeoPackage. When it is a large
> GeoPackage > 1 GB in size this is quite problematic and causes a cascade of
> uploads and downloads.
> 
> 
> 
> Over the last few months I have tried to suppress this modification for the
> large GeoPackages in our system in the following ways: (1) make the layer
> read only in the QGIS project (no effect); (2) make the Dropbox folder read
> only for all users except the Dropbox administrator (this is effective for
> the large base layers); (3) change the properties of the file in windows to
> read only (effective for the Dropbox administrator); and (4) add a user
> variable to the windows system environment variables: OGR_SQLITE_JOURNAL =
> DELETE (effective for QGIS but not for Intramaps Roam)
> 
> 
> 
> Combining 2, 3, and 4 allows me to use GeoPackage for layers I can ensure
> are placed in the read-only Dropbox folders.
> 
> 
> 
> It is likely that my organization will never need the simultaneous editing
> of datasets that is facilitated by Write-Ahead Logging, but I would really
> like to be able to replace shapefiles and File geodatabases with GeoPackages
> in all of our projects going forward.
> 
> 
> 
> Therefore, is there a simpler way to suppress the write ahead logging
> behaviour that causes the files to believe they have been updated for all
> GeoPackages created in QGIS or used in a QGIS project?

My memory is already fading about the details, but the reason for the current 
behaviour is that if you open a GeoPackage database in read-only non-WAL mode 
and someone else try to edit it concurrently, unpleasant locking will occur. 
Actually due to how QGIS works, locking could occur even in the same session 
(since background rendering threads can open read-only connections while the 
main connection is writing data). See https://issues.qgis.org/issues/15351

So the solution was:

commit b6b8759efbeb833d0d3dbf6df008087701361ad3
Author: Even Rouault <even.roua...@spatialys.com>
Date:   Thu Sep 22 23:23:00 2016 +0200

Fix database locking when editing GeoPackage

Concurrent read and write can lock a GeoPackage database given
the default journaling mode of SQLite (delete). Use WAL when
possible to avoid that.

Fixes #15351

I understand your use case as well, but I don't have a solution in mind that 
would fit all. 
I don't remember if I tried scenarios where one connection would open in 
DELETE mode and another one would try to turn to WAL on. 
A quick try shows that the a connection that opens in DELETE mode doesn't see 
the change to WAL done by another one. And WAL documentation doesn't describe 
this scenario, so it is quite unknown what would happen. So the current 
behaviour is to be on the safe side and always turn on WAL on opening.

Or perhaps have some option in a dialog box somewhere! "I promise that I will 
not edit this file in any concurrent session, or in this session, so you can 
safely open it in DELETE mode" ? Or spend more time digging again into the 
details of sqlite locking to see if we could make a clever use of it, and/or 
possibly improving it in sqlite itself (rather involved!)

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] GeoJSON Dropping Nulls

2017-07-07 Per discussione Even Rouault

On vendredi 7 juillet 2017 16:16:18 CEST Michele M Tobias wrote:
> I'm seeing some strange behavior working with geojson files in newer
> versions of QGIS (2.18.9 and 2.18.10, but not 2.18.6) with GDAL 2.2.0 (but
> not 2.1.3). When I use "save as" saving to geojson or edit a geojson file
> and save the edits, any null values disappear.  This is particularly a
> problem for a project I'm working on where I need to keep a standard
> attribute table and fill in nulls for any missing values to be updated in
> later.  I have more details about the problem here, including example
> files: https://github.com/UCDavisLibrary/ava/issues/270 It's not GDAL alone
> that's causing the problem - I tried saving a file into geojson in the
> osgeo4w shell with ogr2ogr and it didn't cause the same problem.
> 

We discussed about that with Alex Mandel on QGIS IRC a couple hours ago.
He filed https://issues.qgis.org/issues/16812 about that

> Thanks for your help!
> 
> Michele
> 
> Michele Tobias, PhD
> GIS Data Curator
> Data Management Program
> UC Davis Library
> 
> 161 Shields Library
> (530)752-7532
> mmtob...@ucdavis.edu
> https://www.library.ucdavis.edu/service/data-management/


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] [QGIS-Developer] Geodesic Measure Tool With Bearing Information

2017-06-14 Per discussione Even Rouault

On jeudi 15 juin 2017 07:41:06 CEST Nyall Dawson wrote:
> On 13 June 2017 at 00:20, C Hamilton  wrote:
> >> > Right now it could use the Vincinty algorithms in QGIS to calculate the
> >> > distances, but eventually I would like to see the geographiclib
> >> > algorithms
> >> > used.
> >> 
> >> I'd start with the bearing changes (since they are self contained, and
> >> the methods for calculating bearing are already present in the QGIS
> >> api). Then using geographiclib could be discussed as a separate set of
> >> changes, to be decided on its own merit. If you mix these two changes
> >> it's likely to get bogged down as one or other part of the changes are
> >> debated.
> > 
> > Ok.
> 
> This looks relevant:
> 
> https://github.com/sourcepole/kadas-albireo/commit/df64318e5198f40f0e96d27fc
> f9e179fbbf682fd
> 
> There's a lot of unrelated changes in there, but it also shows the
> basics for using geographiclib within the QGIS c++ code.
> 

Note that at least since proj 4.9.0, Charles Karney has ported its geodesic 
code from 
geographiclib in C (since proj.4 is C only) and the header is normally 
installed with proj-dev 
package (I've just checked on my Ubuntu 16.04 with proj 4.9.2. I see it is also 
available in 
OSGeo4W):

https://github.com/OSGeo/proj.4/blob/4.9.0/src/geodesic.h

Even
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Confirming Problem Importing Geopackage Data?

2017-06-07 Per discussione Even Rouault

On mercredi 7 juin 2017 13:09:53 CEST Tobias Wendorff wrote:
> Am Di, 6.06.2017, 22:11 schrieb Even Rouault:
> > Renaming the layers in ArcMap probably fixed the case inconsistency.
> > 
> > Fixed per https://trac.osgeo.org/gdal/ticket/6916
> 
> I don't know, if .gpkg allows UTF8 characters in column name.
> Doesn't lower() introduce problems with UTF8? Some characters might
> turn into a different notation? Don't know exactly, just want to
> give a hint ;)

Indee from https://sqlite.org/lang_corefunc.html#lower , 
"The lower(X) function returns a copy of string X with all ASCII characters 
converted to lower 
case. The default built-in lower() function works for ASCII characters only."
it is not very clear what happens when non ASCII characters are encountered.

Little experiment :

$ sqlite3
sqlite> select lower('EéE');
eée

So apparently it lets non-ASCII UTF-8 characters unchanged, so I think using 
lower() should be 
safe

Even


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Confirming Problem Importing Geopackage Data?

2017-06-06 Per discussione Even Rouault

On mardi 6 juin 2017 15:13:20 CEST Michael Treglia wrote:
> Thanks Even.
> 
> I renamed layers within Arc, just for anonymity of the data, and then tried
> to open the data again in QGIS 2.18.9, and it worked, so not sure what was
> going on.
> 
> If I can get another demonstrable and replicable example, I'll file a
> ticket.

I've fixed the issue with the input provided by Michael. The issue was that 
there was an 
inconsistency in the case of the table name between the table_name columns of 
gpkg system 
tables, and the table name in sqlite_master system table.
I'm not completely sure if such .gpkg are strictly valid, but as SQLite3 
doesn't accept both foo 
and Foo as different tables, comparisons can be made case insensitive to avoid 
such issues.

Renaming the layers in ArcMap probably fixed the case inconsistency.

Fixed per https://trac.osgeo.org/gdal/ticket/6916

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Confirming Problem Importing Geopackage Data?

2017-06-06 Per discussione Even Rouault

On mardi 6 juin 2017 14:41:19 CEST Michael Treglia wrote:
> Hi All,
> 
> I have a gpkg file into which I exported a bunch of layers from ArcGIS.
> When drag the gpkg file into the QGIS desktop (2.18.8 & 2.18.9), most
> layers show up as not having a geometry type, though the geometries appears
> fine in both Arc (10.3.1) and QGIS 2.14.3, on a different PC.
> 
> Can anybody confirm that there's an issue?

Please file a ticket about that with a dataset demonstrating the issue attached 
(or a link to it)

Even

> 
> I'll note my 2.18.9 was installed via OSGeo4w, compiled against GDAL 2.2.0;
> I'm not sure how my 2.14.3 was installed, as it's on a separately-managed
> server, but it's compiled against GDAL 2.0.2. (SpatialLite version 4.3 in
> both cases).
> 
> I'm not seeing this specific issue reported in the bug tracker, but might
> be missing it, and wanted to make sure I'm not missing something else
> obvious.
> 
> Thanks!
> Mike


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Unable to open XTF (INTERLIS 2) file in QGIS

2017-03-17 Per discussione Even Rouault

On vendredi 17 mars 2017 07:48:05 CET Germán Carrillo wrote:
> Hi All,
> 
> I'm attempting to visualize the data contained in a XTF file in QGIS (OGR
> provider). I'm using these data [1], specifically the Lake polygon layer,
> which has 3 features.
> 
> QGIS displays the geometries but cannot distinguish among features, i.e.,
> if I select one, all 3 lakes are selected. Moreover, all 3 records in the
> attribute table have the same attributes. Strangely, printing each
> feature's id using PyQGIS gives me -1.
> 
> One could think it's an OGR issue; however, the following OGR command is
> able to distinguish the 3 features properly, each one with its own id and
> attributes:
> 
> $ ogrinfo -ro data/CHBaseEx_WaterNet_V1.xtf CHBaseEx_WaterNet_V1.WaterNet_
> WithModificationObjects.Lake
> 
> 
> I've confirmed this issue in QGIS v2.14.x, 2.18.x, and 2.99 (built on
> Monday). In the latter, the attribute table shows only one record, although
> QGIS displays 3 polygons.
> 
> 
> Is this a known problem? Should I create an issue in QGIS Issue Tracking?

This is an issue in the OGR Interlis 2 driver that didn't set a feature id to 
the features, so it 
was left to its -1 default value (you could see it on the ogrinfo report)

Fixed per https://trac.osgeo.org/gdal/ticket/6839

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Limit to Geodatabase File Size?

2017-03-08 Per discussione Even Rouault

> I've never used GPKG but I do use spatialite all the time.  I always create
> without the index then create the index after I'm done.  It makes sense
> because if you start with the spatial index then it has to reshuffle the
> RTree with each new item.
> 
> In c++ code I set SPATIAL_INDEX=DEFFERED so that the index is built at the
> end.  Even would this work from the command line too?

The GPKG driver implements this deffered behaviour when you set SPATIAL_INDEX 
to YES (or 
doesn't specifiy it since it is the default behaviour). 

Basically it is 2 statements :

CREATE VIRTUAL TABLE rtree_tablename_geomcol
USING rtree(id, minx, maxx, miny, maxy)

--> instantaneous


Populate the spatial index from the bounding box of geometries:

INSERT OR REPLACE INTO rtree_tablename_geomcol
 SELECT fid, st_minx(geom), st_maxx(geom), st_miny(geom), st_maxy(geom)
FROM tablename

--> takes hours. Had to stop it since it made my computer slow.

I'm not sure of the reason of this very slow performance. Perhaps the I/O 
performance of my 
disk sucks, or there's something intrinsically slow to build the spatial index 
given the dataset.

> Also maybe this
> isn't available anymore I'm still using 1.11.x.

Yes spatial index support in GPKG requires at least GDAL 2.0

Even

> 
> 
> > Jonathon
> > 
> > On Tue, Mar 7, 2017 at 9:38 PM, Even Rouault <even.roua...@spatialys.com>
> > 
> > wrote:
> > > On mardi 7 mars 2017 14:30:13 CET Parker, Jonathon wrote:
> > > > New user to QGIS here. Went through some tutorials, had ESRI
> > > > experience
> > > 
> > > in
> > > 
> > > > distance past. Starting a demonstration project using data from the US
> > > > 
> > > > Census Bureau https://www.census.gov/geo/maps-data/data/tiger.html
> > > > 
> > > > 
> > > > 
> > > > I successfully added shapefiles and small geodatabases
> > > > 
> > > > ftp://ftp2.census.gov/geo/tiger/TGRGDB16// as layers. When I tried the
> > > > 
> > > > same process with a large geodatabase (tlgdb_2016_a_us_edges.gdb), it
> > > > did
> > > > 
> > > > not render. The feature count was reported as 0. I contacted the US
> > > 
> > > Census
> > > 
> > > > bureau and they verified it was not corrupted and worked on ESRI
> > > 
> > > products.
> > > 
> > > > Is there a size limit that prevents QGIS from working with large
> > > > files?
> > > > 
> > > > The table in the geodatabase is 13GB.
> > > 
> > > I've downloaded the file and reproduced the issue. The OGR OpenFileGDB
> > > driver had a bug in the particular situation where the attribute
> > > description section of the .gdbtable is beyond the 4 GB offset. I've
> > > just
> > > fixed it per
> > > 
> > > https://trac.osgeo.org/gdal/ticket/6830
> > > 
> > > 
> > > 
> > > The feature count now reported is ~68 million and it displays fine in
> > > QGIS
> > > (against the fixed GDAL version)
> > > 
> > > 
> > > 
> > > Best regards,
> > > 
> > > 
> > > 
> > > Even
> > > 
> > > 
> > > 
> > > --
> > > 
> > > Spatialys - Geospatial professional services
> > > 
> > > http://www.spatialys.com
> 
>  
>  


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Limit to Geodatabase File Size?

2017-03-08 Per discussione Even Rouault

On mercredi 8 mars 2017 12:46:54 CET Parker, Jonathon wrote:
> Even,
> 
> Thanks for the quick turn around.  How do I get the bug fix?  Do I download
> and install 2.18.4 or is there an easier way to get the bug fix?

You'll need to wait for the fix to land in an official GDAL release, so not 
immediately.

Or you can use the today's gdal-dev daily build from OSGeo4W that must have the 
fix, but you'll
need to convert the file to something else before with ogr2ogr
For example:
ogr2ogr -f GPKG tlgdb_2016_a_us_edges.gpkg tlgdb_2016_a_us_edges.gdb -lco 
SPATIAL_INDEX=NO

The SPATIAL_INDEX=NO option is due to my experimentation. There's a strong 
performance issue
when creating a SQLite spatial index on such a big table. Which is a bit sad 
given that it would
be the most useful in that situation...

> 
> Jonathon
> 
> On Tue, Mar 7, 2017 at 9:38 PM, Even Rouault <even.roua...@spatialys.com>
> 
> wrote:
> > On mardi 7 mars 2017 14:30:13 CET Parker, Jonathon wrote:
> > > New user to QGIS here. Went through some tutorials, had ESRI experience
> > 
> > in
> > 
> > > distance past. Starting a demonstration project using data from the US
> > > 
> > > Census Bureau https://www.census.gov/geo/maps-data/data/tiger.html
> > > 
> > > 
> > > 
> > > I successfully added shapefiles and small geodatabases
> > > 
> > > ftp://ftp2.census.gov/geo/tiger/TGRGDB16// as layers. When I tried the
> > > 
> > > same process with a large geodatabase (tlgdb_2016_a_us_edges.gdb), it
> > > did
> > > 
> > > not render. The feature count was reported as 0. I contacted the US
> > 
> > Census
> > 
> > > bureau and they verified it was not corrupted and worked on ESRI
> > 
> > products.
> > 
> > > Is there a size limit that prevents QGIS from working with large files?
> > > 
> > > The table in the geodatabase is 13GB.
> > 
> > I've downloaded the file and reproduced the issue. The OGR OpenFileGDB
> > driver had a bug in the particular situation where the attribute
> > description section of the .gdbtable is beyond the 4 GB offset. I've just
> > fixed it per
> > 
> > https://trac.osgeo.org/gdal/ticket/6830
> > 
> > 
> > 
> > The feature count now reported is ~68 million and it displays fine in QGIS
> > (against the fixed GDAL version)
> > 
> > 
> > 
> > Best regards,
> > 
> > 
> > 
> > Even
> > 
> > 
> > 
> > --
> > 
> > Spatialys - Geospatial professional services
> > 
> > http://www.spatialys.com


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Limit to Geodatabase File Size?

2017-03-07 Per discussione Even Rouault

On mardi 7 mars 2017 14:30:13 CET Parker, Jonathon wrote:
> New user to QGIS here.  Went through some tutorials, had ESRI experience in
> distance past.  Starting a demonstration project using data from the US
> Census Bureau https://www.census.gov/geo/maps-data/data/tiger.html
> 
> I successfully added shapefiles and small geodatabases
> ftp://ftp2.census.gov/geo/tiger/TGRGDB16// as layers.  When I tried the
> same process with a large geodatabase (tlgdb_2016_a_us_edges.gdb), it did
> not render. The feature count was reported as 0.  I contacted the US Census
> bureau and they verified it was not corrupted and worked on ESRI products.
> 
> Is there a size limit that prevents QGIS from working with large files?
> The table in the geodatabase is 13GB.

I've downloaded the file and reproduced the issue. The OGR OpenFileGDB driver 
had a bug in 
the particular situation where the attribute description section of the 
.gdbtable is beyond 
the 4 GB offset. I've just fixed it per
https://trac.osgeo.org/gdal/ticket/6830

The feature count now reported is ~68 million and it displays fine in QGIS 
(against the fixed 
GDAL version)

Best regards,

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] [gdal-dev] converting ECW to GTiff for Linux

2017-03-04 Per discussione Even Rouault

On samedi 4 mars 2017 18:55:03 CET Micha Silver wrote:
> Hello:
> (I'm double posting as this problem is somewhere between GDAL and QGIS.)
> 
> I have converted an ECW compressed aerial photo to geotiff with JPEG
> compression on a Windows computer with gdal_translate. The resulting tif
> image displays fine with QGIS on the Windows machine.  GDAL version 2.1.3
> on the Windows computer, and QGIS 2.18.4 My command was similar to:
> gdal_translate -of GTiff -projwin . -co TILED=YES -co TFW=YES -co
> COMPRESS=JPEG -co PHOTOMETRIC=YCBCR in.ecw out.tif
> 
> I then copy the tif (and tfw) to my Linux machine, and try to open with
> QGIS, but QGIS crashes immediately. I should note that gdalinfo shows full
> information on both the Windows and Linux systems. On the Linux computer
> GDAL is 2.1.2 and QGIS is 2.18.4 .

On the Linux system where this crashes with QGIS, could you do :

gdalinfo -checksum -stats the.tif

If this crashes, the issue is clearly in GDAL. Otherwise it is likely in QGIS.

> 
> Any ideas what could be the cause of QGIS crashing with the gdal_translate
> converted tif image? Is there something that needs to be done on windows to
> create a Linux usable geotiff??
> 
> Thanks


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Something to think about GeoPackage and WAL

2017-01-12 Per discussione Even Rouault

On jeudi 12 janvier 2017 10:17:52 CET Jukka Rahkonen wrote:
> Hi,
> 
> I did not manage to send mail to QGIS users, but see this:
> http://gis.stackexchange.com/questions/224188/geopackage-error-is-mounted-an
> d-in-wal-mode-this-combination-is-not-allowed
> 
> Could it mean issues which end users can't understand if GeoPackage is
> left to WAL mode and user tries to use it from network drive or as
> read-only? Is it even possible to change the journal mode in this case
> without moving the db into local disk?

Jukka,

(Adding qgis-user in CC)

WAL will not enabled by QGIS if the GeoPackage is located on a network driver 
(as best as we 
can detect that situation !) since the SQLite doc mentions that WAL isn't safe 
in that mode 
(see https://www.sqlite.org/wal.html)
And I don't think WAL can be enabled either on a read-only file/directory (the 
opening will go 
on even if we can't turn WAL on)

If moving a WAL enabled DB (the .gpkg plus the -wal and -shm) on a network 
share, then this 
should probably work, but not in a reliable way in a concurrent use case. If 
moving on a read-
only location, according to the doc, the database shouldn't be openable. If 
moving only the 
.gpkg file without the -wal and -shm, you'll probably get an outdated version 
of the database, 
or will not be able to open it at all. Not sure.

WAL is turned off by QGIS on layer removal, but QGIS must be the last program 
to have the 
file opened so that operation to be successful. If a GeoPackage is left in WAL 
mode (on a 
local disk), a workaround is to close all connections to it, and open / close 
it again with QGIS.

If there are interoperability problems, people can either define the
OGR_SQLITE_JOURNAL environment variable to DELETE or set the QGIS setting 
"/qgis/
walForSqlite3" (in advanced mode) to false, and this will prevent QGIS from 
enabling WAL on 
opening. The drawback is potential deadlocks in some situations where a reader 
and writer 
would run concurrently.

If the current situation causes more harm than good, we can change the  "/qgis/
walForSqlite3" setting default value to false. Together with enabling WAL, I 
made another 
change in the OGR provider that prevented a reader that has finished from 
iterating over 
features from being left in a ghost active state, but I guess there will be 
real read/write 
concurrent situations were deadlocks will appear.

I'm afraid there's no perfect solution. SQLite is a file-based database. One 
cannot ask for it 
the full power of server-based databases regarding concurrent use. 

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS MBTiles

2016-11-23 Per discussione Even Rouault

> No, the error message you got shouldn't normally come up if 
minzoom and
> maxzoom are available. Can you share your dataset (possibly 
privately) so I
> can investigate?
> 

With the file privately sent by Stefan, both gdalinfo and QGIS work fine 
for me. Do you use the latest QGIS 2.14.X point release that ship with a 
recent GDAL ? There might have been fixes in the MBTiles driver.

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS MBTiles

2016-11-23 Per discussione Even Rouault

On mercredi 23 novembre 2016 14:15:58 CET Stefan Löhr wrote:
> Hi,
> 
> thanks for your fast answer.
> 
> I opened the dataset with "spatialite-gui.exe" and it seems to be
> populated with a lot of tiles.
> The minzoom and maxzoom entries in "metadata" are also present:
> ##
> namevalue
> formatpng
> version1.0.0
> center16.500071,63.322846,4
> namedoete
> tagsasdf
> bounds-20.999853,37.645693,53.95,89.00
> legend[]
> creditsasdf
> maxzoom19
> typeoverlay
> attribution
> minzoom0
> descriptionad
> use_constraintsasdf
> ##
> 
> The coordinates seem to be somewhat strange, but the rest looks ok.

What's strange in the coordinates ?
But indeed the rest should be OK

> The strange thing is:
> It says maxzoom=19, but in the "tiles" table, there are only tiles up to
> zoom-level 17.
> Could that be the problem?

No, the error message you got shouldn't normally come up if minzoom and maxzoom 
are 
available. Can you share your dataset (possibly privately) so I can 
investigate? 

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS MBTiles

2016-11-23 Per discussione Even Rouault

Stefan

> 
> when I try to add a raster layer based on a local MBTiles file in QGis
> 2.14 I get this error:
> "Nicht unterstütztes Datenformat: D:\Files\Testdata.mbtiles ist keine
> unterstützte Rasterdatenquelle Cannot find min and max tile numbers"
> 
> Did someone successfully use MBTiles with QGis?

Yes sure.

This error messages comes from the underlying GDAL MBTiles driver. There is 
something in 
your dataset that it doesn't like. Either the driver isn't robust to this case, 
either the mbtiles 
isn't conformant. But hard to tell without access to the file. 

> 
> Where is QGIS searching for the "min and max tile numbers"? Does it have
> to be in one of the metadata tables?
> Or is it somehow calcluated on the fly?

It looks at the metadata table for 'minzoom' and 'maxzoom' entries. If not 
found, it then 
looks at the tiles table to query the min and max values of zoom_level.
The error could occur if the MBTiles dataset is completely empty. Otherwise I'm 
not sure 
without inspecting it.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGis and Mac hardware

2016-11-07 Per discussione Even Rouault

Le lundi 07 novembre 2016 19:18:33, Nicolas Cadieux a écrit :
> Hi,
> 
> Currently GPU acceleration is not used.  Maybe some plugins can use it but
> I have not seen anything yet.  I have put in a request to use GPU for the
> qgis 3.0 raster calculator as it needs to be upgraded for huge rasters. 
> Currently , the qgis raster calculator can not handle big rasters.  The
> new version needs to break up the raster into smaller chucks and if you
> have done that, you might as well use the multi CPU ou GPU to do the rest.
>  I an hopping this will be fixed for 3.0.  Rendering is multithreaded and
> will be faster the more cores you have.  Personally, I would go for the
> quad core without even thinking for this and other softwares.  I would
> think of getting an ssd also.  I have not used qgis on a MAC (PC and
> ubuntu only) but I expect this comment  to be valid for Macs also.

If you use gdalwarp, with GDAL configured with OpenCL (which is rarely the case 
for most standard binary builds), you can get GPU acceleration for non nearest 
neighbour resampling. Although currently the gain is generally modest (or even 
slower) w.r.t multi-threaded acceleration, and you may face issues with buggy 
OpenCL drivers. Might be that the OpenCL code path isn't completely optimized 
too as it is still rather experimental.

> 
> Nicolas
> 
> > Le 7 nov. 2016 à 11:12, Matthias Kuhn-2 [via OSGeo.org]
> >  a écrit :
> > 
> > Hi
> > 
> > On 11/07/2016 05:02 PM, Square Cube wrote:
> > >> I am looking to update some hardware and am trying to evaluate what
> > >> performance difference in QGis might be gained from the options
> > >> currently on offer from Apple.  I have two specific questions:
> > >> 
> > >> - Does QGis have GPU acceleration features that leave the graphics
> > >> card to perform certain calculations, instead of the CPU?
> > 
> > No
> > 
> > >> - Does QGis take advantage of multiple cores and / or multithreaded
> > >> rendering?
> > 
> > Yes, and it is expected that it will make even more use of multiple
> > cores in the upcoming QGIS 3.0 version.
> > 
> > Matthias
> > 
> > >> I am trying to understand if it would be worth my while investing in a
> > >> 15" MacBookPro (with quad core CPU and discreet graphics card) over a
> > >> 13" MacBook Pro (with dual core CPU and integrated graphics).
> > >> 
> > >> Very grateful for any advice.
> > > 
> > > ___
> > > Qgis-user mailing list
> > > [hidden email]
> > > List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
> > > Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
> > 
> > ___
> > Qgis-user mailing list
> > [hidden email]
> > List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
> > Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user
> > 
> > If you reply to this email, your message will be added to the discussion
> > below:
> > http://osgeo-org.1560.x6.nabble.com/QGis-and-Mac-hardware-tp5294661p5294
> > 662.html To start a new topic under Quantum GIS - User, email
> > ml-node+s1560n4125267...@n6.nabble.com To unsubscribe from Quantum GIS -
> > User, click here.
> > NAML
> 
> --
> View this message in context:
> http://osgeo-org.1560.x6.nabble.com/QGis-and-Mac-hardware-tp5294661p529468
> 2.html Sent from the Quantum GIS - User mailing list archive at Nabble.com.

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Error when adding new field to GeoPackage

2016-10-31 Per discussione Even Rouault

Le lundi 31 octobre 2016 14:37:49, Guelat Jerome a écrit :
> Hi,
> 
> I get an error when saving edits after adding a new column to an existing
> GeoPackage table (see the ticket: http://hub.qgis.org/issues/15614). This
> happens for most of the available field types, but text columns seem to
> work OK. I can add columns using ogrinfo but it would be nice to be able
> to do that directly in QGIS.
> 
> Thanks for your help!

I've just fixed it.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] DB manager error

2016-10-29 Per discussione Even Rouault

Le vendredi 28 octobre 2016 21:09:26, vous avez écrit :
> Hi, 
> 
> I'm using QGIS 2.99.0-Master rev. b2587b7 on Ubuntu 16.04 64bits.
> 
> When I tried to open a SQL query on a POSTGIS database I've got the
> following error while trying to open the saved connection to the
> database: 
> 
> 2016-10-28T15:58:35   1   Traceback (most recent call last):
> 
> File 
> "/usr/share/qgis/python/plugins/db_manager/db_model.py", 
line
> 440, in rowCount
> 
>   self._refreshIndex(parent, True)
> 
> File 
> "/usr/share/qgis/python/plugins/db_manager/db_model.py", 
line
> 491, in _refreshIndex
> 
>   if item.populate():
> 
> File 
> "/usr/share/qgis/python/plugins/db_manager/db_model.py", 
line
> 166, in populate
> 
>   if not connection.connect():
> 
> File
> "/usr/share/qgis/python/plugins/db_manager/db_plugins/postgis/plugin.py",
> line 99, in connect
> 
>   return self.connectToUri(uri)
> 
> File 
"/usr/share/qgis/python/plugins/db_manager/db_plugins/plugin.py",
> line 108, in connectToUri
> 
>   self.db = self.databasesFactory(self, uri)
> 
> File
> "/usr/share/qgis/python/plugins/db_manager/db_plugins/postgis/plugin.py",
> line 69, in databasesFactory
> 
>   return PGDatabase(connection, uri)
> 
> File
> "/usr/share/qgis/python/plugins/db_manager/db_plugins/postgis/plugin.py",
> line 107, in __init__
> 
>   Database.__init__(self, connection, uri)
> 
> File 
"/usr/share/qgis/python/plugins/db_manager/db_plugins/plugin.py",
> line 217, in __init__
> 
>   self.connector = self.connectorsFactory(uri)
> 
> File
> "/usr/share/qgis/python/plugins/db_manager/db_plugins/postgis/plugin.py",
> line 110, in connectorsFactory
> 
>   return PostGisDBConnector(uri)
> 
> File
> "/usr/share/qgis/python/plugins/db_manager/db_plugins/postgis/connector.py
> ", line 67, in __init__
> 
>   self.connection = 
psycopg2.connect(expandedConnInfo.encode('utf-8'))
> 
> File 
> "/usr/lib/python3/dist-packages/psycopg2/__init__.py", 
line 164,
> in connect
> 
>   conn = _connect(dsn, 
> connection_factory=connection_factory,
> async=async)
> 
>   TypeError: argument 1 must be str, not bytes
> 
> 
> I can load postgis layers from the same connection without problems.
> What can I do?

I see a pull request about that https://github.com/qgis/QGIS/pull/3635, but it 
is still tagged work in progress.

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] csv/csvt import

2016-10-29 Per discussione Even Rouault

Le samedi 29 octobre 2016 10:18:57, Siki Zoltan a écrit :
> Hi,
> 
> I tried to use csv/csvt files to add Delimited Text Layer.
> The widths of the fields are not considered by QGIS when saving into other
> formats. In the layer properties dialog the width and precision of fields
> are zeros after loading csv/csvt. I suppose the width and precision of
> fields should be preserved after reading of csv/csvt files.
> 
> I tested it in QGIS 2.2/2.8/2.14 but I got maximal width for fields
> (depending on the field type, character 254, float 23.15) when I save the
> file as a shapefile.
> 
> My csvt file:
> "String(20)","Real(11.2)","Real(11.2)","Real(7.2)","String(20)","String(20)
> "
> 
> My csv file:
> point_id,e,n,z,pc,pt
> 2002,650769.92,237368.68,104.560,,
> 2004,650730.28,237418.54,104.520,,
> 2005,650706.42,237412.43,104.270,,
> 2006,650707.76,237449.22,104.590,,
> 
> Is it a feature or bug?
> Any workaround?

If you open through the OGR provider (Add vector layer), this will work. Note 
that the OGR CSV driver will warn (if you do ogrinfo for example) that the 
values in the z field do not respect the Real(7.2) formatting

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Split features using Geopackage layer?

2016-10-19 Per discussione Even Rouault

Le mercredi 19 octobre 2016 14:14:40, Chris Buckmaster a écrit :
> Hi
> 
> I have a geopackage polygon layer and I am trying to split a simple feature
> within it.
> 
> When I try to do it, I get an error message 'failed to insert, UNIQUE
> constraint failed'.
> 
> Is this a known issue with splitting features using the geopackage format,
> it could it be how the layer is set up?

No, it is due to the fact that when splitting the feature, the newly created 
feature receives the same 'fid' value (the first column in the attribute table) 
as the other part, and that when creating the feature, the fid, when set, is 
honoured. Hence given the duplication of fid, the failure. Before saving the 
edits, you can manually unset the fid.

Although something should probably be done either in the splitting tool to 
unset the fid column of the new feature. Or perhaps the provider could be made 
more tolerant in the case of duplicated FID to assign a new one.

This could happen for other formats handled by OGR like spatialite (if not 
going throug the QGIS spatialite provider), FileGDB, etc...

Would probably deserve a ticket in the bug tracker.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS/GDAL GML flipped coordinate order with URN/URI srsName

2016-10-11 Per discussione Even Rouault

Le mardi 11 octobre 2016 12:21:32, Ilkka Rinne a écrit :
> Ok, thanks for the quick answer and bug fix!
> 
> Any idea about the planned version of QGIS that would include GDAL 2.1.2?

I'm not sure what Jürgen's plans are, but once GDAL 2.1.2 will be released and 
available in OSGeo4W, I think that should probably affect all QGIS versions 
that currently use GDAL 2.1.X ?

> 
> 
> Ilkka
> 
> On Tue, Oct 11, 2016 at 12:59 PM, Even Rouault <even.roua...@spatialys.com>
> 
> wrote:
> > Le mardi 11 octobre 2016 10:47:09, Ilkka Rinne a écrit :
> > > Hi,
> > > 
> > > We are experiencing an unwanted flipping of coordinate order
> > 
> > interpretation
> > 
> > > when adding GML 3.2 files as new vector layers in QGIS 2.2 Valmiera.
> > > 
> > > The GML data is has been encoded as a wfs:featureCollection with EU
> > 
> > INSPIRE
> > 
> > > ProtectedSite 4.0 features as collection members (XML Schema:
> > > http://inspire.ec.europa.eu/schemas/ps/4.0/ProtectedSites.xsd). When
> > 
> > using
> > 
> > > the URI-type srsName "http://www.opengis.net/def/crs/EPSG/0/3035; the
> > > coordinates seem to be interpreted (by GDAL?) in order (easting,
> > 
> > northing),
> > 
> > > while the official axis order (according to the EPSG-registry) is
> > > (northing, easting). If we change the srsName to the URN-type format
> > > "urn:ogc:def:crs:EPSG::3035" in the GML file, the import seems to use
> > > the correct coordinate order.
> > > 
> > > To comply with the INSPIRE requirements, we have to use the URI format
> > > srsNames, but this is likely to cause a lot of confusion for the QGIS &
> > > GDAL users of the data if the coordinates are interpreted in the wrong
> > > order.
> > > 
> > > I'm fully aware of the massively confusing lat/lon (x/y) coordinate
> > > order interpretation misalignment between the geographers and the
> > > (GIS)
> > 
> > software
> > 
> > > implementors, but was not aware that this issues also shows it's ugly
> > 
> > head
> > 
> > > for other CRSes than CRS 84/EPSG:4326.
> > 
> > It's a mess indeed, reinforced by constant changes in srsName syntaxes.
> > 
> > > Question 1: Is this an intentional feature or a bug?
> > 
> > This was a bug. I've just fixed it per https://trac.osgeo.org/gdal/
> > ticket/6678
> > (in time for the upcoming GDAL 2.1.2)
> > 
> > > Question 2: Is it possible to force GDAL from the QGIS UI to interpret
> > 
> > the
> > 
> > > coordinate order of these GML files in the official order?
> > 
> > No workaround I can think of, apart the manual replacing as you did.
> > 
> > Even
> > 
> > --
> > Spatialys - Geospatial professional services
> > http://www.spatialys.com

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS/GDAL GML flipped coordinate order with URN/URI srsName

2016-10-11 Per discussione Even Rouault

Le mardi 11 octobre 2016 10:47:09, Ilkka Rinne a écrit :
> Hi,
> 
> We are experiencing an unwanted flipping of coordinate order interpretation
> when adding GML 3.2 files as new vector layers in QGIS 2.2 Valmiera.
> 
> The GML data is has been encoded as a wfs:featureCollection with EU INSPIRE
> ProtectedSite 4.0 features as collection members (XML Schema:
> http://inspire.ec.europa.eu/schemas/ps/4.0/ProtectedSites.xsd). When using
> the URI-type srsName "http://www.opengis.net/def/crs/EPSG/0/3035; the
> coordinates seem to be interpreted (by GDAL?) in order (easting, northing),
> while the official axis order (according to the EPSG-registry) is
> (northing, easting). If we change the srsName to the URN-type format
> "urn:ogc:def:crs:EPSG::3035" in the GML file, the import seems to use the
> correct coordinate order.
> 
> To comply with the INSPIRE requirements, we have to use the URI format
> srsNames, but this is likely to cause a lot of confusion for the QGIS &
> GDAL users of the data if the coordinates are interpreted in the wrong
> order.
> 
> I'm fully aware of the massively confusing lat/lon (x/y) coordinate order
> interpretation misalignment between the geographers and the (GIS) software
> implementors, but was not aware that this issues also shows it's ugly head
> for other CRSes than CRS 84/EPSG:4326.

It's a mess indeed, reinforced by constant changes in srsName syntaxes.

> 
> Question 1: Is this an intentional feature or a bug?

This was a bug. I've just fixed it per https://trac.osgeo.org/gdal/ticket/6678 
(in time for the upcoming GDAL 2.1.2)

> 
> Question 2: Is it possible to force GDAL from the QGIS UI to interpret the
> coordinate order of these GML files in the official order?

No workaround I can think of, apart the manual replacing as you did.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Attribut table in 2.16

2016-10-07 Per discussione Even Rouault

Le vendredi 07 octobre 2016 17:27:43, Pathin a écrit :
> Hello
> 
> 
> i have updated from 2.14.3 to 2.16 and got an issue i cant understand:
> 
> 
> I have imported a SQLite database to a vector layer what has worked
> fine. On this layer i have worked with the attribut table and select a
> specific line and zoomed in to the map.
> 
> But now after the update the layer is available by the table, but not
> shown as a map anymore so i cant select a specific line anymore
> 
> 
> downgrade to 2.14.3 "solved" it so i can work again but i would like to
> share to it to the group

Mind sharing your database ? I've just tried with one of mine and that works 
fine when using a Spatialite connection or when opening it throug Add vector 
layer.

> 
> 
> thx

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] I cannot open this JSON file

2016-10-07 Per discussione Even Rouault

Le vendredi 07 octobre 2016 16:56:21, Willem Korevaar a écrit :
> According to th eopen data from teh Duch government this should be a JSON
> file:
> http://opendata.cbs.nl/ODataApi/OData/81884NED
> From :
> https://data.overheid.nl/data/dataset/onveiligheid-buurtkenmerken-2008-2009
> 
> I cannot open it in QGIS with Protocol > url , there is a faillure and QGIS
> says it doesnt recognize it. Anyone?

This is a JSON file, but definitely not a GeoJSON one.

The content is:
{
  
"odata.metadata":"http://opendata.cbs.nl/ODataApi/OData/81884NED/$metadata","value":
[
{
  
"name":"TableInfos","url":"http://opendata.cbs.nl/ODataApi/OData/81884NED/TableInfos;
},{
  
"name":"UntypedDataSet","url":"http://opendata.cbs.nl/ODataApi/OData/81884NED/UntypedDataSet;
},{
  
"name":"TypedDataSet","url":"http://opendata.cbs.nl/ODataApi/OData/81884NED/TypedDataSet;
},{
  
"name":"DataProperties","url":"http://opendata.cbs.nl/ODataApi/OData/81884NED/DataProperties;
},{
  
"name":"Perioden","url":"http://opendata.cbs.nl/ODataApi/OData/81884NED/Perioden;
}
  ]
}

which seems to be pointers to other files. And from a quick look at a couple of 
them, they don't seem to be GeoJSON files either.


> Thanks in advance,
> Willem Korevaar

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] WFS layer list not updated?

2016-10-01 Per discussione Even Rouault

Hi Ari,

> 
> I'm using 2.16.2 and when I connect to a WFS provider the layer list is
> not updated after the first time I press Connect in the dialog. This
> means that when I publish or unpublish layers at the server those
> changes are not shown in the client.
> 
> I can force the update by removing the client cache (rm -rf
> .qgis2/cache/data7).
> 
> I can fix the problem by changing 'false' to 'true' in the sendGet call
> in QgsWFSCapabilities::requestCapabilities (that is the bool
> forceRefresh parameter). However, changing code seems rather drastic, so
> what am I doing wrong?

This is probably a shortcoming of my recent refactoring in the WFS provider 
where I borrowed some things from the WMS provider. Instead of full caching as 
currently done vs no caching as before, perhaps there could be a medium way: 
if capabilities are needed due to pushing the Connect button, then we should 
re-issue a network request. In other cases, we could use the cached version. 
Actually looking at the code, I think this is what the WMS provider does.

Could you file a ticket about that (if there isn't already one) ?

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS won't open Martian .img files

2016-07-10 Per discussione Even Rouault

Le dimanche 10 juillet 2016 20:33:12, Adam Neather a écrit :
> Thanks to Even Rouault and Nicolas Cadieux for the replies. I think we
> can say that I need the GDAL driver to open the Mars data .img files
> with QGIS.
> 
> Next question - how do I install the driver? Thanks to Even and Nicolas
> I've found the download page and have the software, but I'm not sure how
> to use it with QGIS. Do I have to load it as a custom plug in? If so,
> how? If I follow Even correctly, then I can load it at the same time as
> the .img; is that simply a matter of opening both files (the .img and
> the GDAL driver) at the same time, when loading a raster?

You don't have to download anything. The driver is built in base GDAL. Just 
open the .lbl file in QGIS with "Add raster layer" or drag the .lbl and 
that should work.


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Shapefile with file .cpg(codepage)

2016-07-07 Per discussione Even Rouault

Le jeudi 07 juillet 2016 12:16:17, Andrea Peri a écrit :
> Hi,
> I have some shapefiles with an extra file having extension .cpg.
> The .cpg extension is a optional extension for declaration of Code-Page.
> Is this file knowed and supported from QGIS 2.14. ?
> 
> I don't know if the qgis when loading a shapefile is using the ogr or
> instead is using an own shapefile provider, but however I don't know
> if gdal is knowing and using this .cpg file.

Andrea,

OGR does use the .cpg file when present (I'm just updating 
http://gdal.org/drv_shapefile.html since it only mentions the reading of the 
codepage byte in the DBF header, but the .cpg presence overrides that) and 
QGIS uses OGR to read shapefiles.

But, in QGIS, the default behaviour is to make OGR *not* use the encoding 
detected by the OGR shapefile driver and use instead the user defined encoding 
in the GUI. Unless you go to Settings / Options / Data Sources and uncheck the 
"Ignore shapefile encoding declaration", in which case OGR will manage the 
transcoding itself.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Gdal dev 2.1 in OSGeo prompt.

2016-05-04 Per discussione Even Rouault

Le mercredi 04 mai 2016 21:21:32, Nicolas Cadieux a écrit :
> Hi,
> What is the best way of getting the OSGeoW4 prompt to use the gdal dev
> command line tools rather than the current stable distro? Thanks!

https://lists.osgeo.org/pipermail/gdal-dev/2016-April/044043.html

"""There's a nighly build of GDAL in OSGeo4W now.   Invoke gdal-dev-env from 
the osgeo4w shell to setup it's paths."""

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Open geometryless tables from GeoPackage

2016-04-28 Per discussione Even Rouault

Le jeudi 28 avril 2016 07:03:45, Stefan Ziegler a écrit :
> Yes, compiled and running against gdal master.
> 
> Do you just drag 'n' drop the geopackage file into QGIS? Do you mind trying
> this one: http://sogeo.services/tmp/av_2495.gpkg
> 
> I checked it with "ogrinfo -so av_2495.gpkg" and it only shows the tables
> with geometry (gdal 1.11 and gdal master).

Stefan,

OK, I see. This database hasn't been produced by GDAL. GDAL only produces
and recognizes non-spatial tables that conform to 
http://gdal.org/geopackage_aspatial.html ,
which is a GDAL-specific extension that was created to address the support of
such tables : https://trac.osgeo.org/gdal/ticket/5521 ,
http://lists.opengeospatial.org/pipermail/geopackage/2014-June/thread.html#22

You can upgrade your DB to aspatial by doing the following :

$ sqlite3 av_2495.gpkg

INSERT INTO gpkg_contents (table_name, data_type, last_change)
 SELECT name, 'aspatial', strftime('%Y-%m-%dT%H:%M:%fZ', 'now' ) FROM
 sqlite_master WHERE type = 'table' AND name NOT LIKE 'gpkg_%' AND
 name NOT LIKE 'sqlite_%' AND name NOT LIKE 'rtree_%' AND
 name NOT IN (SELECT table_name FROM gpkg_contents);

INSERT INTO gpkg_extensions
  (table_name, column_name, extension_name, definition, scope)
VALUES
  (
NULL,
NULL,
'gdal_aspatial',
'http://gdal.org/geopackage_aspatial.html',
'read-write'
  );

The original point for the aspatial extension was to be compatible with 
requirements from metadata,
but the latest revision of the GPKG spec has demoted metadata support itself as 
an extension.
So perhaps GDAL could be more lax on opening.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Open geometryless tables from GeoPackage

2016-04-27 Per discussione Even Rouault

Le mercredi 27 avril 2016 21:25:04, Stefan Ziegler a écrit :
> Hi
> 
> is there a working approach to open geometryless tables from a geopackage
> file?

Works for me with https://github.com/qgis/QGIS/pull/3031 applied, but there's 
nothing in it specifically to handle that.
Is your QGIS build against GDAL 2.0 ? (1.11 cannot handle GPKG geometry less 
tables)

> 
> When opening a geopackage file with the standard dialog, there are only the
> tables shown with geometries.
> 
> In the db browser the geometryless tables are shown but cannot be added:
> 
> "/home/stefan/Projekte/qgis-ili2pg-workshop/av_2495.gpkg|layername=Einzelob
> jekte_Einzelobjekt is an invalid layer - not loaded"
> 
> best regards
> Stefan

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] precision Identify Features tool

2016-04-13 Per discussione Even Rouault

Le mercredi 13 avril 2016 20:05:27, Paulo van Breugel a écrit :
> I just opened a raster layer (asci layer), with real numbers with 6
> decimals.
> 
> When click with the 'Identify Features' tool on a raster cell it shows many
> more decimals, e.g., 83.32058715820312 for a raster cell that has value
> 83.320587 (as can be easily checked by opening the layer as text file).
> 
> Any idea what I am missing here?

The raster is probably recognized as a Float32 raster (or is exposed as 
Float64), but the Identify tool  displays with the precision of a 64-bit 
floating point number. And due to how floating point computation works, the 
extension from 32 to 64 bit is done with working with the IEEE-754 binary 
representation and not the decimal representation, hence 83.320587f gets 
converted as 83.32058715820312.

For example with Python :
print struct.unpack('f', struct.pack('f', 83.320587))
(83.32058715820312,)

If the raster is exposed as a Float32 (which is likely if it is recognized by 
the AAIGRID driver : http://gdal.org/frmt_various.html#AAIGrid ), then the 
Identify tool could be improved to respect this precision.

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Problem with QGIS 2.14.1 and MapInfo 12.0.3 coexistence and sharing the same data

2016-04-12 Per discussione Even Rouault

Claus,

This is known as https://hub.qgis.org/issues/14378 . I'm currently working on 
a fix for that.

As a workaround (not necessarily practical), you can try setting the .tab, 
.dat, .map and .id files to be read-only in the Explorer

Even

> Hi everybody
> 
> I have encountered a problem with QGIS and MapInfo coexistence.
> I have no problems opening the same MapInfo tab file (native mapinfo 7.8
> format) from 2 pc's running either QGIS 2.14.1 or Mapinfo 12.0.3. at the
> same time.
> 
> BUT
> If I style the Mapinfo layer in QGIS 2.14.1, categorized from values in a
> column, then MapInfo cannot open the table, or using the info button in an
> already opened table and sends an error :
> Translated from Danish is says: Cannot acces the xxx.DAT for reading. Check
> that the file is not opened from another program.
> 
> Closing QGIS I can again access the file from Mapinfo.
> 
> Anybody has a workaround from that!
> 
> It makes sharing the same data impossible between QGIS 2.14.1 and MapInfo
> users.
> 
> NB! No such problems with QGIS 2.8.5 and MapInfo 12.0.3
> 
> Tecnical details
> PC1 Win8.1 64bit, QGIS 2.14.1 64bit,  MapInfo 12.0.3 32bit
> PC2 ditto
> PC3 Win7 32bit - QGIS 2.8.5 32 bit, MapInfo 12.0.3 32bit
> 
> Regards
> 
> Claus Sørensen
> Viborg Kommune,
> Denmark

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] How do I create a geo-referenced .KAP or .BSB file

2016-04-10 Per discussione Even Rouault

Le dimanche 10 avril 2016 15:14:00, justin bonar a écrit :
> Hello,
> 
> I would like to create geo-referenced RNC files (.KAP or .BSB) from image
> taken from earth.nullschool.net.
> 
> Can someone advise how I would do this?

Justin,

Assuming you're not afraid of self compiling and some tiny source code 
modification...

GDAL has code for write support for BSB, but this is disabled in standard 
builds because of a patent claim on write support.

Some background on this in GDAL's frmts/bsb/bsb_read.c 

If you think you're fine with that, it's nevertheless possible to enable write 
support by uncommenting the //#define BSB_CREATE line at beginning of 
frmts/bsb/bsbdataset.cpp and re-compiling.

Note: you need to use the latest checkout of GDAL trunk since I've just tried 
the above and fixed a couple compilation issues. It may also be not completely 
functional, being disabled by default and little used.

There's also http://libbsb.sourceforge.net/ that might perhaps be used. Not 
sure if it handles georeferencing however.

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Broken Wiki ArcGIS REST example?

2016-04-06 Per discussione Even Rouault

Le mercredi 06 avril 2016 17:12:13, Randal Hale a écrit :
> It hasn't worked for me in a while - I think the break started sometime
> around gdal 1.10 - I *think*. For a while (and I haven't checked in
> quite some time) I couldn't get gdal_translate to produce an xml.

Perhaps related to https://trac.osgeo.org/gdal/ticket/6186

You need GDAL 1.11.4, 2.0.2 or 2.1.0beta1


> 
> 
> 
> Randy
> 
> On 04/06/2016 10:52 AM, Anita Graser wrote:
> > Hi,
> > 
> > I just wanted to try the example listed on
> > https://hub.qgis.org/wiki/17/Arcgis_rest:
> > 
> > qgis.utils.iface.addRasterLayer("http://server.arcgisonline.com/ArcGIS/re
> > st/services/World_Imagery/MapServer?f=json=true","raster")
> > 
> > 
> > ... but it does not seem to work anymore: "GDAL provider cannot open
> > GDAL dataset ..."
> > 
> > Does anyone have a working example?
> > 
> > Best wishes,
> > Anita
> > 
> > 
> > ___
> > Qgis-user mailing list
> > Qgis-user@lists.osgeo.org
> > List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
> > Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] FW: Inconsistent WFS layer requests

2016-03-09 Per discussione Even Rouault

Le mercredi 09 mars 2016 09:22:57, Thomas Colley a écrit :
> Hi
> 
> I mistakenly posted this to developer list instead of user yesterday, hope
> it's ok to post here too.
> 
> I am trying to troubleshoot a problem we're having with a WFS service in
> QGIS 2.14. I am loading in a layer with 'Only request features overlapping
> the current view extent' and have tried cached and uncached. The layer
> loads only 1 feature and that is not within the map canvas. I have used
> Wireshark to capture the requests and can see that QGIS is making the
> GetFeature request with no bbox parameter and a maxfeatures parameter set
> to 1.
> 
> In QGIS 2.12.2 the GetFeature request with no bbox and maxfeatures=1 is
> still made but then it is followed up with another GetFeature request with
> the correct bbox and all features are loaded onto the layer as expected. I
> am wondering if the maxfeatures=1 request is a probe to test that valid
> features are being returned before making the proper request but in 2.14
> the proper request never gets sent. Looking at the layer metadata in QGIS
> the layer source shows the correct request URL.
> 
> It's confusing as we have other WFS services which are working fine in
> 2.14. If I can narrow the issue down to being a problem with the service
> itself I can approach the supplier.

The request with maxfeatures=1 is done when the response to the 
DescribeFeatureType request couldn't determine if the geometry type is point, 
line or polygon. So it downloads a single feature and use its geometry type as 
the layer type.

So I guess there's a logic problem in the code in that probably rarely used 
code path.

Did you try doing a zoom in or zoom out or pan operation ? This might perhaps 
help triggering a new request (not sure though)

> 
> It's a commercial WFS but if anyone is able to test I can provide the
> connection URL.

I'm currently doing a major overhaul of the WFS provider ( 
https://github.com/qgis/QGIS-Enhancement-Proposals/issues/53 ) . That would be 
good if I could test your use case at some point. But not before a few weeks 
as the code isn't ready yet.

> 
> Below are details of the packets being sent from QGIS.
> 
> Thanks very much
> 
> Tom
> 
> WFS connection string:
> http://www.getmapping.com/web.feature.service/FeatureService.svc//x
> /mastermap?
> 
> Layer source from Layer Properties > Metadata tab > Properties
> http://www.getmapping.com/web.feature.service/FeatureService.svc//x
> /mastermap?SERVICE=WFS=1.0.0=GetFeature=topographi
> carea=EPSG:27700=274257.18611745559610426,195451.6365115281660
> 1098,274325.73425826861057431,195546.70830481656594202
> 
> Requests sent when adding WFS layer (caught by Wireshark):
> http://www.getmapping.com/web.feature.service/FeatureService.svc//x
> /mastermap?SERVICE=WFS=1.0.0=topographicarea=274257.1
> 8611745559610426,195451.63651152816601098,274325.73425826861057431,195546.7
> 0830481656594202=DescribeFeatureType
> 
> http://www.getmapping.com/web.feature.service/FeatureService.svc//x
> /mastermap?SERVICE=WFS=1.0.0=GetFeature=topographi
> carea=EPSG:27700&=1
> 
> http://www.getmapping.com/web.feature.service/FeatureService.svc//x
> /mastermap?SERVICE=WFS=1.0.0=GetCapabilities=topog
> raphicarea=EPSG:27700=274257.18611745559610426,195451.63651152
> 816601098,274325.73425826861057431,195546.70830481656594202

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] GDAL/OGR 2.0 in OSGeo4W

2016-02-20 Per discussione Even Rouault

Le dimanche 21 février 2016 00:25:57, Nathan Woodrow a écrit :
> Hey Jürgen,
> 
> Thanks very much.  Interested to see if the TAB editing is supported now.

TAB editing in GDAL was developed () in particular for QGIS use case.

> 
> Regards,
> 
> On Sun, Feb 21, 2016 at 2:50 AM, Jürgen E.  wrote:
> > Hi,
> > 
> > 2.0.2 is now the default GDAL in OSGeo4W and the nightlies (qgis-ltr-dev,
> > qgis-rel-dev and qgis-dev) already picked it up.
> > 
> > With the next release the regular packages (2.14 and 2.8) will also be
> > updated
> > to use it.
> > 
> > Please test and report any issues.
> > 
> > 
> > Jürgen
> > 
> > --
> > Jürgen E. Fischer   norBIT GmbH Tel.
> > +49-4931-918175-31
> > Dipl.-Inf. (FH) Rheinstraße 13  Fax.
> > +49-4931-918175-50
> > Software Engineer   D-26506 Norden
> > http://www.norbit.de
> > QGIS release manager (PSC)  GermanyIRC: jef on
> > FreeNode
> > 
> > ___
> > Qgis-user mailing list
> > Qgis-user@lists.osgeo.org
> > List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
> > Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] OpenJPEG

2016-01-25 Per discussione Even Rouault

Le lundi 25 janvier 2016 17:18:18, Aaron Boxer a écrit :
> Hello,
> 
> I have made improvements to OpenJPEG speed and memory usage,
> and I am releasing the code under the Affero GPL license.

Aaron,

You're of course the ultimate voice for licensing choices of your code but 
given that you've already submitted code as PR to the main OpenJPEG project, 
peple could reasonably assume it was intended to be under the same BSD 2-
clause license. 
Especially since https://github.com/CodecCentral/ronin mentions BSD-2 in the 
main README and LICENSE files.
I find it not very respectful for people included myself who have, apparently 
wrongly, tested or reviewed it under that assumption.
I don't think it will be mergable into mainline openjpeg if you keep it under 
Affero GPL, but that's perhaps intended.

Others: this relates to discussions on gdal-dev (& openjpeg list/tracker).  
Look for openjpeg in subjects at
https://lists.osgeo.org/pipermail/gdal-dev/2016-January/date.html

Even

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS big fail in deleting Shapes!!

2016-01-20 Per discussione Even Rouault

Le mercredi 20 janvier 2016 09:04:52, Neumann, Andreas a écrit :
> Hi,
> 
> Regarding GDAL update to version 2.0: if your platform is Windows, I
> already asked Jürgen if he can update to version 2.0 - as there are
> other issues/new functionality that would be enabled by moving to GDAL
> 2.0 (like better Geopackage support, curve support (for Interlis and
> later also Geopackage) and probably some other improvements). So there
> is a good chance that QGIS 2.14 on Windows will ship with GDAL 2.0.

FYI, GDAL 2.0.2 RC1 should be out for the end of the week.

> 
> @Uwe: while I totally agree that this is a severe error that should be
> fixed, apparently it isn't that easy. And I can only recommend that you
> should not build your GIS projects around Shapefiles. There are other
> options around, like Postgis, SpatiaLite or Geopackage that also allows
> to share data with ArcGIS. And you could also store and edit your data
> in SpatiaLite, Postgis or Geopackage and then only export to Shapefiles
> if you need to deliver your data to your client. Shapefiles are from the
> past century and should be avoided - for many reasons.
> 
> I also remember that Marco Hugentobler (Sourcepole) was contracted to
> fix these Shapefile issues for the province of Vorarlberg - though I may
> be wrong. So there is already someone investing to fix this bug. Marco -
> correct me if I am wrong.
> 
> Hopefully, there can be a fix about this Shapefile issue in QGIS 2.14
> ...
> 
> Andreas
> 
> On 2016-01-20 07:34, Matthias Kuhn wrote:
> > Hi all,
> > 
> > Here is some more eplanation of the issue:
> > 
> > http://gis.stackexchange.com/questions/58101/how-to-avoid-creating-corrup
> > t-shapefiles-during-editing/79677#79677
> > 
> > I spent quite some time to fix this issue as much as possible from QGIS
> > side (Hint: "they" are not just lazy and ignore problems ;-) ) It seems
> > that there is also a second (closely related) problem left to be solved
> > and it looks that it can be solved by shipping a newer GDAL version.
> > Apparently the issue has not been completely solved yet and apparently
> > it is considered to be from "critical" over "big fail" to "fatal" from
> > the audience. It would be awesome if you could ask the maintainer of
> > your operation system's builds for a quote to ship QGIS with GDAL 2.0 on
> > your operation system. We will gladly provide you with email addresses.
> > 
> > Best regards,
> > Matthias
> > 
> > On 01/19/2016 10:54 PM, Mats Elfström wrote:
> > Hi!
> > While all that Nicolas and Uwe says is generally true and reasonable, it
> > is also true that bugs and errors can be more or less dangerous. With
> > limited resources, it is even more important that bugs are handled
> > according to a simple risk analysis. In this case, it would go something
> > like this:
> > What is the probability that users will edit shape files (high, moderate,
> > low): High What is the probability that QGIS will damage shape files:
> > Moderate (as I understand) What is the severity of the potential damage
> > (fatal, reparable, negligible): Fatal And High * Moderate * Fatal equals
> > a probable risk for fatal damage, which ought to prompt for a rapid fix
> > of ths bug, regardless of where it resides in the QGIS ecosystem. But
> > oddly enough, I have browsed through many of the posts on the shape file
> > issue, but have yet to find a technical description of the exact cause
> > and origin of this bug.
> > 
> > Regards, Mats.E
> > 
> > 2016-01-19 22:14 GMT+01:00 Saber Razmjooei
> > :
> > 
> > Dear Uwe,
> > 
> > There are some work-arounds and solutions provided within the ticket.
> > 
> > Also, work is underway to resolve the issue in the next LTR release. Last
> > time I checked (5 minutes ago) it is still present in master.
> > 
> > I totally agree with you that there should be a warning to users to make
> > them aware of this limitation.
> > 
> > In general, QGIS and other GIS/non-GIS proprietary software come with
> > bugs. The reason developers haven't fixed this bug and other bugs is due
> > to lack of resources (time, money).
> > 
> > The project steering committee is committed to resolving the bugs and
> > they have already allocated a big proportion of their income for bug
> > fixing. You can see last year's financial report hear:
> > 
> > http://qgis.org/en/_downloads/PublicQGISfinancialreport2014.pdf
> > 
> > If you think, there should be more work done to increase the quality of
> > QGIS, you can get involved or contribute to the project:
> > 
> > http://qgis.org/en/site/getinvolved/governance/sponsorship/sponsorship.ht
> > ml
> > 
> > Regards,
> > 
> > Saber
> > 
> > FROM: Qgis-user [mailto:qgis-user-boun...@lists.osgeo.org] ON BEHALF OF
> > Uwe Fischer SENT: 19 January 2016 19:10
> > TO: 'qgis-user ML'
> > SUBJECT: [Qgis-user] QGIS big fail in deleting Shapes!!
> > 
> > Hello all,
> > 
> > today I was working on a shape dataset and I had to delete some features.
> >

Re: [Qgis-user] QGIS big fail in deleting Shapes!!

2016-01-20 Per discussione Even Rouault

Le mercredi 20 janvier 2016 11:14:49, Jürgen E. Fischer a écrit :
> Hi Andreas,
> 
> On Wed, 20. Jan 2016 at 09:04:52 +0100, Neumann, Andreas wrote:
> > I also remember that Marco Hugentobler (Sourcepole) was contracted to
> > fix these Shapefile issues for the province of Vorarlberg - though I may
> > be wrong. So there is already someone investing to fix this bug. Marco -
> > correct me if I am wrong.
> 
> Yes, please.  I thought this was a fact - and didn't look at all into
> the problem.
> 
> Not sure if QGIS always does a REPACK (it used to after each delete - but
> that changes feature ids and interferes with internal stuff).

I've had a look at current master : it sets a "mShapefileMayBeCorrupted" flag 
when SetFeature() or DeleteFeature() are called, and if that flag is set at 
QgsOgrProvider::syncToDisc() time, it calls repack(). So it looks OK .

Just tested with QGIS master that deleting a feature leads to repacking when 
saving the edition. That said, that should work with GDAL 1.11 as well. The 
GDAL 2.0 fix is about also make REPACK repacking shapefiles whose geometries 
have been edited.

> I'm also
> not sure that GDAL 2 fixes the problem - the comments to the ticket are
> not conclusive in that regard 

Agreed...

> (GRASS is probably not a good test as it
> also uses GDAL and hence also should follow it's interpretation of what's
> deleted; see also
> https://lists.osgeo.org/pipermail/gdal-dev/2016-January/043483.html).
> 
> 
> Jürgen
> 
> 
> PS: @Saber may your GDAL 2.0.0dev come from experiments with this?  OSGeo4W
> has 1.11.3 (or 2.0.0 as experimental - which QGIS doesn't use - but no
> *dev).

-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Qgis 2.13 on windows 7: ogr2ogr oracle to shape does not keep decimals.

2016-01-11 Per discussione Even Rouault

Le lundi 11 janvier 2016 11:29:59, Lopez Agostini, Diana a écrit :
> Hi,
> 
> We've got a vector layer in a spatial oracle repository.
> We need to export it (both geometry and attributes) to a shape file.
> 
> We cannot use "save as shape" directly from qgis, because the oracle layer
> has a blob field that prevent this operation. And we cannot drop the blob
> field because it is needed for other purposes.
> 
> So we're using OGR2OGR instead.
> The oracle layer have several fields  tipus 'double':
> 
> [cid:image001.png@01D14C63.67A05890]
> 
> After de conversion, the shape file keeps  the double format, but we've
> lost the decimal data. For exemple:
> 
> Before conversion: AREA_HA = 3379.88378172
> After conversion:AREA_HA=  3379
> 
> We have repeated the conversion applying different "locale" configurations,
> in case the problem was originated by the decimal separation character
> used: "." or ",". However, the results are equally bad.
> 
> Do you know any workaround for this issue ?

From http://gdal.org/drv_oci.html :

"""It might be necessary to define the environment variable NLS_LANG to 
"American_America.UTF8" to avoid issues with floating point numbers being 
truncated to integer on non-English environments."""


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user

1 2 >

1 - 100 di 153 matches

Mail list logo