Re: [Qgis-user] Shapefiles

[Nicolas Cadieux]

I agree,

Unless someone actually when true the trouble of story the geometry in the 
database.  Look for x an y columns for point or a geometry column or a WKT 
format.  If so, you could recreate the geometry from that information.  But 
chances are you have no such data in the databas.

Nicolas Cadieux
https://gitlab.com/njacadieux

> Le 8 avr. 2021 à 08:14, Andreas Neumann  a écrit :
> 
> 
> Hi Paddy,
> 
> The .shp file contains the geometries. If your data contains no geometries, 
> then you might be fine with the .dbf file that contains the no-geometry part 
> of your data (attribute columns). No - it is not possible to recreate the 
> .shp file if it is missing.
> 
> Andreas
> 
>> On 2021-04-08 14:07, Pat Brown wrote:
>> 
>> Hi,
>>So, just to confirm what I already suspect? If you have the .shx and .dbf 
>> files but the .shp file is missing then you might as well get rid of them? 
>> It is not possible to recreate the .shp file?
>>  
>> Thanks,
>>  
>> Paddy 
>> 
>> ___
>> 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
___
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] Something to use found - RE: How to "Affine Transform" and rotate at a given point

[Jeff Sonnentag]

Something to use was found based on this:

https://gis.stackexchange.com/questions/322262/rotating-shape-around-defined-point-instead-of-centroid-in-qgis

It's handy to know you can hold CTRL and click to set the point rotation occurs 
around.


From: Qgis-user  On Behalf Of Jeff Sonnentag
Sent: Friday, April 9, 2021 10:21 AM
To: qgis-user 
Subject: [Qgis-user] How to "Affine Transform" and rotate at a given point

So basically we have some "engineering" somebody just made up their own 
coordinate system for that arrived as a DWG and was converted to a DXF in order 
to get it to appear in QGIS (3.16.5).  At first I thought it was in feet but 
aligning a corner showed that wasn't right and it wasn't in inches either.

After figuring out and setting a Scale Factor and then re-aligning to a corner 
things are OK except for a rotation that needs to occur.  It needs to occur 
around the corner match point.  Is that even possible in QGIS?  It is all built 
into the transform function/tool in ArcMap so all you have to do is set where X 
and Y are supposed to come from and their final destination and the rotation 
occurs around that reference point.  Is there a way to do that in QGIS?

Setting a rotation "around the z-axis" sends things off into the wild blue 
yonder. . . .
___
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] How to "Affine Transform" and rotate at a given point

[Jeff Sonnentag]

So basically we have some "engineering" somebody just made up their own 
coordinate system for that arrived as a DWG and was converted to a DXF in order 
to get it to appear in QGIS (3.16.5).  At first I thought it was in feet but 
aligning a corner showed that wasn't right and it wasn't in inches either.

After figuring out and setting a Scale Factor and then re-aligning to a corner 
things are OK except for a rotation that needs to occur.  It needs to occur 
around the corner match point.  Is that even possible in QGIS?  It is all built 
into the transform function/tool in ArcMap so all you have to do is set where X 
and Y are supposed to come from and their final destination and the rotation 
occurs around that reference point.  Is there a way to do that in QGIS?

Setting a rotation "around the z-axis" sends things off into the wild blue 
yonder. . . .
___
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] Generate the bearing in the attribute table

[C Hamilton]

There is a somewhat obscure feature to the Shape Tools plugin
 that will calculate the
distance of each segment in a line and it also calculates the azimuth as
well, but you will not automatically see it as it is in the attribute
table. I don't know if it will work for you or not.

Once you install the Shape Tools plugin, run this algorithm from the QGIS
menu "Vector->Shape Tools->Geodesic measurement layer". Make sure "Measure
total length rather than each line segment" is NOT checked then run it. It
will display measurements of each segment in your lines, but if you open up
the attribute table, you will see a "heading_to" attribute which is your
azimuth. You could right-mouse click on your new layer and change the label
to the heading_to attribute to have it displayed on the map. I probably
should add a feature to that algorithm that allows you to display either
the distance or azimuth.

Regards,

Calvin



On Fri, Apr 9, 2021 at 8:22 AM carlos h  wrote:

> Hello, I think I made a conceptual error, it is not a polygon, but a
> polyline that has "n" number of vertices.
>
> https://ibb.co/vsGF22W
>
> I think that when you see the attached image you will understand what I
> want, what I want is to generate the bearing in the attribute table.
>
> Thanks
>
> CA
>
>
>
> El jue, 8 abr 2021 a las 22:17, chris hermansen ()
> escribió:
>
>> David & list,
>>
>> On Thu, Apr 8, 2021, 20:02 David Strip  wrote:
>>
>>> On 4/8/2021 6:32 PM, chris hermansen wrote:
>>>
>>> Please correct me if I am wrong but are not the start and end points of
>>> a polygon the same?
>>>
>>> I believe he's referring to the start and end point of each edge.
>>>
>>
>> As far as I know a polygon only has one edge, unless it has interior
>> holes. That's what the $geometry would be, anyway.
>>
>>> ___
> 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

