Ortofoto a colori ad alta risoluzione !!!!!!

[tb14]

Salve ragazzi, la discussione si sta facendo veramente interessante, ma mi si è presentato un problema nell'utilizzo della routine merge per unire più ecw in un unico file. La ho utilizzato in questo modo:

>gdal_merge.bat -o 564.ecw -of ECW -co ="LARGE_OK=TRUE" 564010.ecw 564020.ecw 564030.ecw 564040.ecw 564050.ecw 564060.ecw 564070.ecw 564080.ecw

ed ho ottenuto il seguente messaggio di errore:

ERROR 6: GDALDriver::Create() ... no create method implemented for this format.
Traceback (most recent call last):
File "C:\PROGRA~1\FWTOOL~1.0A8\bin\gdal_merge.py", line 453, in ?
t_fh.SetGeoTransform( geotransform )
AttributeError: 'NoneType' object has no attribute 'SetGeoTransform'

Ora, non sono ancora andato a spulciare il sorgente in py, in quanto non particolarmente esperto. Qualcuno ha un suggerimento da darmi?

Hai ragione l'ECW ha solo il CreateCopy e non il Create diretto.
Vuol solo dire che devi passare attraverso un altro formato ( per esempio GTiff )

Prova :

> gdal_merge.bat -o 564.tif -of GTiff -co ="LARGE_OK=TRUE" 564010.ecw 564020.ecw 564030.ecw 564040.ecw 564050.ecw 564060.ecw 564070.ecw 564080.ecw

e poi trasforma il tiff in ecw :

> gdal_translate -of ECW co ="LARGE_OK=TRUE" 564.tif 564.ecw

Fammi sapere

Prodotti MTB scontati

 

[Excalib65]

Il merge sembra funzionare perfettamene (sta elaborando ed i tempi non sembrano brevi!!!)

Il problema ora è un altro: la dimensione del file che sto ottenendo è di ben 35250 x 11380 pixel, che credo sia superiore al limite imposto dal compressore ECW (500Mb), quindi mi vedo costretto ad utilizzare un formato definitivo diverso dall'ECW (ad es. il formato ozfx3 di OziExplorer, visto che utilizzo OziExplorer).

Qualcuno ha risolto il problema in altro modo? E' possibile comunque comprimere files più grossi di 500 Mb in ECW?

 

[tb14]

Il merge sembra funzionare perfettamene (sta elaborando ed i tempi non sembrano brevi!!!)

Il problema ora è un altro: la dimensione del file che sto ottenendo è di ben 35250 x 11380 pixel, che credo sia superiore al limite imposto dal compressore ECW (500Mb), quindi mi vedo costretto ad utilizzare un formato definitivo diverso dall'ECW (ad es. il formato ozfx3 di OziExplorer, visto che utilizzo OziExplorer).

Qualcuno ha risolto il problema in altro modo? E' possibile comunque comprimere files più grossi di 500 Mb in ECW?

