<div dir="ltr">Hi Jukka<div><br><div>Sorry about that, I forgot that yesterday.</div></div><div><br></div><div>I did a test now and that actually works!</div><div><br></div><div>Changing DATA to the vrt-file and setting </div><div><br></div><div>PROJECTION</div><div> "init=epsg:3857"<br></div><div>END</div><div><br></div><div>creates an image with sat.tif placed correctly!</div><div><br></div><div>$ shp2img -m icesat.map -o /data/temp/Jesper/ttt.png -all_debug 5<br>msLoadMap(): 0.014s<br>msDrawMap(): rendering using outputformat named png (AGG/PNG).<br>msDrawRasterLayerLow(sat): entering.<br>msResampleGDALToMap in effect: cellsize = 4202.211585<br>msDrawGDAL(sat): using RAW_WINDOW=0 0 12531 12770, dst=0,0,279,284<br>msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,2,3,4<br>msDrawMap(): Layer 0 (sat), 557.242s<br>msDrawMap(): Drawing Label Cache, 0.000s<br>msDrawMap() total time: 557.243s<br>msSaveImage(/data/temp/Jesper/ttt.png) total time: 0.119s<br>msFreeMap(): freeing map at 0x7b1f50.<br>freeLayer(): freeing layer at 0x7ba080.<br>freeLayer(): freeing layer at 0x7bb760.<br>shp2img total time: 557.378s<br></div><div><br></div><div>BUT, note 557.378s on the last row...</div><div><br></div><div>Also a WMS-request works but takes about 4 minutes to complete!</div><div><br></div><div>It looks like I need to reproject all files new GeoTIFFs to get a decent speed.</div><div><br></div><div>Thanks for all the replies!</div><div><br></div><div>Jesper</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Den tors 3 sep. 2020 kl 13:31 skrev Rahkonen Jukka (MML) <<a href="mailto:jukka.rahkonen@maanmittauslaitos.fi">jukka.rahkonen@maanmittauslaitos.fi</a>>:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div lang="FI">
<div class="gmail-m_-163699377220558419WordSection1">
<p class="MsoNormal"><span>Hi,<u></u><u></u></span></p>
<p class="MsoNormal"><span><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">I would appreciate if you could have a try with<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">"gdalwarp -of vrt -t_srs "EPSG:3857" sat.tif sat_3857.vrt"</span><span lang="EN-US"><u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-US">-Jukka Rahkonen-<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p>
<p class="MsoNormal"><b>Lähettäjä:</b> mapserver-users <<a href="mailto:mapserver-users-bounces@lists.osgeo.org" target="_blank">mapserver-users-bounces@lists.osgeo.org</a>>
<b>Puolesta </b>Jesper Kempe<br>
<b>Lähetetty:</b> torstai 3. syyskuuta 2020 13.57<br>
<b>Vastaanottaja:</b> Jeff McKenna <<a href="mailto:jmckenna@gatewaygeomatics.com" target="_blank">jmckenna@gatewaygeomatics.com</a>><br>
<b>Kopio:</b> Mapserver <<a href="mailto:mapserver-users@lists.osgeo.org" target="_blank">mapserver-users@lists.osgeo.org</a>><br>
<b>Aihe:</b> Re: [mapserver-users] Use Ground Control Points from GeoTIFF<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">shp2img creates an empty blank image with a mapfile like this:<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">MAP<br>
NAME "map"<br>
WEB<br>
METADATA<br>
"wms_title" "WMS mapserver"<br>
"wms_srs" "EPSG:4326 EPSG:3021 EPSG:9999 EPSG:3006 EPSG:900913"<br>
"wms_extent" "1225000 6130000 1910000 7675000"<br>
"wms_enable_request" "*"<br>
"wfs_enable_request" "*"<br>
END<br>
END<br>
<br>
PROJECTION<br>
"init=epsg:4326"<br>
END<br>
<br>
EXTENT 15 60 30 68<br>
SIZE 200 200<br>
<br>
LAYER<br>
NAME "sat"<br>
DATA "/data/temp/Jesper/sat.tif"<br>
TYPE RASTER<br>
STATUS DEFAULT<br>
<br>
METADATA<br>
"wms_title" "sat"<br>
END<br>
END<u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">The output from shp2img looks like this:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">$ shp2img -m ttt.map -o /data/temp/Jesper/ttt.png -all_debug 5<br>
msLoadMap(): 0.003s<br>
msDrawMap(): rendering using outputformat named png (AGG/PNG).<br>
msDrawRasterLayerLow(sat): entering.<br>
msDrawRasterLayerGDAL(): Entering transform.<br>
msDrawRasterLayerGDAL(): src=15,10645,15,15, dst=0,0,200,200<br>
msDrawRasterLayerGDAL(): source raster PL (14.962,10645.462) for dst PL (0,0).<br>
msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,2,3,4<br>
msDrawMap(): Layer 0 (sat), 0.011s<br>
msDrawMap(): Drawing Label Cache, 0.000s<br>
msDrawMap() total time: 0.011s<br>
msSaveImage(/data/temp/Jesper/ttt.png) total time: 0.003s<br>
msFreeMap(): freeing map at 0x1d14f50.<br>
freeLayer(): freeing layer at 0x1d1d080.<br>
freeLayer(): freeing layer at 0x1d1e1b0.<br>
shp2img total time: 0.017s<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">If I add:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> PROJECTION<br>
"init=epsg:4326"<br>
END<br>
EXTENT 17.0473154408269387 62.1606759790530177 27.5792082876237146 66.7744046276809797<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">shp2img outputs:<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">msLoadMap(): 0.003s<br>
msDrawMap(): rendering using outputformat named png (AGG/PNG).<br>
msDrawRasterLayerLow(sat): entering.<br>
msDrawRasterLayerGDAL(): Entering transform.<br>
msDrawRasterLayerGDAL(): src=0,0,10474,10717, dst=27,63,140,61<br>
msDrawRasterLayerGDAL(): source raster PL (-49.554,-33.774) for dst PL (27,63).<br>
msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,2,3,4<br>
msDrawMap(): Layer 0 (sat), 0.025s<br>
msDrawMap(): Drawing Label Cache, 0.000s<br>
msDrawMap() total time: 0.025s<br>
msSaveImage(/data/temp/Jesper/ttt.png) total time: 0.004s<br>
msFreeMap(): freeing map at 0x2453f50.<br>
freeLayer(): freeing layer at 0x245c080.<br>
freeLayer(): freeing layer at 0x245d430.<br>
shp2img total time: 0.033s<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Then the image is shown in ttt.png but is not in the correct place and not rotated correctly.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Any help or ideas are appreciated!<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Thanks<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Jesper<u></u><u></u></p>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<p class="MsoNormal">Den ons 2 sep. 2020 kl 18:23 skrev Jeff McKenna <<a href="mailto:jmckenna@gatewaygeomatics.com" target="_blank">jmckenna@gatewaygeomatics.com</a>>:<u></u><u></u></p>
</div>
<blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin-left:4.8pt;margin-right:0cm">
<p class="MsoNormal">I'm curious what happens when you try to generate a map image from the
<br>
commandline: shp2img -m test.map -o ttt.png -all_debug 5<br>
<br>
-jeff<br>
<br>
<br>
<br>
-- <br>
Jeff McKenna<br>
MapServer Consulting and Training Services<br>
co-founder of FOSS4G<br>
<a href="http://gatewaygeo.com/" target="_blank">http://gatewaygeo.com/</a><br>
<br>
<br>
<br>
On 2020-09-02 4:51 a.m., Jesper Kempe wrote:<br>
> Hi<br>
> <br>
> This is the mapfile:<br>
> <br>
> MAP<br>
> NAME "map"<br>
> WEB<br>
> METADATA<br>
> "wms_title" "WMS mapserver"<br>
> "wms_srs" "EPSG:4326 EPSG:3021 EPSG:9999 <br>
> EPSG:3006 EPSG:900913"<br>
> "wms_extent" "1225000 6130000 1910000 7675000"<br>
> "wms_enable_request" "*"<br>
> "wfs_enable_request" "*"<br>
> END<br>
> END<br>
> <br>
> PROJECTION<br>
> "init=epsg:4326"<br>
> END<br>
> <br>
> EXTENT -180 -90 180 90<br>
> <br>
> LAYER<br>
> NAME "sat"<br>
> DATA "sat.tif"<br>
> TYPE RASTER<br>
> STATUS ON<br>
> <br>
> METADATA<br>
> "wms_title" "sat"<br>
> END<br>
> <br>
> PROJECTION<br>
> "init=epsg:4326"<br>
> END<br>
> EXTENT 17.0473154408269387 62.1606759790530177 <br>
> 27.5792082876237146 66.7744046276809797<br>
> <br>
> END<br>
> <br>
> END<br>
> <br>
> I have also tried both without PROJECTION in LAYER and with<br>
> PROJECTION<br>
> AUTO<br>
> END<br>
> <br>
> None of them works.<br>
> <br>
> EXTENT in the mapfile is taken from Qgis metadata of the layer.<br>
> <br>
> <br>
> Thanks<br>
> Jesper<br>
> <br>
> Den tis 1 sep. 2020 kl 23:41 skrev Steve Lime <<a href="mailto:sdlime@gmail.com" target="_blank">sdlime@gmail.com</a>
<br>
> <mailto:<a href="mailto:sdlime@gmail.com" target="_blank">sdlime@gmail.com</a>>>:<br>
> <br>
> What does the mapfile look like?<br>
> <br>
> On Tue, Sep 1, 2020 at 10:26 AM Jesper Kempe <<a href="mailto:jesper.kempe@gmail.com" target="_blank">jesper.kempe@gmail.com</a><br>
> <mailto:<a href="mailto:jesper.kempe@gmail.com" target="_blank">jesper.kempe@gmail.com</a>>> wrote:<br>
> <br>
> Hi<br>
> <br>
> I have a GeoTiff file that gdalinfo presents like this:<br>
> <br>
> gdalinfo sat.tif<br>
> Driver: GTiff/GeoTIFF<br>
> Files: sat.tif<br>
> Size is 10474, 10717<br>
> Coordinate System is `'<br>
> GCP Projection =<br>
> GEOGCS["WGS 84",<br>
> DATUM["WGS_1984",<br>
> SPHEROID["WGS 84",6378137,298.257223563,<br>
> AUTHORITY["EPSG","7030"]],<br>
> AUTHORITY["EPSG","6326"]],<br>
> PRIMEM["Greenwich",0],<br>
> UNIT["degree",0.0174532925199433],<br>
> AUTHORITY["EPSG","4326"]]<br>
> GCP[ 0]: Id=1, Info=<br>
> (0,0) -><br>
> (27.5783052509812,65.9184735836942,229.992792810313)<br>
> GCP[ 1]: Id=2, Info=<br>
> (524,0) -><br>
> (27.1416007200269,65.9731480840222,172.994255514815)<br>
> GCP[ 2]: Id=3, Info=<br>
> (1048,0) -><br>
> (26.7005564512604,66.0268770755367,152.994545626454)<br>
> GCP[ 3]: Id=4, Info=<br>
> (1572,0) -><br>
> (26.257537545879,66.0793551447287,131.994995356537)<br>
> GCP[ 4]: Id=5, Info=<br>
> (2096,0) -><br>
> (25.8102107755995,66.1308402121171,163.993212723173)<br>
> ....<br>
> <br>
> with 482 GCPs.<br>
> <br>
> <br>
> The file opens ok in Qgis and is located in the right location.<br>
> But I can´t in any way get it to show correctly in Mapserver.<br>
> <br>
> Is this possible to do? Or do I have to gdalwarp the file to a<br>
> known projection and set the PROJECTION of the layer?<br>
> <br>
> Thanks<br>
> Jesper<br>
> <br>
> <br>
<br>
<br>
<br>
_______________________________________________<br>
mapserver-users mailing list<br>
<a href="mailto:mapserver-users@lists.osgeo.org" target="_blank">mapserver-users@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/mapserver-users" target="_blank">https://lists.osgeo.org/mailman/listinfo/mapserver-users</a><u></u><u></u></p>
</blockquote>
</div>
</div>
</div>
</blockquote></div>