MML:n JPEG2000 ortoilmakuvat ja niiden käyttö

JRA · 2012-05-15 17:39:32

JPEG2000 on upea tiedostomuoto, mutta valitettavasti niin mutkikas, että maailmasta ei löydy monta hyvää ja nopeaa ohjelmakirjastoa niiden käsittelemiseen. Yksikään todella nopea kirjasto ei ole vapaa. Kuitenkin netistä löytyy ilmaisia ja laillisia ohjelmia tehokkaaseen JPEG2000-käsittelyyn vähintään omaan henkilökohtaiseen käyttöön. Tarkistakaa lisenssit ennen kuin alatte käydä kauppaa kuvanmuunnoksilla.

Kirjoitan aikanaan lisää JPEG2000:sta, mutta koska niitä nyt on jo MML:ltä saatavissa, niin annetaan tässä lyhyt ensiapu.

Nämä kannattaa ladata ja pitää tallessa ikuisesti, sillä jonain päivänä ne saattavat kadota netistä.

http://home.gdal.org/fwtools/
FWTools 2.0.4 (Kakadu ja ECW-JPEG2000)
FWTools 2.4.7 (ECW-JPEG2000 ja MrSID-JPEG2000)

Kakadu-demo-ohjelmat
http://www.kakadusoftware.com/index.php … &Itemid=22

Tuoreista GDAL-versioissa on oletuksena vain avoin OpenJPEG-JPEG2000 ajuri, eikä sillä saa MML:n tiedostoja edes puretuksi tiffeiksi. Ja koska mm. Quantum GIS käyttää GDAL:ia, niin vakio-JPEG2000 ajurilla tulee saamaan kurjan kokemuksen.

Tässä pienet nopeusvertailut MML:n mustavalkokuvan ja värikuvan purkamisesta eri GDAL-ajureilla geotiffiksi:

GDAL-JP2KAK (Kakadu) (FWtools 2.0.4) mustavalko 100 sekuntia, värikuva 340 sekuntia
GDAL-JP2ECW (ECW-JPEG2000) FWTools 2.4.7 mustavalko 160 sekuntia, värikuva 514 sekuntia
GDAL-JP2MrSID FWTools 2.4.8 mustavalko 130 sekuntia, värikuva 263 sekuntia

MrSID-ajurin käyttämiseksi täytyy antaa FWTools 2.4.7:lle ympäristömuuttuja GDAL_SKIP=JP2ECW
Winsdowsissa siis komentoriviltä:
set GDAL_SKIP=JP2ECW

Vertailuaika Kakadun kdu_expand-ohjelmalla mustavalko 71 sekuntia, värikuva 194 sekuntia

Yksi toisto, mutta tulokset vastaavat omia kokemuksia. GDAL-Kakadun värikuva-ajan hitaus on ainoa, mikä kummastuttaa.

Ohjelma, joka kannattaa myös pitää tallessa, jos on sellaisen joskus ER Mapperilta ladannut, on nimeltään
ECW JPEG 2000 Compressor 7.0
Se oli aikanaan vapaasti ladattavissa ja sen käyttö oli ja on täysin vapaata, mutta sitä ei ole koskaan saanut jakaa eteenpäin.

Sitten vielä luultavasti nämä OSSIM-versiot tulevat Kakadu-laajennoksen kanssa, mutta en pysty nyt tarkistamaan:
http://download.osgeo.org/ossim/install … taller.exe (Win 32-bit)
http://download.osgeo.org/ossim/install … 8_12-4.exe (Win 64-bit)

OSSIM:in kehittäjillä on täysi Kakadu-lisenssi, ja siksi OSSIM:in Kakadu-ominaisuuksien käytöllä ei ole mitään rajaa, eli OSSIM:lla saa tehdä JPEG2000-bisnestä aivan vapaasti. Tosin sen Kakadu-ohjain ei ole vielä ihan huippuluokkaa.

OpenJUMP Plus-versio avaa hienosti MML:n JPEG2000 kuvat. Juuri nyt hyvä versio ladattavaksi on 1.5.2 RC3, joka löytyy omasta kansiostaan täältä
http://sourceforge.net/projects/jump-pi … snapshots/