Re: [Qgis-user] Generate the bearing in the attribute table

[Kirk Schmidt]

Hi Carlos;

It seems you have the calculation of azimuth figured out.  I would 
generate a new filed in the attribute table and use the field calculator 
to calculate your bearings. If you are comfortable with python, you can 
create a script as well.


If you are using a projected coordinate system, take note that North 
would be actually be grid north and your calculate azimuth would be in 
relation to it.  If you want to reference actual true north, you would 
have to calculate the convergence between true and grid north for your 
project area.  If you want the coordinates referencing magnetic north 
(perhaps for field work), then you would also have to calculate the 
magnetic declination for your project.


    Kirk Schmidt

On 4/8/2021 6:30 PM, Sebastian Gutwein wrote:

Carlos-
Interesting problem-
You might get better responses on stack exchange.
I have the same confusion as Ujaval. Each side of your polygon will 
have a different bearing.
If all you care about is the average angle then you can use this 
expression  main_angle( $geometry)

This will return azimuth of the polygon in degrees from north.
I don't know how to convert this to bearing in QGIS.



On Thu, Apr 8, 2021 at 3:28 PM carlos h > wrote:



  Hello, for example have this in the attribute table:

https://ibb.co/wWn4PfM 

What I need is to have the bearing in the polygon's attribute table.

  Thanks

CA




image.gif
El jue, 8 abr 2021 a las 10:49, Ujaval Gandhi
(mailto:uja...@spatialthoughts.com>>) escribió:

Can you elaborate with an example? It's not clear what the
bearing would be for a polygon.
Logo    
Ujaval Gandhi
Spatial Thoughts
mobile: +91-8095684687
email: uja...@spatialthoughts.com

LinkedIn icon
 Twitter
icon 





On Thu, Apr 8, 2021 at 8:50 PM carlos h
mailto:carloslda...@gmail.com>>
wrote:


Hello, nice to write here.

How can I generate the bearing in the attribute table
of a polygon? the bearing not in azimuth, but for
example N30E etc

Thanks

CA


image.gif
___
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



___
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


--
Kirk Schmidt, MScF, BScF, RPF
General Manager
Nortek Resource Solutions Inc.
RR # 1
Thorburn, NS
B0K 1W0
Tel (902) 922.3607
Email: k...@nortekresources.com
Web: www.nortekresources.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-it-user] recuperare uno shapefile poligonale corrotto

[Andrea Giudiceandrea]

Alessio Degioannini wrote
>> - il plugin core "Geometry Checker"
> In pratica, più tenti di correggerli, più gli 
> errori aumentano: effettivamente il plugin non pare funzionare come 
> dovrebbe. 


Ciao Alessio,
come ti accennavo il plugin "Geometry Checker" purtroppo attualmente
presenta vari bug relativi alle funzionalità di soluzioni dei problemi
geometrici/topologici soprattutto per quanto riguarda le sovrapposizioni e i
gap, come puoi notare anche leggendo le varie segnalazioni di bug in
proposito. In particolare, come ti dicevo, le sovrapposizioni non vengono
corrette. https://github.com/qgis/QGIS/issues/29556

Inoltre, per esempio, la tabella riassuntiva degli errori "corretti" e dei
nuovi errori trovati è spesso totalmente errata.
https://github.com/qgis/QGIS/issues/34664

Non noto problemi invece per le funzionalità di controllo dei problemi
geometrici/topologici.

Ho verificato lo shapefile poligoni.shp che mi hai inviato con 756 feature.
Tutti e tre i tool in QGIS non riportano alcun errore nelle geometrie. Ho
effettuato la verifica anche con ArcMap 9.3.1 e anch'esso non riporta alcun
errore geometrico.

Non ho effettuato verifiche né tramite gli appositi tool dei geo database
PostGIS/PostgreSQL e Spatialite/SQLite, né tramite quelli eventualmente
esistenti in GRASS o SAGA.

Comunque, per riparare errori delle geometrie, il tool di processing QGIS
"Fix geometries" funziona abbastanza bene ed è progettato per usare lo
stesso algoritmo di soluzione degli errori geometrici usato da PostGIS.


Per quanto riguarda invece la presenza di gap o overlaps, effettivamente lo
shapefile ne contiene svariati (circa un migliaio tra gaps e overlaps
secondo Topology Checker, circa 500 secondo Geometry Checker con tolleranza
impostata di default a 1E-8), ma, da quello che ho potuto verificare
controllando il valore delle coordinate memorizzato nello shapefile, quasi
tutti sono dovuti al fatto che vari vertici di poligoni adiacenti
differiscono di pochissimo fra loro.

Per farti un esempio, due punti corrispondenti di due poligoni adiacenti
presentano lo stesso identico valore della coordinata Y mentre i valori
della coordinata X sono:
1: 409478,23033759469...
2: 409478,23033759475,,,
e quindi differiscono di circa 0,006 metri = 0,06 nanometri, che è
circa l'ordine di grandezza delle dimensioni degli atomi.

