Vairāk

Rgba ģeotifu flīzēšana, izmantojot vrt, nevis gdal_warp

Rgba ģeotifu flīzēšana, izmantojot vrt, nevis gdal_warp


Esmu uzmanīgi nogriezis apkakli no USGS ģeotifu komplekta, un es vēlētos tos izlaist ar gdal2tiles, bet izveidot vrt, kas pasaka alfa slāni, vai arī tas, ka gdal2tiles pārklājas divus avota failus vienā flīzē, nedarbojas. Dažas lietas, kuras esmu izmēģinājis, padarīja visus melnos pikseļus caurspīdīgus. Citi mēģinājumi padarīja manu avota failu caurspīdīgās robežas tā, it kā tās būtu vienkrāsaini baltas, tādēļ, ja viena flīze ir izgatavota no diviem blakus esošiem tifiem, vienas caurspīdīgā apmale nosedz citas noderīgo saturu. Šis ir mans process. Varbūt kāds var pateikt, kāds reklāmguvuma solis man pietrūkst:

  • Lejupielādējiet parasto USGS ģeotifu vietnē Topoquest
  • gdal_translate -expand rgba -co COMPRESS = LZW orginal.tif transl.tif
  • Programmā Photoshop atveriet Translate.tif, nogrieziet apkakli, saglabājiet.
  • gdalcopyproj.py tulkots.tif sagatavots.tif

Dariet to visiem

  • gdalbuildvrt temp.vrt * .tif
  • gdal2tiles.py temp.vrt flīzes

P.S. Es tos visus varu vienmērīgi apvienot ar gdal_warp, bet šai 30 cilvēku grupai tas veido kopējo ģeotifu 51223 x 65163, un citos komplektos būs daudz izšķērdētas vietas, ja avota tifs atrodas pa diagonāli, bet gdal_warp izveidos taisnstūrveida izeja.


GDAL atbalsta maskas slāni, bet tas atbalsta arī krāsu "nav datu", kuru es mēdzu izmantot, jo tā samazina gala attēla izmēru par 25%. Tātad darbības, ko es izmantotu ar jūsu datiem, pēc apkakles nogriešanas ir:

gdal_translate -izvērst rgb -gcp  -co TILED = JĀ -a_srs  -a_nodata 255 255 255 oriģināls.tif tulkots.tif

Piezīme, es esmu ievietojis Zemes kontroles punkti (GSP) un projekcija, kas jums jānosaka empīriski. Ņemiet vērā arī, ka nodata vērtībai esmu piešķīris tīru baltu krāsu.

Šie attēli joprojām nav pareizi ģeoreferēti, jo tiem joprojām ir tikai GCP, nevis pareizas koordinātas. Tātad jums vajadzēs skrietgdalwarp:

gdalwarp -t_srs EPSG: 4326 -co TILED = YES -srcnodata 255 255 255 -dstnodata 255 255 255 Translate.tif warped.tif

Vēlreiz ņemiet vērā mezglu vērtību precīzu iestatījumu. Tas, iespējams, ir pārspīlēts, taču tas nesāp, un man tas ir bijisgdalwarpignorēt, manuprāt, netiešos nodata vales.

Tālāk jums jāskriengdalbuildvrttāpat kā jums tas ir, bet es parasti vēlreiz skaidri norādīju mezglu vērtības, lai būtu drošībā:

gdalbuildvrt -srcnodata 255 255 255 -vrtnodata 255 255 255 temp.vrt .tif

Rediģēt

Tas var būt tuvāk vēlamajam:

gdalbuildvrt -srcnodata 255 255 255 -hidenodata temp.vrt .tif

Visbeidzot, jūs varat piezvanītgdal2tiles.pykā jums ir.


Skatīties video: Plytelių klijavimas