Re: [mapserver-users] How to dynamically change MapServer's GRID MININTERVAL?

2020-09-03 Thread Lime, Steve D (MNIT)
Hi David: That functionality was geared towards common use cases so I’d bet 
money it’s not supported but I’ll check the code to be sure. If you’re really 
just after 5 vs. 10 then using multiple layers would be an easy work around 
(e.g. grid5 vs grid10). Note that we’re always interested in 500 errors so if 
you are able to pass along how you triggered one that would be helpful. --Steve

From: mapserver-users  On Behalf Of 
David Miller
Sent: Thursday, September 3, 2020 8:08 AM
To: mapserver-users@lists.osgeo.org
Subject: [mapserver-users] How to dynamically change MapServer's GRID 
MININTERVAL?

I started looking into using something similar to this example of the GRID 
layer definition in a MapServer mapfile that will draw lat/lon line grid at 
least every 5 degrees:

LAYER

  NAME "grid"

  METADATA

"DESCRIPTION" "Grid"

  END

  TYPE LINE

  STATUS ON

  CLASS

NAME "Graticule"

COLOR 0 0 0

LABEL

  COLOR  255 0 0

  FONT "sans"

  TYPE truetype

  SIZE 8

  POSITION AUTO

  PARTIALS FALSE

  BUFFER 2

  OUTLINECOLOR 255 255 255

END

  END

  PROJECTION

"init=epsg:4326"

  END

  GRID

LABELFORMAT '%g°'

MININTERVAL 5

  END

END # Layer

However, I don't want the MININTERVAL value to be hardcoded and would like to 
be able to vary that value by entering a different value in the URL call from a 
web page using the map_object notation. But, I've been unsuccessful. I have 
tried this in the URL call to MapServer:

_layer[grid].grid[0]=MININTERVAL+10 #changing to 10 degrees instead of 5

using the second example at this link: 
https://mapserver.org/cgi/controls.html#cgi-mapfile-change-parameters.
 But, I either get a parsing error or an Internal Server error. I've tried 
turning debug level to 5 and looking at logs but they don't offer any help. 
Perhaps I'm just not using the correct map_object notation, but I'm not sure 
what else to use given the examples in the link. I also attempted to look at 
the MapServer code to no avail.

Perhaps I just can't change this parameter via the URL?

I'm using MapServer v 7.2.2 but it seems like the GRID layer has been available 
for a while now.

Thanks for any help/suggestions/thoughts!

Dave M
___
mapserver-users mailing list
mapserver-users@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/mapserver-users

[mapserver-users] How to dynamically change MapServer's GRID MININTERVAL?

2020-09-03 Thread David Miller
I started looking into using something similar to this example of the GRID
layer definition in a MapServer mapfile that will draw lat/lon line grid at
least every 5 degrees:

LAYER
  NAME "grid"
  METADATA
"DESCRIPTION" "Grid"
  END
  TYPE LINE
  STATUS ON
  CLASS
NAME "Graticule"
COLOR 0 0 0
LABEL
  COLOR  255 0 0
  FONT "sans"
  TYPE truetype
  SIZE 8
  POSITION AUTO
  PARTIALS FALSE
  BUFFER 2
  OUTLINECOLOR 255 255 255
END
  END
  PROJECTION
"init=epsg:4326"
  END
  GRID
LABELFORMAT '%g°'
MININTERVAL 5
  END
END # Layer

However, I don't want the MININTERVAL value to be hardcoded and would like
to be able to vary that value by entering a different value in the URL call
from a web page using the map_object notation. But, I've been unsuccessful.
I have tried this in the URL call to MapServer:

_layer[grid].grid[0]=MININTERVAL+10 #changing to 10 degrees instead of 5

using the second example at this link:
https://mapserver.org/cgi/controls.html#cgi-mapfile-change-parameters. But,
I either get a parsing error or an Internal Server error. I've tried
turning debug level to 5 and looking at logs but they don't offer any help.
Perhaps I'm just not using the correct map_object notation, but I'm not
sure what else to use given the examples in the link. I also attempted to
look at the MapServer code to no avail.

Perhaps I just can't change this parameter via the URL?

I'm using MapServer v 7.2.2 but it seems like the GRID layer has been
available for a while now.

Thanks for any help/suggestions/thoughts!

Dave M
___
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-03 Thread Rahkonen Jukka (MML)
Hi Jesper,