Non so come sono state create queste piccolissime discrepanze, ma in
generale tali livelli di precisione delle coordinate in ambito GIS è
solitamente eccessivo e inutile.

Solitamente è meglio fare in modo che le coordinate dei vertici delle
geometrie siano memorizzate con una precisione accettabile, per esempio
0,001 metri o 0,005 metri, di fatto facendo in modo che ogni vertice possa
ricadere solo nei nodi di una griglia di 0,001 x 0,001 metri o 0,005 x 0,005
metri.

Per fare questo puoi usare il tool di processing QGIS "Snap points to grid"
impostando per esempio il valore di 0,005 per X e 0,005 per Y.

Se effettui a questo punto i controlli di validità, vedrai che il layer
risultante non avrà errori di geometria, come il layer di partenza, ma
verranno riscontrati solo 72 gaps e 0 overlaps sia dal Topology Checker che
dal Geometry Checker (lasciando a 0 il valore di soglia per l'area di gaps e
overlaps).

Dalla tabella dei risultati del controllo effettuato da Geometry Checker,
ordinando l'elenco degli errori per la colonna Value vedrai che ci sono solo
due gaps con superficie inferiore a 0,01 metri quadrati mentre gli altri 70
hanno superfici superiori a 8 metri quadrati.

Questo significa che solo i primi due gaps con superficie inferiore a 0,01
metri quadrati potrebbero essere effettivi gaps da correggere manualmente,
mentre gli altri 70 in realtà corrispondono a "buchi" interni ai poligoni.

Infatti solitamente si fissa una valore di area di soglia piccola (a seconda
dei propri dati) ma maggiore di 0 per l'individuazione dei gaps. Fissandola
per esempio a 1 metro quadrato, Geometry Checker elencherà solo due gaps.

Per correggere questi unici due gaps, basta editare il layer e rimuovere in
entrambi i casi il vertice in più.

Per imporre uno "snapping" fisso su griglia anche durante la digitazione o
modifica di geometrie, puoi impostare per ogni layer la precisione scelta
per la memorizzazione delle coordinate con l'opzione "Geometry precision"
nella tab Digitizing delle proprietà del layer. Nella stessa finestra di
impostazioni possono essere impostate altre verifiche geometriche e
topologiche che vengono effettuate automaticamente prima del salvataggio del
layer.

Non riscontrato crash, né "propagazione degli errori" facendo delle prove di
editazione del layer: forse sei stato fuorviato dalla errata tabella
riassuntiva degli errori "corretti" e dei nuovi errori trovati del Geometry
Checker.

Spero di esserti stato utile e spero magari che qualche altro utente più
esperto di PostGIS e Spatialite o GRASS e SAGA possa proporre altre
soluzioni.

Comunque penso che lo snapping alla griglia sia 

Re: [Qgis-user] New Map View / New 3D Map View not opening

[Nicole Stoffels]

Hi Charles,

thank you! That was the problem. With a new userprofile my new map opened!

Best wishes
Nicole

Am Fr., 9. Apr. 2021 um 14:37 Uhr schrieb Charles Dixon-Paver <
char...@kartoza.com>:

> It may be a configuration problem. Can you try to open a 3D view with a
> new QGIS userprofile?
>
> On Fri, 9 Apr 2021 at 14:32, Nicole Stoffels 
> wrote:
>
>> Hi Saber,
>>
>> thank you for your fast reply and the tutorial!
>>
>> Unfortunately I never get to that point. Some colleagues still use QGIS
>> 3.14 and there the new window with the map appears, when they click on View
>> > New 3D Map View. However since my upgrade to QGIS 3.16 no new window
>> opens. In fact there is no reaction at all, when I click View > New 3D Map
>> View.
>>
>> Best regards
>> Nicole
>>
>> Am Fr., 9. Apr. 2021 um 12:26 Uhr schrieb Saber Razmjooei <
>> saber.razmjo...@lutraconsulting.co.uk>:
>>
>>> Hi Nicole,
>>>
>>> You may try this tutorial:
>>>
>>> https://www.lutraconsulting.co.uk/blog/2018/03/01/working-with-qgis-3d-part-1/
>>>
>>> The interface has slightly changed/improved, but in general it gives you
>>> an idea how to view your terrain rasters in 3D.
>>>
>>> Kind regards
>>> Saber
>>>
>>>
>>> On Fri, 9 Apr 2021 at 09:10, Nicole Stoffels 
>>> wrote:
>>>
 Dear all,

 I would like to create a 3D map view of DEM data with QGIS 3.16.5 on a
 Windows 10. I used this dataset for a test:


 https://www.lgln.niedersachsen.de/download/121709/ATKIS-DGM_DGM1_DGM5_DGM10_DGM25_DGM50_aus_dem_aktuellen_Laserdatenbestand_.zip

 I load the data as a raster which works just fine (EPSG: 25832). I can
 produce hillshades and everything. But when I go to View -> New 3D Map View
 absolutely nothing happens. No new window opens. I don't even get an error
 message. The same happens with New Map View.

 At gis.stackexchange I can't find entries with similar problems.

 Does anyone have any ideas, solutions, the same problems?

 Thank you very much in advance!

 Best regards
 Nicole

 ___
 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