Aiheeseen liittyviä vanhoja tekstiluonnoksia
http://latuviitta.org/documents/OpenJUM … EG2000.pdf
http://latuviitta.org/documents/OpenJUM … r_tool.pdf

Kosmo GIS ja gvSIG pärjäävät oletettavasti yhtä hyvin, koska ne käyttävät samoja JPEG2000-komponentteja.

Ja ettei tarvitse ruveta väittämään, että JPEG2000 on hidas tiedostomuoto, niin kannattaa katsella kuvia Kakadu-demopakettiin kuuluvalla kdu_show-ohjelmalla.

Zer · 2012-05-30 06:46:38

Windowsin QGis osaa näyttää MML:n JPEG2000-ilmakuvat oikein ja oikeassa paikassa, ainakin OSGEO4W-asennuksena. Tämä GDAL-versiolla 1.9.1.

En kylläkään muista, mitä kaikkea lisäpaketteja asensin (valinnan Advanced Install kautta).

Sitten on mahdollista omiin tarpeisiin muuntaa häviöllisesti pakatuksi.

* * *

Jos on tarpeen, niin jp2-tiedostoista on mahdollista irroittaa metadata ExifToolilla. (Windowsille on versio, joka ei vaadi perliä - en ole testannut.) Esimerkiksi Savonlinnan keskustan ilmakuvan olennainen data näyttää tältä:

File Name                       : N5311A.jp2
Image Height                    : 12000
Image Width                     : 12000
Pixel Scale                     : 0.5 0.5 0
Model Tie Point                 : 0 0 0 596000 6864000 0
Geo Tiff Version                : 1.1.0
GT Model Type                   : Projected
GT Raster Type                  : Pixel Is Area
GT Citation                     : ETRS89 / TM35FIN(E,N)
Geog Citation                   : ETRS89
Geog Angular Units              : Angular Degree
Projected CS Type               : Unknown (3067)
Proj Linear Units               : Linear Meter
Image Size                      : 12000x12000

JRA · 2012-05-30 14:33:53

Zer wrote:

Windowsin QGis osaa näyttää MML:n JPEG2000-ilmakuvat oikein ja oikeassa paikassa, ainakin OSGEO4W-asennuksena. Tämä GDAL-versiolla 1.9.1.
En kylläkään muista, mitä kaikkea lisäpaketteja asensin (valinnan Advanced Install kautta).

Pitää paikkansa. Näyttäisi siltä, että kuvat saa nykyviin asentamalla tämän:
gdal-mrsid: MrSID Raster Plugin for GDAL

Asennus ei vaadi mitään muuta, esimerkiksi erikseen hankittavien dll:ien laittamista oikeaan paikkaan kuten ECW Raster Plugin vaatii. Tämä oli hyvä tieto.

Gdalinfo listaa enemmänkin metadataa kaikenlaisista tukemistaan kuvista.

