Ocena wątku:
  • 1 Głosów - 2 Średnio
  • 1
  • 2
  • 3
  • 4
  • 5
VRT dla różnych kolorystycznie plików
17-06-2011, 14:57,
#1
VRT dla różnych kolorystycznie plików
Witam wszystkich,
używam Quantuma już od jakiegoś czasu i dzięki poradom z forum i osobistym borysa mam już pewne osiągnięcia w serwowaniu rastrów (ortofoto) przez WMS w LAN-ie.

Mój problem polega na tym, że z geodezji dostałem pliki dla danego obszaru (geotiff) o różnym poziomie aktualności (od 2004 do 2009 roku - B&W, i chyba RGB). Przygotowałem dla wszystkich piramidy, ale przy próbie stworzenia pliku VRT dla powiedzmy 3 rodzajów GeoTiff-ów wtyczka do tworzenia wirtualnych rastrów przyjmuje tylko jeden rodzaj z zadanych plików i dla niego tworzy wirtualny raster. pozostałe pliki są pomijane. Może ktoś zna jakiś sposób na stworzenie pojedynczego VRT lub ujednolicenie plików.

Pozdrawiam
17-06-2011, 15:47,
#2
RE: VRT dla różnych kolorystycznie plików
Witam
Spróbuj zaznaczyć opcję 'jako osobne kanały'.
Wg pomocy polecenia gdalbuildvrt sprawdza poprawność wszystkich warstw pod względem kolorów, rozdzielczości, liczby kanałów itp. Jeśli któreś warstwy odbiegają charakterystyką są pomijane. Można tego uniknąć dodając parametr -separate do lini poleceń.
17-06-2011, 15:51,
#3
RE: VRT dla różnych kolorystycznie plików
A jak nie pójdzie to przekonwertuj B&W do RGB i pójdzie.
27-06-2011, 13:15,
#4
RE: VRT dla różnych kolorystycznie plików
Z parametrem - separate pokazywały się różne, bliżej dla mnie nie jasne, błędy, poszedłem więc drugą drogą Odoakra i faktycznie udało mi się przekonwertować B&W na RGB, ale ich rozmiar wyjściowy, biorąc pod uwagę rozdzielczość jest powalający - dla pojedynczego pliku skok z 3 MB przed konwersją i piramidowaniem do 23 MB po wszystkich zabiegach. Całość obszaru mapy rozrósł się do 17 GB, ale nie to jest istotne. Plik VRT się tworzy i działa tak jak powinien.

Kolejnym problemem jest to, że mam 3 pliki (z ponad 300) , które mają zapisaną projekcję w inny sposób niż pozostałe - w obu przypadkach jest to WGS84, ale w zapisany inaczej w parametrach pliku. No i znowu ten sam problem - gdalbuildvrt pomija właśnie te 3. Przy próbie dostosowania układu wsp. dla tych 3, nadal są różnice i pliki są pomijane w VRT. Natomiast gdy użyję gdalwarp dla wszystkich kostek wyrzuca mi obszar mapy poza wszelkie układy i nakłada jakąś bardzo dziwną georeferencję pomimo, iż gdalinfo twierdzi, że jest to WGS84.

W tej chwili próbuję konwersji całości na układ 92
27-06-2011, 14:08,
#5
RE: VRT dla różnych kolorystycznie plików
Też się z tym kiedyś spotkałem.
W reprojekcji rastra gdalem możesz wymusić wejściowy układ współrzędnych: "s_srs"
Spróbuj te trzy pliki zapisać w np 1992 (gdalwarp s_srs EPSG:4326 t_srs EPSG: 2180 wejscie wyjscie)
a potem z powrotem do 4326 - wtedy definicja powinna zostać naprawiona...
27-06-2011, 14:23,
#6
RE: VRT dla różnych kolorystycznie plików
Dzięki za wskazówkę - przemieliłem wszystko do 92 i puściłem przez mapserver - jest ok. Jak będę miał więcej czasu popracuję nad optymalizacją (i Twoim sposobem na naprawę definicji) - najbardziej chciałbym w jakiś sposób ograniczyć wielkość tych plików B&W, bo takie ilości danych zabiją mi sieć, a z czasem pewnie ruszę z kafelkami.
27-06-2011, 15:15,
#7
RE: VRT dla różnych kolorystycznie plików
sieć Ci nie padnie, bo przez sieć nie idą surowe pliki. Co innego serwer i przestrzeń dyskowaSmile
Sam nie wiem jak z pliku 24bit (RGB) przejść do 8 bit w taki sposób, żeby przeszło to przez gdalvrt (czyli zgadzały się palety wszystkich plików...) może połącz wszystko w jedno i przekonwertuj do 8bit? Albo użyj ECW zamiast geotiff. Ale z tym trzeba powalczyć i chyba licencja pozwala używać go tylko na desktopach.
01-07-2011, 19:35,
#8
RE: VRT dla różnych kolorystycznie plików
(27-06-2011, 13:15)vansen napisał(a): dla pojedynczego pliku skok z 3 MB przed konwersją i piramidowaniem do 23 MB po wszystkich zabiegach.

A jaką kompresję stosujesz? Ustawienia kompresji mogą bardzo duuużo zmienić.
04-07-2011, 07:57, (Ten post był ostatnio modyfikowany: 04-07-2011, 08:18 {2} przez vansen.)
#9
RE: VRT dla różnych kolorystycznie plików
Moje ustawienia są takie:

Generowanie piramid:
Resampling - średnia
Poziomy 8 - 2 4 8 16 32 64 128 256
Bez podglądu w formacie TIFF z kompresją JPG
i bez formatu imagine (nawet nie wiem co to jest)