>>>
>>>
>>> --
>>> Saber Razmjooei
>>> www.lutraconsulting.co.uk
>>>
>>> ___
>> 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

Re: [Qgis-user] New Map View / New 3D Map View not opening

[Charles Dixon-Paver]

It may be a configuration problem. Can you try to open a 3D view with a new
QGIS userprofile?

On Fri, 9 Apr 2021 at 14:32, Nicole Stoffels 
wrote:

> Hi Saber,
>
> thank you for your fast reply and the tutorial!
>
> Unfortunately I never get to that point. Some colleagues still use QGIS
> 3.14 and there the new window with the map appears, when they click on View
> > New 3D Map View. However since my upgrade to QGIS 3.16 no new window
> opens. In fact there is no reaction at all, when I click View > New 3D Map
> View.
>
> Best regards
> Nicole
>
> Am Fr., 9. Apr. 2021 um 12:26 Uhr schrieb Saber Razmjooei <
> saber.razmjo...@lutraconsulting.co.uk>:
>
>> Hi Nicole,
>>
>> You may try this tutorial:
>>
>> https://www.lutraconsulting.co.uk/blog/2018/03/01/working-with-qgis-3d-part-1/
>>
>> The interface has slightly changed/improved, but in general it gives you
>> an idea how to view your terrain rasters in 3D.
>>
>> Kind regards
>> Saber
>>
>>
>> On Fri, 9 Apr 2021 at 09:10, Nicole Stoffels 
>> wrote:
>>
>>> Dear all,
>>>
>>> I would like to create a 3D map view of DEM data with QGIS 3.16.5 on a
>>> Windows 10. I used this dataset for a test:
>>>
>>>
>>> https://www.lgln.niedersachsen.de/download/121709/ATKIS-DGM_DGM1_DGM5_DGM10_DGM25_DGM50_aus_dem_aktuellen_Laserdatenbestand_.zip
>>>
>>> I load the data as a raster which works just fine (EPSG: 25832). I can
>>> produce hillshades and everything. But when I go to View -> New 3D Map View
>>> absolutely nothing happens. No new window opens. I don't even get an error
>>> message. The same happens with New Map View.
>>>
>>> At gis.stackexchange I can't find entries with similar problems.
>>>
>>> Does anyone have any ideas, solutions, the same problems?
>>>
>>> Thank you very much in advance!
>>>
>>> Best regards
>>> Nicole
>>>
>>> ___
>>> 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
>>>
>>
>>
>> --
>> Saber Razmjooei
>> www.lutraconsulting.co.uk
>>
>> ___
> 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

Re: [Qgis-user] New Map View / New 3D Map View not opening

[Nicole Stoffels]

Hi Saber,

thank you for your fast reply and the tutorial!

Unfortunately I never get to that point. Some colleagues still use QGIS
3.14 and there the new window with the map appears, when they click on View
> New 3D Map View. However since my upgrade to QGIS 3.16 no new window
opens. In fact there is no reaction at all, when I click View > New 3D Map
View.

Best regards
Nicole

Am Fr., 9. Apr. 2021 um 12:26 Uhr schrieb Saber Razmjooei <
saber.razmjo...@lutraconsulting.co.uk>:

> Hi Nicole,
>
> You may try this tutorial:
>
> https://www.lutraconsulting.co.uk/blog/2018/03/01/working-with-qgis-3d-part-1/
>
> The interface has slightly changed/improved, but in general it gives you
> an idea how to view your terrain rasters in 3D.
>
> Kind regards
> Saber
>
>
> On Fri, 9 Apr 2021 at 09:10, Nicole Stoffels 
> wrote:
>
>> Dear all,
>>
>> I would like to create a 3D map view of DEM data with QGIS 3.16.5 on a
>> Windows 10. I used this dataset for a test:
>>
>>
>> https://www.lgln.niedersachsen.de/download/121709/ATKIS-DGM_DGM1_DGM5_DGM10_DGM25_DGM50_aus_dem_aktuellen_Laserdatenbestand_.zip
>>
>> I load the data as a raster which works just fine (EPSG: 25832). I can
>> produce hillshades and everything. But when I go to View -> New 3D Map View
>> absolutely nothing happens. No new window opens. I don't even get an error
>> message. The same happens with New Map View.
>>
>> At gis.stackexchange I can't find entries with similar problems.
>>
>> Does anyone have any ideas, solutions, the same problems?
>>
>> Thank you very much in advance!
>>
>> Best regards
>> Nicole
>>
>> ___
>> 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
>>
>
>
> --
> Saber Razmjooei
> www.lutraconsulting.co.uk
>
>
___
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] Generate the bearing in the attribute table

[carlos h]

Hello, I think I made a conceptual error, it is not a polygon, but a
polyline that has "n" number of vertices.

https://ibb.co/vsGF22W

I think that when you see the attached image you will understand what I
want, what I want is to generate the bearing in the attribute table.