Anche a questo pensa il grande Frank War-----m ( FWTools ) e la sua Gdal. ( sono stato censurato
L'opzione -co ="LARGE_OK=TRUE" ( -co stà per create-option ) permette di comprimere qualunque dimensione di file. Resta poi a l'utilizzatore assicurarsi di avere la licenza per poterlo fare. siccome stiamo utilizzando la compressione solo per hobby non credo che questo sia più di tanto lesivo dei diritti di autore .

TB14

 

[Excalib65]

Grazie del suggerimento tb14, ma non ottengo il risultato sperato:

C:\Programmi\FWTools1.0.0a8>gdal_translate -of ecw -co="LARGE_OK=TRUE" S_Antioco.tif S_Antioco.ecw
Input file size is 10975, 18341
0.ERROR 1: ECW SDK 500MB compress limit exceeded.

E' un problema di sequenza di parametri?

Grazie per la pazienza!!!

 

[tb14]

Grazie del suggerimento tb14, ma non ottengo il risultato sperato:

C:\Programmi\FWTools1.0.0a8>gdal_translate -of ecw -co="LARGE_OK=TRUE" S_Antioco.tif S_Antioco.ecw
Input file size is 10975, 18341
0.ERROR 1: ECW SDK 500MB compress limit exceeded.

E' un problema di sequenza di parametri?

Grazie per la pazienza!!!

No il problema è che andavo a memoria e ... stò invecchiando

Il parametro corretto è :
-co="LARGE_OK=YES"

e non true

 

[Excalib65]

Non funge nemmeno così, sto provando a dare un'occhiata alla documentazione di gdal...

Chi prima trova la soluzione... non vince nulla ma almeno riesce a creare degli ecw di dimensione accetabile...

Grazie

 

[tb14]

Non funge nemmeno così, sto provando a dare un'occhiata alla documentazione di gdal...

Chi prima trova la soluzione... non vince nulla ma almeno riesce a creare degli ecw di dimensione accetabile...

Grazie

Ti assicuro che ho compresso file di 40000x40000 ( 4-5 Gb di file non compressi )
Stai utilizzando l'ultima versione ??

 

[Excalib65]

La versione installata è la 1.0.0a8, ma la cosa strana è che sul sito ufficiale l'ultima versione rilasciata risulta la 1.0.0a7
Ho già da tempo installato sul mio pc l' ECW Compressor 2.6, può avere a che fare col mio problema?

Booohhh???

 

[Excalib65]

La sintassi non era corretta: tu riportavi: -co ="LARGE_OK=YES", mentre la sintassi corretta è: -co "LARGE_OK=YES"

Grazie comunque tb14, sei stato grande.

 

[letdarri]

Effettivamente la discussione è molto interessante!!!
Seguendo i Vs suggerimenti sono riuscito a creare un file di unione .tiff di oltre 800 Mb, ma quando eseguo la conversione mi dà il seguente errore:

Input file size is 30000, 9500
0.ERROR 1: ECW SDK 500MB compress limit exceeded

sicuramente sbaglio qualcosa ?!?!

 

[letdarri]

Come non detto, come al solito ho letto superficialmente le risposte date da tb14!!!!!

Ho saltato un passaggio, ma perchè me ne accorgo sempre dopo azzz..!!!
Che figura solita che faccio!!!!!!

 

[geppo960]

Un grosso applauso a TB14 per questo topic e' d'obbligo.
una domandina pero'. Dove sbaglio in questa richiesta?

C:\Programmi\FWTools1.0.0a8>gdal_translate -of GTiff -co "TILED=YES" -co "TFW=YE
S" -projwin 733000 4724000 735000 4726000 ecwp://ww3.pcn.minambiente.it/Cartogra
fia/Topografia/IGM_25000/wgs84_33/12.ecw gradoli.tif
Input file size is 88000, 96000
Computed -srcwin 213200 10400 800 -799 from projected window.
Computed -srcwin falls outside raster size of 88000x96000.
C:\Programmi\FWTools1.0.0a8>

Mi torna sempre l'errore, ho anche diminuito l'area di richiesta.
grazie

 

[sembola]

Un grosso applauso a TB14 per questo topic e' d'obbligo.
una domandina pero'. Dove sbaglio in questa richiesta?

C:\Programmi\FWTools1.0.0a8>gdal_translate -of GTiff -co "TILED=YES" -co "TFW=YE
S" -projwin 733000 4724000 735000 4726000 ecwp://ww3.pcn.minambiente.it/Cartogra
fia/Topografia/IGM_25000/wgs84_33/12.ecw gradoli.tif
Input file size is 88000, 96000
Computed -srcwin 213200 10400 800 -799 from projected window.
Computed -srcwin falls outside raster size of 88000x96000.
C:\Programmi\FWTools1.0.0a8>

Mi torna sempre l'errore, ho anche diminuito l'area di richiesta.
grazie

Direi che la seconda coordinata nord (4726000) è superiore alla prima (4724000) mentre dev'essere inferiore. La prima coppie di coordinate è l'angolo nordovest, il secondo l'angolo sudest. Siccome la distanza dall' equatore aumenta con la latitudine, la prima coordinata nord dev' essere per forza superiore alla seconda. Infatti se guardi bene le dimensioni della "finestra" estratta sono 800 di longitudine e -799 di latitudine.

 

[geppo960]

Direi che la seconda coordinata nord (4726000) è superiore alla prima (4724000) mentre dev'essere inferiore. La prima coppie di coordinate è l'angolo nordovest, il secondo l'angolo sudest. Siccome la distanza dall' equatore aumenta con la latitudine, la prima coordinata nord dev' essere per forza superiore alla seconda. Infatti se guardi bene le dimensioni della "finestra" estratta sono 800 di longitudine e -799 di latitudine.

Ma difatti la cosa sembrava pure a me ma ci deve essere un problema di lettura da compe, perche infatti il 2^ punto che richiedo è piu' a sud-est del primo ma mi dice proprio 4726000. Bho strano.
Cmq grazie.

 

[douphen]

Hai ragione l'avevo dimenticato .L'ho appena fatto per il 25000. il link per l'umbria è:

ecwp://ww3.pcn.minambiente.it/Cartografia/Topografia/IGM_25000/wgs84_33/10.ecw

Per le altre regioni valgono le stesse regole delle ortofoto. Per le altre scale ... non ho provato ... ma credo che basti un pò di immaginazione

PS. La funzione di decodifica non si chiama decode ma "num_to_str" e per ottenere l'indirizzo basta modificarla come :

function num_to_str(strcri){

str_out = "";
num_out = strcri;

for(i = 0; i < num_out.length; i += 2) {
num_in = parseInt(num_out.substr(i,[2])) + 23;
num_in = unescape('%' + num_in.toString(16));
str_out += num_in;
}
alert( str_out); strdec = unescape(str_out);
return strdec;
}

Perchè non riesco a importare il 100000, il 250000 e le altre per il Trentino anche mettendo il 04. Sono riuscito solo con le ortofoto e il 25000.
Grazie

 

[douphen]

Fatto, spero di essere stato utile:

Foto B/N
ecwp://ww3.pcn.minambiente.it/Cartografia/Fotogrammetria/Regioni/BN/wgs84_32/03.ecw

IGM 1:25.000
ecwp://ww3.pcn.minambiente.it/Cartografia/Topografia/IGM_25000/wgs84_32/03.ecw

IGM 1:100.000
ecwp://ww3.pcn.minambiente.it/Cartografia/Topografia/IGM_100000/wgs84_32/03.ecw

IGM 1:250.000
ecwp://ww3.pcn.minambiente.it/Cartografia/Topografia/IGM_250000/wgs84_32/03.ecw

IGM 1:500.000
ecwp://ww3.pcn.minambiente.it/Cartografia/Topografia/IGM_500000/wgs84_32/03.ecw

Carta Politico/Fisica
ecwp://ww3.pcn.minambiente.it/Cartografia/Topografia/Carta_fisico_Politica/wgs84_32/03.ecw

Atlante De Agostini
ecwp://ww3.pcn.minambiente.it/Cartografia/Topografia/dago/wgs84_32/03.ecw


Ovviamente bisogna cambiare la regione e se necessario il fuso.
Ciao!!!!

Scusa TB14 dovevo chiederlo a lui perchè non funziona per il 100000, 250000, 500000, politica e dago anche se sostituisco il numero 04 per il Trentino.

Ciao

Non sono stato molto attento a quello che avete detto sull'importazione dei file perchè non ci capisco niente!!!!!!!
Qualcuno sa passarmi il file da mettere sul mio notebok del Trentino perchè io non sono in grado di fare tutto quello che avete detto.

Hekp me please!!!!!!!!!!!!!!!!!!!!1

 

[KORKI]

Perchè non riesco a importare il 100000, il 250000 e le altre per il Trentino anche mettendo il 04. Sono riuscito solo con le ortofoto e il 25000.
Grazie

Fatto, spero di essere stato utile:

Foto B/N
ecwp://ww3.pcn.minambiente.it/Cartografia/Fotogrammetria/Regioni/BN/wgs84_32/03.ecw

IGM 1:25.000
ecwp://ww3.pcn.minambiente.it/Cartografia/Topografia/IGM_25000/wgs84_32/03.ecw

IGM 1:100.000
ecwp://ww3.pcn.minambiente.it/Cartografia/Topografia/IGM_100000/wgs84_32/03.ecw

IGM 1:250.000
ecwp://ww3.pcn.minambiente.it/Cartografia/Topografia/IGM_250000/wgs84_32/03.ecw

IGM 1:500.000
ecwp://ww3.pcn.minambiente.it/Cartografia/Topografia/IGM_500000/wgs84_32/03.ecw

Carta Politico/Fisica
ecwp://ww3.pcn.minambiente.it/Cartografia/Topografia/Carta_fisico_Politica/wgs84_32/03.ecw

Atlante De Agostini
ecwp://ww3.pcn.minambiente.it/Cartografia/Topografia/dago/wgs84_32/03.ecw


Ovviamente bisogna cambiare la regione e se necessario il fuso.
Ciao!!!!

Attento al fuso come viene detto alla fine. Se la zona scelta ricade nel fuso 33, sostituisci wgs84_32 con wgs84_33.

 

[HDOP]

Attento al fuso come viene detto alla fine. Se la zona scelta ricade nel fuso 33, sostituisci wgs84_32 con wgs84_33.

Il trentino cade nel fuso 32.
Non è un problema di fusi.
Neanche io riesco ad accedere ai file: Atlante De Agostini, Carta Politico/Fisica, ecc... Per qualsiasi regione.
Riesco solo a importare le ortofoto e il 25000 IGM.
Per la Calabria e la Sicilia non riesco a importare niente!!
Qualcuno è riuscito a importare le ortofoto della Calabria?

 

[KORKI]

Attento al fuso come viene detto alla fine. Se la zona scelta ricade nel fuso 33, sostituisci wgs84_32 con wgs84_33.

Il trentino cade nel fuso 32.
Non è un problema di fusi.
Neanche io riesco ad accedere ai file: Atlante De Agostini, Carta Politico/Fisica, ecc... Per qualsiasi regione.
Riesco solo a importare le ortofoto e il 25000 IGM.
Per la Calabria e la Sicilia non riesco a importare niente!!
Qualcuno è riuscito a importare le ortofoto della Calabria?

Il trentino ricade sotto 2 fusi. Ad esempio le 3 Cime di Lavaredo ricadono nel 33. Questa è la zona di passaggio.



Però l'indirizzo ricade sotto wgs84_32.

 

[HDOP]

Il trentino ricade sotto 2 fusi. Ad esempio le 3 Cime di Lavaredo ricadono nel 33. Questa è la zona di passaggio.

Vedi l'allegato 14498

Però l'indirizzo ricade sotto wgs84_32.

Hai perfettamente ragione Il Trentino, il Veneto, l'Emilia Romagna e tante altre regioni ricadono su due fusi. Ma l'ortofoto del Trentino è proiettatta sul fuso 32.
La tabella seguente riporta i fusi rispetto ai quali sono state proiettate le diverse regioni Italiane:
1 Abruzzo 33
2 Basilicata 33
3 Calabria 33
4 Campania 33
5 Emilia Romagna 32
6 Friuli Venezia Giulia 33
7 Lazio 33
8 Liguria 32
9 Lombardia 32
10 Marche 33
11 Molise 33
12 Piemonte 32
13 Puglia 33
14 Sardegna 32
15 Sicilia 33
16 Toscana 32
17 Trentino Alto Adige 32
18 Umbria 33
19 Valle D&#8217;Aosta 32
20 Veneto 32

A parte i fusi, sei riuscito ad accedere ai file: Atlante De Agostini, Carta Politico/Fisica, ecc... ?