[MapServer-users] Problems with vrt created from GCP raster

Scott public at postholer.com
Tue Dec 6 10:19:19 PST 2022


The only thing I can think of is create an entirely new .tif from your 
gcp .tif, then it shouldn't be a problem.

However, this could result a very large set of new sentinel tifs. Bad.

On 12/6/22 10:07, Marcin Niemyjski wrote:
> Scott,
> 
>  1. gdalbuildvrt doesnt work with raster georeferenced by ground control
>     points. Thats why I am using gdal warp ->
>     https://gis.stackexchange.com/questions/283595/transform-rectify-unprojected-raster-using-target-extent-co-ordinates <https://gis.stackexchange.com/questions/283595/transform-rectify-unprojected-raster-using-target-extent-co-ordinates>
>  2. Source files are already in cloud optimized format
> 
> Do you have any other ideas?
> 
> Best,
> Marcin
> 
> Wysyłane z aplikacji Outlook dla systemu iOS <https://aka.ms/o0ukef>
> ------------------------------------------------------------------------
> *Od:* MapServer-users <mapserver-users-bounces at lists.osgeo.org> w 
> imieniu użytkownika Scott <public at postholer.com>
> *Wysłane:* Tuesday, December 6, 2022 6:39:00 PM
> *Do:* mapserver-users at lists.osgeo.org <mapserver-users at lists.osgeo.org>
> *Temat:* Re: [MapServer-users] Problems with vrt created from GCP raster
> My "first is worst" suggestion would be to create a .vrt for each day:
> 
> gdalbuildvrt Day1.vrt [date]*.tif
> gdalbuildvrt Day2.vrt [date]*.tif
> ...
> 
> Next, in your .map file use runtime substitution:
> 
> DATA "/file/path/Day%day%.vrt"
> 
> If you really want to make this fast, convert all your source .tif's to
> format COG. All you tif's will be tiled:
> 
> gdalwarp -f COG -co COMPRESS=DEFLATE -co BIGTIFF=YES sent.tif sent_cog.tif
> 
> 
> 
> On 12/6/22 09:18, Marcin Niemyjski via MapServer-users wrote:
>> Hello,
>> 
>> up to now, I have been working with VRT as a carrier for the rasters I 
>> wanted to use in the mapfile. Now in order to use the rasters that were 
>> provided to me (they are GCP georeferenced rasters) I have to convert 
>> them to VRT(via gdal warp), hence my question - how can I mosaic the 
>> rasters within one layer in Mapserver (one layer = one day of photos 
>> from Sentinel-2 for example). I can't create vrt to existing vrt, and 
>> creating a tileindex for each layer generates too many files. I wanted 
>> to create one vrt gathering together all of raster but gdal_warp uses 
>> only the first one as data source.
>> 
>> What can I do?
>> 
>> Best regards,
>> Marcin
>> 
>> 
>> 
>> _______________________________________________
>> MapServer-users mailing list
>> MapServer-users at lists.osgeo.org
>> https://lists.osgeo.org/mailman/listinfo/mapserver-users 
> <https://lists.osgeo.org/mailman/listinfo/mapserver-users>
> _______________________________________________
> MapServer-users mailing list
> MapServer-users at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/mapserver-users 
> <https://lists.osgeo.org/mailman/listinfo/mapserver-users>



More information about the MapServer-users mailing list