Thanks

CA



El jue, 8 abr 2021 a las 22:17, chris hermansen ()
escribió:

> David & list,
>
> On Thu, Apr 8, 2021, 20:02 David Strip  wrote:
>
>> On 4/8/2021 6:32 PM, chris hermansen wrote:
>>
>> Please correct me if I am wrong but are not the start and end points of a
>> polygon the same?
>>
>> I believe he's referring to the start and end point of each edge.
>>
>
> As far as I know a polygon only has one edge, unless it has interior
> holes. That's what the $geometry would be, anyway.
>
>>
___
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-it-user] recuperare uno shapefile poligonale corrotto

[Alessio Degioannini]

Ciao Andrea,

tutto molto chiaro ed esauriente!

Procedo come mi hai spiegato.

Grazie e cordiali saluti.

Alessio

Il 09/04/2021 11:23, Andrea Giudiceandrea ha scritto:

Alessio Degioannini wrote

- il plugin core "Geometry Checker"

In pratica, più tenti di correggerli, più gli
errori aumentano: effettivamente il plugin non pare funzionare come
dovrebbe.


Ciao Alessio,
come ti accennavo il plugin "Geometry Checker" purtroppo attualmente
presenta vari bug relativi alle funzionalità di soluzioni dei problemi
geometrici/topologici soprattutto per quanto riguarda le sovrapposizioni e i
gap, come puoi notare anche leggendo le varie segnalazioni di bug in
proposito. In particolare, come ti dicevo, le sovrapposizioni non vengono
corrette. https://github.com/qgis/QGIS/issues/29556

Inoltre, per esempio, la tabella riassuntiva degli errori "corretti" e dei
nuovi errori trovati è spesso totalmente errata.
https://github.com/qgis/QGIS/issues/34664

Non noto problemi invece per le funzionalità di controllo dei problemi
geometrici/topologici.

Ho verificato lo shapefile poligoni.shp che mi hai inviato con 756 feature.
Tutti e tre i tool in QGIS non riportano alcun errore nelle geometrie. Ho
effettuato la verifica anche con ArcMap 9.3.1 e anch'esso non riporta alcun
errore geometrico.

Non ho effettuato verifiche né tramite gli appositi tool dei geo database
PostGIS/PostgreSQL e Spatialite/SQLite, né tramite quelli eventualmente
esistenti in GRASS o SAGA.

Comunque, per riparare errori delle geometrie, il tool di processing QGIS
"Fix geometries" funziona abbastanza bene ed è progettato per usare lo
stesso algoritmo di soluzione degli errori geometrici usato da PostGIS.


Per quanto riguarda invece la presenza di gap o overlaps, effettivamente lo
shapefile ne contiene svariati (circa un migliaio tra gaps e overlaps
secondo Topology Checker, circa 500 secondo Geometry Checker con tolleranza
impostata di default a 1E-8), ma, da quello che ho potuto verificare
controllando il valore delle coordinate memorizzato nello shapefile, quasi
tutti sono dovuti al fatto che vari vertici di poligoni adiacenti
differiscono di pochissimo fra loro.

Per farti un esempio, due punti corrispondenti di due poligoni adiacenti
presentano lo stesso identico valore della coordinata Y mentre i valori
della coordinata X sono:
1: 409478,23033759469...
2: 409478,23033759475,,,
e quindi differiscono di circa 0,006 metri = 0,06 nanometri, che è
circa l'ordine di grandezza delle dimensioni degli atomi.

Non so come sono state create queste piccolissime discrepanze, ma in
generale tali livelli di precisione delle coordinate in ambito GIS è
solitamente eccessivo e inutile.

Solitamente è meglio fare in modo che le coordinate dei vertici delle
geometrie siano memorizzate con una precisione accettabile, per esempio
0,001 metri o 0,005 metri, di fatto facendo in modo che ogni vertice possa
ricadere solo nei nodi di una griglia di 0,001 x 0,001 metri o 0,005 x 0,005
metri.

Per fare questo puoi usare il tool di processing QGIS "Snap points to grid"
impostando per esempio il valore di 0,005 per X e 0,005 per Y.

Se effettui a questo punto i controlli di validità, vedrai che il layer
risultante non avrà errori di geometria, come il layer di partenza, ma
verranno riscontrati solo 72 gaps e 0 overlaps sia dal Topology Checker che
dal Geometry Checker (lasciando a 0 il valore di soglia per l'area di gaps e
overlaps).

Dalla tabella dei risultati del controllo effettuato da Geometry Checker,
ordinando l'elenco degli errori per la colonna Value vedrai che ci sono solo
due gaps con superficie inferiore a 0,01 metri quadrati mentre gli altri 70
hanno superfici superiori a 8 metri quadrati.

Questo significa che solo i primi due gaps con superficie inferiore a 0,01
metri quadrati potrebbero essere effettivi gaps da correggere manualmente,
mentre gli altri 70 in realtà corrispondono a "buchi" interni ai poligoni.

Infatti solitamente si fissa una valore di area di soglia piccola (a seconda
dei propri dati) ma maggiore di 0 per l'individuazione dei gaps. Fissandola
per esempio a 1 metro quadrato, Geometry Checker elencherà solo due gaps.

Per correggere questi unici due gaps, basta editare il layer e rimuovere in
entrambi i casi il vertice in più.

Per imporre uno "snapping" fisso su griglia anche durante la digitazione o
modifica di geometrie, puoi impostare per ogni layer la precisione scelta
per la memorizzazione delle coordinate con l'opzione "Geometry precision"
nella tab Digitizing delle proprietà del layer. Nella stessa finestra di
impostazioni possono essere impostate altre verifiche geometriche e
topologiche che vengono effettuate automaticamente prima del salvataggio del
layer.

Non riscontrato crash, né "propagazione degli errori" facendo delle prove di
editazione del layer: forse sei stato fuorviato dalla errata tabella
riassuntiva degli errori "corretti" e dei nuovi errori trovati del Geometry
Checker.

Spero di esserti stato utile e spero 

Re: [Qgis-user] New Map View / New 3D Map View not opening

[Saber Razmjooei]

Hi Nicole,

You may try this tutorial:
https://www.lutraconsulting.co.uk/blog/2018/03/01/working-with-qgis-3d-part-1/

The interface has slightly changed/improved, but in general it gives you an
idea how to view your terrain rasters in 3D.

Kind regards
Saber


On Fri, 9 Apr 2021 at 09:10, Nicole Stoffels 
wrote:

> Dear all,
>
> I would like to create a 3D map view of DEM data with QGIS 3.16.5 on a
> Windows 10. I used this dataset for a test:
>
>
> https://www.lgln.niedersachsen.de/download/121709/ATKIS-DGM_DGM1_DGM5_DGM10_DGM25_DGM50_aus_dem_aktuellen_Laserdatenbestand_.zip
>
> I load the data as a raster which works just fine (EPSG: 25832). I can
> produce hillshades and everything. But when I go to View -> New 3D Map View
> absolutely nothing happens. No new window opens. I don't even get an error
> message. The same happens with New Map View.
>
> At gis.stackexchange I can't find entries with similar problems.
>
> Does anyone have any ideas, solutions, the same problems?
>
> Thank you very much in advance!
>
> Best regards
> Nicole
>
> ___
> 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
>


-- 
Saber Razmjooei
www.lutraconsulting.co.uk
___
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-it-user] recuperare uno shapefile poligonale corrotto

[Stefano Campus]

grazie andrea per la tua dettagliata spiegazione.
ricordo una discussione di qualche anno fa sulla lista GFOSS proprio su
questi problemi legati alle coordinate, a cura di sandro furieri.
quello che spesso succede è che passando da una struttura fisica ad
un'altra o da un software ad un altro, si creiino degli errori dovuti
all'approssimazione delle coordinate.
se no ricordo male, ad esempio, arcgis aggancia alla sua griglia interna e
compie già un arrotondamento interno, così un dataset esportato da arcgis a
postgres o viceversa può dare adito a queste incongruenza.

s.

Il giorno ven 9 apr 2021 alle ore 11:23 Andrea Giudiceandrea <
andreaer...@libero.it> ha scritto:

> Alessio Degioannini wrote
> >> - il plugin core "Geometry Checker"
> > In pratica, più tenti di correggerli, più gli
> > errori aumentano: effettivamente il plugin non pare funzionare come
> > dovrebbe.
>
>
> Ciao Alessio,
> come ti accennavo il plugin "Geometry Checker" purtroppo attualmente
> presenta vari bug relativi alle funzionalità di soluzioni dei problemi
> geometrici/topologici soprattutto per quanto riguarda le sovrapposizioni e
> i
> gap, come puoi notare anche leggendo le varie segnalazioni di bug in
> proposito. In particolare, come ti dicevo, le sovrapposizioni non vengono
> corrette. https://github.com/qgis/QGIS/issues/29556
>
> Inoltre, per esempio, la tabella riassuntiva degli errori "corretti" e dei
> nuovi errori trovati è spesso totalmente errata.
> https://github.com/qgis/QGIS/issues/34664
>
> Non noto problemi invece per le funzionalità di controllo dei problemi
> geometrici/topologici.
>
> Ho verificato lo shapefile poligoni.shp che mi hai inviato con 756 feature.
> Tutti e tre i tool in QGIS non riportano alcun errore nelle geometrie. Ho
> effettuato la verifica anche con ArcMap 9.3.1 e anch'esso non riporta alcun
> errore geometrico.
>
> Non ho effettuato verifiche né tramite gli appositi tool dei geo database
> PostGIS/PostgreSQL e Spatialite/SQLite, né tramite quelli eventualmente
> esistenti in GRASS o SAGA.
>
> Comunque, per riparare errori delle geometrie, il tool di processing QGIS
> "Fix geometries" funziona abbastanza bene ed è progettato per usare lo
> stesso algoritmo di soluzione degli errori geometrici usato da PostGIS.
>
>
> Per quanto riguarda invece la presenza di gap o overlaps, effettivamente lo
> shapefile ne contiene svariati (circa un migliaio tra gaps e overlaps
> secondo Topology Checker, circa 500 secondo Geometry Checker con tolleranza
> impostata di default a 1E-8), ma, da quello che ho potuto verificare
> controllando il valore delle coordinate memorizzato nello shapefile, quasi
> tutti sono dovuti al fatto che vari vertici di poligoni adiacenti
> differiscono di pochissimo fra loro.
>
> Per farti un esempio, due punti corrispondenti di due poligoni adiacenti
> presentano lo stesso identico valore della coordinata Y mentre i valori
> della coordinata X sono:
> 1: 409478,23033759469...
> 2: 409478,23033759475,,,
> e quindi differiscono di circa 0,006 metri = 0,06 nanometri, che è
> circa l'ordine di grandezza delle dimensioni degli atomi.
>
> Non so come sono state create queste piccolissime discrepanze, ma in
> generale tali livelli di precisione delle coordinate in ambito GIS è
> solitamente eccessivo e inutile.
>
> Solitamente è meglio fare in modo che le coordinate dei vertici delle
> geometrie siano memorizzate con una precisione accettabile, per esempio
> 0,001 metri o 0,005 metri, di fatto facendo in modo che ogni vertice possa
> ricadere solo nei nodi di una griglia di 0,001 x 0,001 metri o 0,005 x
> 0,005
> metri.
>
> Per fare questo puoi usare il tool di processing QGIS "Snap points to grid"
> impostando per esempio il valore di 0,005 per X e 0,005 per Y.
>
> Se effettui a questo punto i controlli di validità, vedrai che il layer
> risultante non avrà errori di geometria, come il layer di partenza, ma
> verranno riscontrati solo 72 gaps e 0 overlaps sia dal Topology Checker che
> dal Geometry Checker (lasciando a 0 il valore di soglia per l'area di gaps
> e
> overlaps).
>
> Dalla tabella dei risultati del controllo effettuato da Geometry Checker,
> ordinando l'elenco degli errori per la colonna Value vedrai che ci sono
> solo
> due gaps con superficie inferiore a 0,01 metri quadrati mentre gli altri 70
> hanno superfici superiori a 8 metri quadrati.
>
> Questo significa che solo i primi due gaps con superficie inferiore a 0,01
> metri quadrati potrebbero essere effettivi gaps da correggere manualmente,
> mentre gli altri 70 in realtà corrispondono a "buchi" interni ai poligoni.
>
> Infatti solitamente si fissa una valore di area di soglia piccola (a
> seconda
> dei propri dati) ma maggiore di 0 per l'individuazione dei gaps. Fissandola
> per esempio a 1 metro quadrato, Geometry Checker elencherà solo due gaps.
>
> Per correggere questi unici due gaps, basta editare il layer e rimuovere in
> entrambi i casi il vertice in più.
>
> Per imporre uno "snapping" fisso su 

[Qgis-user] New Map View / New 3D Map View not opening

[Nicole Stoffels]

Dear all,

I would like to create a 3D map view of DEM data with QGIS 3.16.5 on a
Windows 10. I used this dataset for a test:

https://www.lgln.niedersachsen.de/download/121709/ATKIS-DGM_DGM1_DGM5_DGM10_DGM25_DGM50_aus_dem_aktuellen_Laserdatenbestand_.zip

I load the data as a raster which works just fine (EPSG: 25832). I can
produce hillshades and everything. But when I go to View -> New 3D Map View
absolutely nothing happens. No new window opens. I don't even get an error
message. The same happens with New Map View.

At gis.stackexchange I can't find entries with similar problems.

Does anyone have any ideas, solutions, the same problems?

Thank you very much in advance!

Best regards
Nicole
___
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] QGIS Server error message - QNetworkDiskCache::prepare()

[Arnaud Vandecasteele]

Hello,

When I check the logs of my QGIS server I always get the same error:
QNetworkDiskCache::prepare() unable to open temporary file

Could you tell me if this is related to an appropriate configuration?
Does it have any performance consequences?

In previous discussions, one solution was to add the parameter
FcgidInitialEnv QGIS_PREFIX_PATH /usr to the virtualhost. However, in my
case, this did not change anything to the display of this error.

Thanks in advance

-- 

Arnaud Vandecasteele
___
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-it-user] Analisi vettori

[Totò]

FabioP wrote
> Buongiorno,
> vorrei sapere come posso ottenere informazioni relative unicamente alle
> porzioni di geometrie contenute in una determinata area.
> Mi spiego meglio.
> Ho un layer che rappresenta le varie tipologie di strade presenti sul
> territorio (autostrada, strada principale, strada secondaria, ecc.)
> Dovrei sapere quanti metri di strada, per ogni topologia, sono presenti in
> una zona (p.es. un quartiere); la zona è definita con un poligono su un
> altro layer.
> Quali strumenti devo usare?

Buongiorno,
se vuoi le info nello stesso layer (senza creare altri layer) usa il
calcolatore di campi;
se la creazione di un nuovo layer non è un problema usa gli strumenti di
processing;