A przy konwersji PCT - RGB zawze idzie jeden kanał (analogowe B&W)

Masz jakieś sugestie?

Edycja:
Poniżej wrzucam wynik gdalinfo dla 2 plików które chciałbym wyświetlić na jednym VRT. Stosowałem różne metody i nie rozumiem powodu dlaczego tylko ten jeden rodzaj nie chce się poprawnie wyświetlić (3 pozostałe po różnych proponowanych na forum obróbkach wyświetlają się razem bez żadnego problemu).

Nie wyświetla się pierwszy z rodzajów plików. (na czerwono)

Driver: GTiff/GeoTIFF
Size is 9082, 9617
Coordinate System is:
PROJCS["unnamed",
GEOGCS["ETRS89",
DATUM["European_Terrestrial_Reference_System_1989",
SPHEROID["GRS 1980",6378137,298.2572221010002,
AUTHORITY["EPSG","7019"]],
AUTHORITY["EPSG","6258"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4258"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",19],
PARAMETER["scale_factor",0.9993],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",-5300000],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (706737.514506234790000,389459.296914019620000)
Pixel Size = (0.250000000000000,-0.250000000000000)
Metadata:
AREA_OR_POINT=Area
TIFFTAG_XRESOLUTION=0
TIFFTAG_YRESOLUTION=0
TIFFTAG_RESOLUTIONUNIT=1 (unitless)
Image Structure Metadata:
COMPRESSION=JPEG
Corner Coordinates:
Upper Left ( 706737.515, 389459.297) ( 21d58'7.69"E, 51d20'2.86"N)
Lower Left ( 706737.515, 387055.047) ( 21d58'2.66"E, 51d18'45.12"N)
Upper Right ( 709008.015, 389459.297) ( 22d 0'4.89"E, 51d19'59.87"N)
Lower Right ( 709008.015, 387055.047) ( 21d59'59.81"E, 51d18'42.13"N)
Center ( 707872.765, 388257.172) ( 21d59'3.76"E, 51d19'22.50"N)
Band 1 Block=256x256 Type=Byte, ColorInterp=Red
Overviews: 4541x4808, 2270x2404, 1135x1202, 567x601, 283x300, 141x150, 70x75, 35x37
Band 2 Block=256x256 Type=Byte, ColorInterp=Green
Overviews: 4541x4808, 2270x2404, 1135x1202, 567x601, 283x300, 141x150, 70x75, 35x37
Band 3 Block=256x256 Type=Byte, ColorInterp=Blue
Overviews: 4541x4808, 2270x2404, 1135x1202, 567x601, 283x300, 141x150, 70x75, 35x37
Band 4 Block=256x256 Type=Byte, ColorInterp=Alpha
Overviews: 4541x4808, 2270x2404, 1135x1202, 567x601, 283x300, 141x150, 70x75, 35x37


Driver: GTiff/GeoTIFF
Size is 9086, 9621
Coordinate System is:
PROJCS["unnamed",
GEOGCS["unnamed",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.2572235630016,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT[,0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",19],
PARAMETER["scale_factor",0.9993],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",-5300000],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (708913.250000000000000,389549.000000000000000)
Pixel Size = (0.250000000000000,-0.250000000000000)
Metadata:
AREA_OR_POINT=Area
TIFFTAG_XRESOLUTION=1
TIFFTAG_YRESOLUTION=1
TIFFTAG_RESOLUTIONUNIT=1 (unitless)
Image Structure Metadata:
COMPRESSION=JPEG
Corner Coordinates:
Upper Left ( 708913.250, 389549.000) ( 22d 0'0.19"E, 51d20'2.90"N)
Lower Left ( 708913.250, 387143.750) ( 21d59'55.11"E, 51d18'45.12"N)
Upper Right ( 711184.750, 389549.000) ( 22d 1'57.44"E, 51d19'59.88"N)
Lower Right ( 711184.750, 387143.750) ( 22d 1'52.31"E, 51d18'42.10"N)
Center ( 710049.000, 388346.375) ( 22d 0'56.26"E, 51d19'22.51"N)
Band 1 Block=256x256 Type=Byte, ColorInterp=Red
Overviews: 4543x4810, 2271x2405, 1135x1202, 567x601, 283x300, 141x150, 70x75, 35x37, 17x18, 8x9
Band 2 Block=256x256 Type=Byte, ColorInterp=Green
Overviews: 4543x4810, 2271x2405, 1135x1202, 567x601, 283x300, 141x150, 70x75, 35x37, 17x18, 8x9
Band 3 Block=256x256 Type=Byte, ColorInterp=Blue
Overviews: 4543x4810, 2271x2405, 1135x1202, 567x601, 283x300, 141x150, 70x75, 35x37, 17x18, 8x9
Band 4 Block=256x256 Type=Byte, ColorInterp=Alpha
Overviews: 4543x4810, 2271x2405, 1135x1202, 567x601, 283x300, 141x150, 70x75, 35x37, 17x18, 8x9


Podobne wątki
Wątek: Autor Odpowiedzi: Wyświetleń: Ostatni post
  Zlecę napisanie skryptu dla QGIS do automatycznej edycji wielu plików SHP Matteusz 3 16 082 30-06-2017, 22:24
Ostatni post: bociusek
  te same atrybuty w różnych warstwach phacops 3 17 175 13-11-2014, 17:47
Ostatni post: martel
  zapisywanie się plików projektu w typie plików innym niż dbf mak 3 19 381 28-07-2011, 11:54
Ostatni post: mak
  Błędy przy tworzeniu piramid i .vrt tomalos 1 11 078 19-02-2011, 20:55
Ostatni post: borys

Skocz do: