date:20200902

Re: [mapserver-users] Use Ground Control Points from GeoTIFF

2020-09-02 Thread 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: 
EPSG:3006 EPSG:900913"

             "wms_extent"           "1225000 613 191 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 >:


What does the mapfile look like?

On Tue, Sep 1, 2020 at 10:26 AM Jesper Kempe mailto:jesper.ke...@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@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/mapserver-users

Re: [mapserver-users] Use Ground Control Points from GeoTIFF

2020-09-02 Thread Rahkonen Jukka (MML)

Try to warp into .VRT, that takes so little time and disk space that it could 
be a workaround if it works.

-Jukka-

Lähettäjä: Jesper Kempe 
Lähetetty: keskiviikko 2. syyskuuta 2020 11.40
Vastaanottaja: Rahkonen Jukka (MML) 
Kopio: Steve Lime ; Mapserver 

Aihe: Re: [mapserver-users] Use Ground Control Points from GeoTIFF

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) 
mailto:jukka.rahko...@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 
mailto:mapserver-users-boun...@lists.osgeo.org>>
 Puolesta Jesper Kempe
Lähetetty: keskiviikko 2. syyskuuta 2020 10.52
Vastaanottaja: Steve Lime mailto:sdl...@gmail.com>>
Kopio: Mapserver 
mailto:mapserver-users@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: EPSG:3006 
EPSG:900913"
"wms_extent"   "1225000 613 191 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 
mailto:sdl...@gmail.com>>:
What does the mapfile look like?

On Tue, Sep 1, 2020 at 10:26 AM Jesper Kempe 
mailto:jesper.ke...@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@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/mapserver-users
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/mapserver-users

Re: [mapserver-users] Use Ground Control Points from GeoTIFF

2020-09-02 Thread 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.rahko...@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  
> *Puolesta
> *Jesper Kempe
> *Lähetetty:* keskiviikko 2. syyskuuta 2020 10.52
> *Vastaanottaja:* Steve Lime 
> *Kopio:* Mapserver 
> *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:
> EPSG:3006 EPSG:900913"
> "wms_extent"   "1225000 613 191 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 :
>
> What does the mapfile look like?
>
>
>
> On Tue, Sep 1, 2020 at 10:26 AM Jesper Kempe 
> 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@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/mapserver-users
>
>
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/mapserver-users

Re: [mapserver-users] Use Ground Control Points from GeoTIFF

2020-09-02 Thread Rahkonen Jukka (MML)

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  Puolesta 
Jesper Kempe
Lähetetty: keskiviikko 2. syyskuuta 2020 10.52
Vastaanottaja: Steve Lime 
Kopio: Mapserver 
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: EPSG:3006 
EPSG:900913"
"wms_extent"   "1225000 613 191 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 
mailto:sdl...@gmail.com>>:
What does the mapfile look like?

On Tue, Sep 1, 2020 at 10:26 AM Jesper Kempe 
mailto:jesper.ke...@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@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/mapserver-users
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/mapserver-users

Re: [mapserver-users] Use Ground Control Points from GeoTIFF

2020-09-02 Thread Jesper Kempe

Hi

This is the mapfile:

MAP
NAME "map"
WEB
METADATA
"wms_title""WMS mapserver"
"wms_srs"  "EPSG:4326 EPSG:3021 EPSG: EPSG:3006
EPSG:900913"
"wms_extent"   "1225000 613 191 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 :

> What does the mapfile look like?
>
> On Tue, Sep 1, 2020 at 10:26 AM Jesper Kempe 
> 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@lists.osgeo.org
>> https://lists.osgeo.org/mailman/listinfo/mapserver-users
>
>
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/mapserver-users

Re: [mapserver-users] Use Ground Control Points from GeoTIFF

Re: [mapserver-users] Use Ground Control Points from GeoTIFF

Re: [mapserver-users] Use Ground Control Points from GeoTIFF

Re: [mapserver-users] Use Ground Control Points from GeoTIFF

Re: [mapserver-users] Use Ground Control Points from GeoTIFF

5 matches

Site Navigation

Mail list logo

Footer information