D:\100t_deflate>gdalinfo  UX51L_RVK_5.tif
Driver: GTiff/GeoTIFF
Files: UX51L_RVK_5.tif
Size is 9600, 9600
Coordinate System is:
PROJCS["ETRS89 / TM35FIN(E,N)",
    GEOGCS["ETRS89",
        DATUM["European_Terrestrial_Reference_System_1989",
            SPHEROID["GRS 1980",6378137,298.2572221010002,
                AUTHORITY["EPSG","7019"]],
            TOWGS84[0,0,0,0,0,0,0],
            AUTHORITY["EPSG","6258"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4258"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",27],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","3067"]]
Origin = (500000.000000000000000,7770000.000000000000000)
Pixel Size = (5.000000000000000,-5.000000000000000)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_XRESOLUTION=254
  TIFFTAG_YRESOLUTION=254
Image Structure Metadata:
  COMPRESSION=DEFLATE
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  500000.000, 7770000.000) ( 27d 0' 0.00"E, 70d 2'13.22"N)
Lower Left  (  500000.000, 7722000.000) ( 27d 0' 0.00"E, 69d36'23.66"N)
Upper Right (  548000.000, 7770000.000) ( 28d15'34.32"E, 70d 1'57.21"N)
Lower Right (  548000.000, 7722000.000) ( 28d14' 2.67"E, 69d36' 8.02"N)
Center      (  524000.000, 7746000.000) ( 27d37'24.25"E, 69d49'14.49"N)
Band 1 Block=9600x1 Type=Byte, ColorInterp=Red
Band 2 Block=9600x1 Type=Byte, ColorInterp=Green
Band 3 Block=9600x1 Type=Byte, ColorInterp=Blue

Last edited by JRA (2012-05-30 19:07:22)

VilleT · 2012-08-20 15:54:19

Kirjoittelenpa tähän kun kummiskin asiaa läheltä liipataan. Elikkäs asiani koskee semmoista kun olen tuolla JOSMilla täydennellyt paikallisia mehtäpolkuja yms. sitä mukaa kun niitä on tullut koluttua ja siinä apuna olen käyttänyt jonkun verran vielä bingin ilmakuvia. Mielenkiinto näitä ortokuvia kohtaan heräsi kun ne näyttäisi olevan hivenen tarkempia ja tuoreempia niin siinäpä sitten aloin tohkeissani OSGeo4W ohjelmalla noita .jp2 tiedostoja muuntamaan .tiffeiksi, että saisi ne JOSMissa näkymään.

Ohjelma teki muunnoksen ja samalla yksittäisen tiedoston kooksi muodostui n. 500Mt ja tämähän ei sitten enää oikein tunnu JOSMin kanssa pelittävän ainakaan omalla koneella kun ei mukamas muisti riitä kuvien latailuun vaikka -Xmx###M-komennolla sitä parhaimmillaan ~1600Mt tarjosin.

Joten kysymykseni lyhykäisyydessään olisi, että mitä minun täytyy tehdä jotta nämä hienot ortokuvat omassa windows koneessani mahdollisimman vähällä kikkailulla saisin JOSMin kanssa toimimaan ?

Kiitokset etukäteen.

edit: Hieman keskustelupalstaa tutkimallahan selvisi, että loistava kapsi tarjoaa moisen palvelun erittäin kätevästi ja ohjeet näiden ortokuvien käyttöön JOSM:lla löytyy sivuston alapuoliskolta.

http://kartat.kapsi.fi/#

Last edited by VilleT (2012-08-20 18:57:41)

JRA · 2012-08-20 19:54:18

Tästä pitäisi varmaan kirjoittaa ohje, mutta GDAL:ia käyttäville ohjelmille (QGis, Mapserver, Merkaartor ym.) saadaan pienet ja nopeasti toimivat ilmakuvat näillä kahdella komennolla:

gdal_translate -of GTiff -co PHOTOMETRIC=YCBCR -co tiled=yes
-co COMPRESS=JPEG input.jp2 output.tif

gdaladdo --config COMPRESS_OVERVIEW JPEG --config PHOTOMETRIC_OVERVIEW YCBCR
output.tif 2 4 8 16 32 64

Ensimmäinen tekee sisäisesti laatoitetun tiff:in, joka käyttää jpeg-pakkausta, toinen tekee kuvan käytöstä nopeampaa kaukaa katsottuna, koska tiff:iin lasketaan 1, 2, 4, 8, 16 ja 32 metrin pikselikoon versiot heti saataville.

OpenStreetMap Forum

#1 2012-05-15 17:39:32

MML:n JPEG2000 ortoilmakuvat ja niiden käyttö

#2 2012-05-30 06:46:38

Re: MML:n JPEG2000 ortoilmakuvat ja niiden käyttö

#3 2012-05-30 14:33:53

Re: MML:n JPEG2000 ortoilmakuvat ja niiden käyttö

#4 2012-08-20 15:54:19

Re: MML:n JPEG2000 ortoilmakuvat ja niiden käyttö

#5 2012-08-20 19:54:18

Re: MML:n JPEG2000 ortoilmakuvat ja niiden käyttö

Board footer