Poor performance may mean that your original tiff is not optimized for server 
usage and you should rewrite it in any case or at least add some overviews. 
Gdalinfo report about the original tiff and .vrt would probably give some 
useful information. If you select to warp images physically (would be my choise 
BTW) please try to do it right. You should find good GDAL commands from these 
slides 
https://www.slideshare.net/geosolutions/geoserver-on-steroids-foss4g-2015.

-Jukka Rahkonen-

Lähettäjä: Jesper Kempe mailto:jesper.ke...@gmail.com>>
Lähetetty: torstai 3. syyskuuta 2020 15.18
Vastaanottaja: Rahkonen Jukka (MML) 
mailto:jukka.rahko...@maanmittauslaitos.fi>>
Kopio: Jeff McKenna 
mailto:jmcke...@gatewaygeomatics.com>>; 
Mapserver 
mailto:mapserver-users@lists.osgeo.org>>
Aihe: Re: [mapserver-users] Use Ground Control Points from GeoTIFF

Hi Jukka

Sorry about that, I forgot that yesterday.

I did a test now and that actually works!

Changing DATA to the vrt-file and setting

PROJECTION
  "init=epsg:3857"
END

creates an image with sat.tif placed correctly!

$ shp2img -m icesat.map -o /data/temp/Jesper/ttt.png -all_debug 5
msLoadMap(): 0.014s
msDrawMap(): rendering using outputformat named png (AGG/PNG).
msDrawRasterLayerLow(sat): entering.
msResampleGDALToMap in effect: cellsize = 4202.211585
msDrawGDAL(sat): using RAW_WINDOW=0 0 12531 12770, dst=0,0,279,284
msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,2,3,4
msDrawMap(): Layer 0 (sat), 557.242s
msDrawMap(): Drawing Label Cache, 0.000s
msDrawMap() total time: 557.243s
msSaveImage(/data/temp/Jesper/ttt.png) total time: 0.119s
msFreeMap(): freeing map at 0x7b1f50.
freeLayer(): freeing layer at 0x7ba080.
freeLayer(): freeing layer at 0x7bb760.
shp2img total time: 557.378s

BUT, note 557.378s on the last row...

Also a WMS-request works but takes about 4 minutes to complete!

It looks like I need to reproject all files new GeoTIFFs to get a decent speed.

Thanks for all the replies!

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-03 Thread Jesper Kempe
Hi Jukka

Sorry about that, I forgot that yesterday.

I did a test now and that actually works!

Changing DATA to the vrt-file and setting

PROJECTION
  "init=epsg:3857"
END

creates an image with sat.tif placed correctly!

$ shp2img -m icesat.map -o /data/temp/Jesper/ttt.png -all_debug 5
msLoadMap(): 0.014s
msDrawMap(): rendering using outputformat named png (AGG/PNG).
msDrawRasterLayerLow(sat): entering.
msResampleGDALToMap in effect: cellsize = 4202.211585
msDrawGDAL(sat): using RAW_WINDOW=0 0 12531 12770, dst=0,0,279,284
msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,2,3,4
msDrawMap(): Layer 0 (sat), 557.242s
msDrawMap(): Drawing Label Cache, 0.000s
msDrawMap() total time: 557.243s
msSaveImage(/data/temp/Jesper/ttt.png) total time: 0.119s
msFreeMap(): freeing map at 0x7b1f50.
freeLayer(): freeing layer at 0x7ba080.
freeLayer(): freeing layer at 0x7bb760.
shp2img total time: 557.378s

BUT, note 557.378s on the last row...

Also a WMS-request works but takes about 4 minutes to complete!

It looks like I need to reproject all files new GeoTIFFs to get a decent
speed.

Thanks for all the replies!

Jesper


Den tors 3 sep. 2020 kl 13:31 skrev Rahkonen Jukka (MML) <
jukka.rahko...@maanmittauslaitos.fi>:

> Hi,
>
>
>
> I would appreciate if you could have a try with
>
> "gdalwarp -of vrt -t_srs "EPSG:3857" sat.tif sat_3857.vrt"
>
>
>
> -Jukka Rahkonen-
>
>
>
> *Lähettäjä:* mapserver-users  
> *Puolesta
> *Jesper Kempe
> *Lähetetty:* torstai 3. syyskuuta 2020 13.57
> *Vastaanottaja:* Jeff McKenna 
> *Kopio:* Mapserver 
> *Aihe:* Re: [mapserver-users] Use Ground Control Points from GeoTIFF
>
>
>
> shp2img creates an empty blank image with a mapfile like this:
>
>
>
> 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 15 60 30 68
> SIZE 200 200
>
> LAYER
>   NAME "sat"
>   DATA "/data/temp/Jesper/sat.tif"
>   TYPE RASTER
>   STATUS DEFAULT
>
>   METADATA
> "wms_title" "sat"
>   END
> END
>
>
>
>
>
> The output from shp2img looks like this:
>
>
>
> $ shp2img -m ttt.map -o /data/temp/Jesper/ttt.png -all_debug 5
> msLoadMap(): 0.003s
> msDrawMap(): rendering using outputformat named png (AGG/PNG).
> msDrawRasterLayerLow(sat): entering.
> msDrawRasterLayerGDAL(): Entering transform.
> msDrawRasterLayerGDAL(): src=15,10645,15,15, dst=0,0,200,200
> msDrawRasterLayerGDAL(): source raster PL (14.962,10645.462) for dst PL
> (0,0).
> msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,2,3,4
> msDrawMap(): Layer 0 (sat), 0.011s
> msDrawMap(): Drawing Label Cache, 0.000s
> msDrawMap() total time: 0.011s
> msSaveImage(/data/temp/Jesper/ttt.png) total time: 0.003s
> msFreeMap(): freeing map at 0x1d14f50.
> freeLayer(): freeing layer at 0x1d1d080.
> freeLayer(): freeing layer at 0x1d1e1b0.
> shp2img total time: 0.017s
>
>
>
>
>
> If I add:
>
>   PROJECTION
> "init=epsg:4326"
>   END
>   EXTENT 17.0473154408269387 62.1606759790530177 27.5792082876237146
> 66.7744046276809797
>
>
>
> shp2img outputs:
>
>
>
> msLoadMap(): 0.003s
> msDrawMap(): rendering using outputformat named png (AGG/PNG).
> msDrawRasterLayerLow(sat): entering.
> msDrawRasterLayerGDAL(): Entering transform.
> msDrawRasterLayerGDAL(): src=0,0,10474,10717, dst=27,63,140,61
> msDrawRasterLayerGDAL(): source raster PL (-49.554,-33.774) for dst PL
> (27,63).
> msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,2,3,4
> msDrawMap(): Layer 0 (sat), 0.025s
> msDrawMap(): Drawing Label Cache, 0.000s
> msDrawMap() total time: 0.025s
> msSaveImage(/data/temp/Jesper/ttt.png) total time: 0.004s
> msFreeMap(): freeing map at 0x2453f50.
> freeLayer(): freeing layer at 0x245c080.
> freeLayer(): freeing layer at 0x245d430.
> shp2img total time: 0.033s
>
>
>
> Then the image is shown in ttt.png but is not in the correct place and not
> rotated correctly.
>
>
>
> Any help or ideas are appreciated!
>
>
>
> Thanks
>
> Jesper
>
>
>
> Den ons 2 sep. 2020 kl 18:23 skrev Jeff McKenna <
> jmcke...@gatewaygeomatics.com>:
>
> 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"  

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

2020-09-03 Thread Rahkonen Jukka (MML)
Hi,

I would appreciate if you could have a try with
"gdalwarp -of vrt -t_srs "EPSG:3857" sat.tif sat_3857.vrt"

-Jukka Rahkonen-

Lähettäjä: mapserver-users  Puolesta 
Jesper Kempe
Lähetetty: torstai 3. syyskuuta 2020 13.57
Vastaanottaja: Jeff McKenna 
Kopio: Mapserver 
Aihe: Re: [mapserver-users] Use Ground Control Points from GeoTIFF

shp2img creates an empty blank image with a mapfile like this:

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 15 60 30 68
SIZE 200 200

LAYER
  NAME "sat"
  DATA "/data/temp/Jesper/sat.tif"
  TYPE RASTER
  STATUS DEFAULT

  METADATA
"wms_title" "sat"
  END
END


The output from shp2img looks like this:

$ shp2img -m ttt.map -o /data/temp/Jesper/ttt.png -all_debug 5
msLoadMap(): 0.003s
msDrawMap(): rendering using outputformat named png (AGG/PNG).
msDrawRasterLayerLow(sat): entering.
msDrawRasterLayerGDAL(): Entering transform.
msDrawRasterLayerGDAL(): src=15,10645,15,15, dst=0,0,200,200
msDrawRasterLayerGDAL(): source raster PL (14.962,10645.462) for dst PL (0,0).
msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,2,3,4
msDrawMap(): Layer 0 (sat), 0.011s
msDrawMap(): Drawing Label Cache, 0.000s
msDrawMap() total time: 0.011s
msSaveImage(/data/temp/Jesper/ttt.png) total time: 0.003s
msFreeMap(): freeing map at 0x1d14f50.
freeLayer(): freeing layer at 0x1d1d080.
freeLayer(): freeing layer at 0x1d1e1b0.
shp2img total time: 0.017s


