[mapserver-users] Use Ground Control Points from GeoTIFF

[Jesper Kempe]

Yes, a gdalwarp on the file works fine but I would really like to avoid
that.
A "gdalwarp -t_srs "EPSG:3857" sat.tif sat_3857.tif" and use the file
sat_3857.tif works.

The file I want to show in Mapserver is available at
https://drive.google.com/file/d/1t3fM0Ej1FGwrdqbvg5EhIFcpnNOfKdQS/view?usp=sharing
if anyone wants to give it a try.

//Jesper

Den ons 2 sep. 2020 kl 10:26 skrev Rahkonen Jukka (MML) <
jukka.rahkonen at maanmittauslaitos.fi>:

> Hi,
>
>
>
> I am not suggesting this as a final solution yet, but gather some
> knowledge and try if image works when a) warped into GeoTIFF with gdalwarp
> b) warped into .vrt
>
>
>
> -Jukka Rahkonen-
>
>
>
> *Lähettäjä:* mapserver-users <mapserver-users-bounces at lists.osgeo.org> *Puolesta
> *Jesper Kempe
> *Lähetetty:* keskiviikko 2. syyskuuta 2020 10.52
> *Vastaanottaja:* Steve Lime <sdlime at gmail.com>
> *Kopio:* Mapserver <mapserver-users at lists.osgeo.org>
> *Aihe:* Re: [mapserver-users] Use Ground Control Points from GeoTIFF
>
>
>
> Hi
>
>
>
> This is the mapfile:
>
>
>
> MAP
>     NAME "map"
>     WEB
>         METADATA
>             "wms_title"            "WMS mapserver"
>             "wms_srs"              "EPSG:4326 EPSG:3021 EPSG:9999
> EPSG:3006 EPSG:900913"
>             "wms_extent"           "1225000 6130000 1910000 7675000"
>             "wms_enable_request"   "*"
>             "wfs_enable_request"   "*"
>         END
>     END
>
>     PROJECTION
>         "init=epsg:4326"
>     END
>
>     EXTENT -180 -90 180 90
>
>     LAYER
>       NAME "sat"
>       DATA "sat.tif"
>       TYPE RASTER
>       STATUS ON
>
>       METADATA
>                 "wms_title" "sat"
>       END
>
>       PROJECTION
>         "init=epsg:4326"
>       END
>       EXTENT 17.0473154408269387 62.1606759790530177 27.5792082876237146
> 66.7744046276809797
>
>     END
>
> END
>
>
>
> I have also tried both without PROJECTION in LAYER and with
>
> PROJECTION
>
>      AUTO
>
> END
>
>
>
> None of them works.
>
>
>
> EXTENT in the mapfile is taken from Qgis metadata of the layer.
>
>
>
>
>
> Thanks
>
> Jesper
>
>
>
> Den tis 1 sep. 2020 kl 23:41 skrev Steve Lime <sdlime at gmail.com>:
>
> What does the mapfile look like?
>
>
>
> On Tue, Sep 1, 2020 at 10:26 AM Jesper Kempe <jesper.kempe at gmail.com>
> wrote:
>
> Hi
>
>
>
> I have a GeoTiff file that gdalinfo presents like this:
>
>
>
> gdalinfo sat.tif
> Driver: GTiff/GeoTIFF
> Files: sat.tif
> Size is 10474, 10717
> Coordinate System is `'
> GCP Projection =
> GEOGCS["WGS 84",
>     DATUM["WGS_1984",
>         SPHEROID["WGS 84",6378137,298.257223563,
>             AUTHORITY["EPSG","7030"]],
>         AUTHORITY["EPSG","6326"]],
>     PRIMEM["Greenwich",0],
>     UNIT["degree",0.0174532925199433],
>     AUTHORITY["EPSG","4326"]]
> GCP[  0]: Id=1, Info=
>           (0,0) -> (27.5783052509812,65.9184735836942,229.992792810313)
> GCP[  1]: Id=2, Info=
>           (524,0) -> (27.1416007200269,65.9731480840222,172.994255514815)
> GCP[  2]: Id=3, Info=
>           (1048,0) -> (26.7005564512604,66.0268770755367,152.994545626454)
> GCP[  3]: Id=4, Info=
>           (1572,0) -> (26.257537545879,66.0793551447287,131.994995356537)
> GCP[  4]: Id=5, Info=
>           (2096,0) -> (25.8102107755995,66.1308402121171,163.993212723173)
>
> ....
>
>
>
> with 482 GCPs.
>
>
>
>
>
> The file opens ok in Qgis and is located in the right location.
>
> But I can´t in any way get it to show correctly in Mapserver.
>
>
>
> Is this possible to do? Or do I have to gdalwarp the file to a known
> projection and set the PROJECTION of the layer?
>
>
>
> Thanks
>
> Jesper
>
>
>
>
>
> _______________________________________________
> mapserver-users mailing list
> mapserver-users at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/mapserver-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20200902/95b4a9b4/attachment.htm>