il tutto è spiegato in questo esempio:
http://hfcqgis.opendatasicilia.it/it/latest/esempi/somma_lunghezze_nel_poligono.html

saluti



-
https://pigrecoinfinito.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html
___
QGIS-it-user mailing list
QGIS-it-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/qgis-it-user

Re: [Qgis-user] WMS Image Projection

[Springfield Harrison]

Hello Andreas & Kirk,

Thanks very much for your suggestions.  I also got help from the 
regional Gov't GIS Manager.  He suggested trying the REST source 
(server?) through the QGIS ArcGIS Map Service 
(https://mapservices.crd.bc.ca/arcgis/rest/services).  It worked and the 
misalignment disappeared!  In QGIS, toggling the two different 2019 
orthophotos showed them jumping by about 1.5 m. Several miles away, the 
result was the same.


The REST source did the trick although the underlying data is the same.

The REST property lines were identical to the WMS version.

Not sure why this works but will use the REST server in the future.

Thanks again for your help . . .

-
Cheers, Spring



On 08/Apr/2021 05:48, kirk wrote:

Hi Springer.
My mistake, you will require a local grid shift file, BC_93_05.gsb. 
The NTv2 is a national NAD27
to NAD83 grid shift file. The BC_93_95 provides localized conversion 
from Nad83 (original) to Nad83 (csrs).



Sent from my Galaxy


 Original message 
From: Springfield Harrison 
Date: 2021-04-08 5:24 a.m. (GMT-04:00)
To: Andreas Neumann 
Cc: qgis-user 
Subject: Re: [Qgis-user] WMS Image Projection

Hi Andreas,

Thanks for the reply, I'll try to clarify:

 1. Version is 3.16.4
 2. The two municipal sources are being viewed in their own web
windows and are displaying the same data
 1. They are the source of the property lines for my GIS
(downloaded shapefile) and the orthophoto (WMS link, not
resident TIFF or other)
 3. The property lines are EPSG:3157 - NAD83(CSRS) / UTM zone 10N -
Projected, and from the municipal source.  I've attached the
shapefile if you have the time.
 4. In QGIS, the ortho is loaded from the WMS source
 1. 
https://mapservices.crd.bc.ca/arcgis/services/OrthoImagery2019/MapServer/WMSserver
 2. EPSG:26910 - NAD83 / UTM zone 10N - Projected
 5. The project CRS is EPSG: 26910 but it doesn't seem to change
anything by changing that or the source projections (at least
within versions of NAD83 UTM Zone 10N)
 6. Thus, the same information displays differently online /versus/ in
QGIS.
 1. The shift seems to be 1.5 m +/-, quite a lot for what I'm
trying to do (building permit with tight setbacks).
 2. These are commonly used projections for our neck of the woods
 3. The two images I sent earlier depict the differences; in QGIS,
earlier orthophotos show a similar or worse displacement. Not
so on the municipal web viewer.
 7. I've read up on projections in the tutorial and Help file, no
revelations there, at least that I could see.

Thanks very much for any thoughts you may have.

-
Cheers, Spring


On 08/Apr/2021 00:47, Andreas Neumann wrote:


Hi,

Are there two different CRS involved here? Are the vectors (parcel 
boundaries) in a different CRS than the orthoimages? This might be a 
problem of a missing grid-shift file in QGIS.


It is also unclear for me if you load your Orthoimages locally (from 
tiff-files as well) or from WMS (or both) and only one of them is 
shifted?


Can you also add information what QGIS version with what proj version 
you use? You can find this information in the about dialogue in the 
help menu.


Greetings,

Andreas

On 2021-04-08 02:37, Springfield Harrison wrote:


Hello All,

Re obtaining orthophotos from a WMS server.

Locally we have orthophotos that can be brought into QGIS using the 
WMS service.  They are also viewable on the local government web map 
service.  In both cases I can compare them to the property lines 
which are from the same source.


On the government web pages (2), the orthophotos and the lot lines 
appear to coincide very well.  In QGIS, the same two layers are 
quite misaligned.


One would suspect a CRS problem but changing the CRS of the Project 
and the WMS image in QGIS has no effect.  Even followed the tutorial 
to no effect.


I'm attaching screen shots of the two results.  Ant thoughts would 
be appreciated, thank you . . . .


-
Cheers, Springfield Harrison



___
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

[QGIS-it-user] Analisi vettori

[FabioP]

Buongiorno,
vorrei sapere come posso ottenere informazioni relative unicamente alle
porzioni di geometrie contenute in una determinata area.
Mi spiego meglio.
Ho un layer che rappresenta le varie tipologie di strade presenti sul
territorio (autostrada, strada principale, strada secondaria, ecc.)
Dovrei sapere quanti metri di strada, per ogni topologia, sono presenti in
una zona (p.es. un quartiere); la zona è definita con un poligono su un
altro layer.
Quali strumenti devo usare?

Grazie mille dell'aiuto
Fabio



--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html
___
QGIS-it-user mailing list
QGIS-it-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/qgis-it-user