If I add:
  PROJECTION
"init=epsg:4326"
  END
  EXTENT 17.0473154408269387 62.1606759790530177 27.5792082876237146 
66.7744046276809797

shp2img outputs:

msLoadMap(): 0.003s
msDrawMap(): rendering using outputformat named png (AGG/PNG).
msDrawRasterLayerLow(sat): entering.
msDrawRasterLayerGDAL(): Entering transform.
msDrawRasterLayerGDAL(): src=0,0,10474,10717, dst=27,63,140,61
msDrawRasterLayerGDAL(): source raster PL (-49.554,-33.774) for dst PL (27,63).
msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,2,3,4
msDrawMap(): Layer 0 (sat), 0.025s
msDrawMap(): Drawing Label Cache, 0.000s
msDrawMap() total time: 0.025s
msSaveImage(/data/temp/Jesper/ttt.png) total time: 0.004s
msFreeMap(): freeing map at 0x2453f50.
freeLayer(): freeing layer at 0x245c080.
freeLayer(): freeing layer at 0x245d430.
shp2img total time: 0.033s

Then the image is shown in ttt.png but is not in the correct place and not 
rotated correctly.

Any help or ideas are appreciated!

Thanks
Jesper

Den ons 2 sep. 2020 kl 18:23 skrev Jeff McKenna 
mailto:jmcke...@gatewaygeomatics.com>>:
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 
> 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,
> 

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

2020-09-03 Thread Jesper Kempe
shp2img creates an empty blank image with a mapfile like this:

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 15 60 30 68
SIZE 200 200

LAYER
  NAME "sat"
  DATA "/data/temp/Jesper/sat.tif"
  TYPE RASTER
  STATUS DEFAULT

  METADATA
"wms_title" "sat"
  END
END


The output from shp2img looks like this:

$ shp2img -m ttt.map -o /data/temp/Jesper/ttt.png -all_debug 5
msLoadMap(): 0.003s
msDrawMap(): rendering using outputformat named png (AGG/PNG).
msDrawRasterLayerLow(sat): entering.
msDrawRasterLayerGDAL(): Entering transform.
msDrawRasterLayerGDAL(): src=15,10645,15,15, dst=0,0,200,200
msDrawRasterLayerGDAL(): source raster PL (14.962,10645.462) for dst PL
(0,0).
msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,2,3,4
msDrawMap(): Layer 0 (sat), 0.011s
msDrawMap(): Drawing Label Cache, 0.000s
msDrawMap() total time: 0.011s
msSaveImage(/data/temp/Jesper/ttt.png) total time: 0.003s
msFreeMap(): freeing map at 0x1d14f50.
freeLayer(): freeing layer at 0x1d1d080.
freeLayer(): freeing layer at 0x1d1e1b0.
shp2img total time: 0.017s


If I add:
  PROJECTION
"init=epsg:4326"
  END
  EXTENT 17.0473154408269387 62.1606759790530177 27.5792082876237146
66.7744046276809797

shp2img outputs:

msLoadMap(): 0.003s
msDrawMap(): rendering using outputformat named png (AGG/PNG).
msDrawRasterLayerLow(sat): entering.
msDrawRasterLayerGDAL(): Entering transform.
msDrawRasterLayerGDAL(): src=0,0,10474,10717, dst=27,63,140,61
msDrawRasterLayerGDAL(): source raster PL (-49.554,-33.774) for dst PL
(27,63).
msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,2,3,4
msDrawMap(): Layer 0 (sat), 0.025s
msDrawMap(): Drawing Label Cache, 0.000s
msDrawMap() total time: 0.025s
msSaveImage(/data/temp/Jesper/ttt.png) total time: 0.004s
msFreeMap(): freeing map at 0x2453f50.
freeLayer(): freeing layer at 0x245c080.
freeLayer(): freeing layer at 0x245d430.
shp2img total time: 0.033s

Then the image is shown in ttt.png but is not in the correct place and not
rotated correctly.

Any help or ideas are appreciated!

Thanks
Jesper

Den ons 2 sep. 2020 kl 18:23 skrev Jeff McKenna <
jmcke...@gatewaygeomatics.com>:

> 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  > > 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) ->
> >