[mapserver-users] Use Ground Control Points from GeoTIFF

[Jeff McKenna]

I'm curious what happens when you try to generate a map image from the 
commandline: shp2img -m test.map -o ttt.png -all_debug 5

-jeff



-- 
Jeff McKenna
MapServer Consulting and Training Services
co-founder of FOSS4G
http://gatewaygeo.com/



On 2020-09-02 4:51 a.m., Jesper Kempe wrote:
> 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 
> <mailto: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
>     <mailto: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
> 
>