[gdal-dev] Image reprojection with GDALWARP

Rahkonen Jukka jukka.rahkonen at maanmittauslaitos.fi
Thu Feb 2 02:26:33 PST 2023


Hi,

There is a PROJ-data repository https://github.com/OSGeo/PROJ-data/ that contains lots of public grid-shift files, including https://github.com/OSGeo/PROJ-data/blob/master/es_ign/es_ign_SPED2ETV2.tif.

You can check from your local Proj data directory if the grid shift file is installed or not. For example OSGeo4W puts they info ...OSGeo4W\share\proj. But grid shift file is used automatically only when Proj thinks that it is used in the most suitable operation that it knows. In this case the operation is placed into 6th place as can be checked with the Proj utility "projinfo". It is possible to force the use of a preferred operation by using the Proj pipeline string.

projinfo -s EPSG:3857 -t EPSG:23030 --spatial-test intersects
...
Operation No. 6:

unknown id, Inverse of Popular Visualisation Pseudo-Mercator + Inverse of ED50 to WGS 84 (41) + UTM zone 30N, 1 m, Spain - mainland, Balearic Islands, Ceuta and Melila - onshore.

PROJ string:
+proj=pipeline
  +step +inv +proj=webmerc +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84
  +step +inv +proj=hgridshift +grids=es_ign_SPED2ETV2.tif
  +step +proj=utm +zone=30 +ellps=intl

WKT2:2019 string:
CONCATENATEDOPERATION["Inverse of Popular Visualisation Pseudo-Mercator + Inverse of ED50 to WGS 84 (41) + UTM zone 30N",
    SOURCECRS[
        PROJCRS["WGS 84 / Pseudo-Mercator",
            BASEGEOGCRS["WGS 84",
                ENSEMBLE["World Geodetic System 1984 ensemble",
                    MEMBER["World Geodetic System 1984 (Transit)"],
                    MEMBER["World Geodetic System 1984 (G730)"],
                    MEMBER["World Geodetic System 1984 (G873)"],
                    MEMBER["World Geodetic System 1984 (G1150)"],
                    MEMBER["World Geodetic System 1984 (G1674)"],
                    MEMBER["World Geodetic System 1984 (G1762)"],
                    MEMBER["World Geodetic System 1984 (G2139)"],
                    ELLIPSOID["WGS 84",6378137,298.257223563,
                        LENGTHUNIT["metre",1]],
                    ENSEMBLEACCURACY[2.0]],
                PRIMEM["Greenwich",0,
                    ANGLEUNIT["degree",0.0174532925199433]],
                ID["EPSG",4326]],
            CONVERSION["Popular Visualisation Pseudo-Mercator",
                METHOD["Popular Visualisation Pseudo Mercator",
                    ID["EPSG",1024]],
                PARAMETER["Latitude of natural origin",0,
                    ANGLEUNIT["degree",0.0174532925199433],
                    ID["EPSG",8801]],
                PARAMETER["Longitude of natural origin",0,
                    ANGLEUNIT["degree",0.0174532925199433],
                    ID["EPSG",8802]],
                PARAMETER["False easting",0,
                    LENGTHUNIT["metre",1],
                    ID["EPSG",8806]],
                PARAMETER["False northing",0,
                    LENGTHUNIT["metre",1],
                    ID["EPSG",8807]]],
            CS[Cartesian,2],
                AXIS["easting (X)",east,
                    ORDER[1],
                    LENGTHUNIT["metre",1]],
                AXIS["northing (Y)",north,
                    ORDER[2],
                    LENGTHUNIT["metre",1]],
            ID["EPSG",3857]]],
    TARGETCRS[
        PROJCRS["ED50 / UTM zone 30N",
            BASEGEOGCRS["ED50",
                DATUM["European Datum 1950",
                    ELLIPSOID["International 1924",6378388,297,
                        LENGTHUNIT["metre",1]]],
                PRIMEM["Greenwich",0,
                    ANGLEUNIT["degree",0.0174532925199433]],
                ID["EPSG",4230]],
            CONVERSION["UTM zone 30N",
                METHOD["Transverse Mercator",
                    ID["EPSG",9807]],
                PARAMETER["Latitude of natural origin",0,
                    ANGLEUNIT["degree",0.0174532925199433],
                    ID["EPSG",8801]],
                PARAMETER["Longitude of natural origin",-3,
                    ANGLEUNIT["degree",0.0174532925199433],
                    ID["EPSG",8802]],
                PARAMETER["Scale factor at natural origin",0.9996,
                    SCALEUNIT["unity",1],
                    ID["EPSG",8805]],
                PARAMETER["False easting",500000,
                   LENGTHUNIT["metre",1],
                    ID["EPSG",8806]],
                PARAMETER["False northing",0,
                    LENGTHUNIT["metre",1],
                    ID["EPSG",8807]]],
            CS[Cartesian,2],
                AXIS["(E)",east,
                    ORDER[1],
                    LENGTHUNIT["metre",1]],
                AXIS["(N)",north,
                    ORDER[2],
                    LENGTHUNIT["metre",1]],
            ID["EPSG",23030]]],
    STEP[
        CONVERSION["Inverse of Popular Visualisation Pseudo-Mercator",
            METHOD["Inverse of Popular Visualisation Pseudo Mercator",
                ID["INVERSE(EPSG)",1024]],
            PARAMETER["Latitude of natural origin",0,
                ANGLEUNIT["degree",0.0174532925199433],
                ID["EPSG",8801]],
            PARAMETER["Longitude of natural origin",0,
                ANGLEUNIT["degree",0.0174532925199433],
                ID["EPSG",8802]],
            PARAMETER["False easting",0,
                LENGTHUNIT["metre",1],
                ID["EPSG",8806]],
            PARAMETER["False northing",0,
                LENGTHUNIT["metre",1],
                ID["EPSG",8807]],
            ID["INVERSE(EPSG)",3856]]],
    STEP[
        COORDINATEOPERATION["Inverse of ED50 to WGS 84 (41)",
            SOURCECRS[
                GEOGCRS["WGS 84",
                    ENSEMBLE["World Geodetic System 1984 ensemble",
                        MEMBER["World Geodetic System 1984 (Transit)"],
                        MEMBER["World Geodetic System 1984 (G730)"],
                        MEMBER["World Geodetic System 1984 (G873)"],
                        MEMBER["World Geodetic System 1984 (G1150)"],
                        MEMBER["World Geodetic System 1984 (G1674)"],
                        MEMBER["World Geodetic System 1984 (G1762)"],
                        MEMBER["World Geodetic System 1984 (G2139)"],
                        ELLIPSOID["WGS 84",6378137,298.257223563,
                            LENGTHUNIT["metre",1]],
                        ENSEMBLEACCURACY[2.0]],
                    PRIMEM["Greenwich",0,
                        ANGLEUNIT["degree",0.0174532925199433]],
                    CS[ellipsoidal,2],
                        AXIS["geodetic latitude (Lat)",north,
                            ORDER[1],
                            ANGLEUNIT["degree",0.0174532925199433]],
                        AXIS["geodetic longitude (Lon)",east,
                            ORDER[2],
                            ANGLEUNIT["degree",0.0174532925199433]],
                    ID["EPSG",4326]]],
            TARGETCRS[
                GEOGCRS["ED50",
                    DATUM["European Datum 1950",
                        ELLIPSOID["International 1924",6378388,297,
                            LENGTHUNIT["metre",1]]],
                    PRIMEM["Greenwich",0,
                        ANGLEUNIT["degree",0.0174532925199433]],
                    CS[ellipsoidal,2],
                        AXIS["geodetic latitude (Lat)",north,
                            ORDER[1],
                            ANGLEUNIT["degree",0.0174532925199433]],
                        AXIS["geodetic longitude (Lon)",east,
                            ORDER[2],
                            ANGLEUNIT["degree",0.0174532925199433]],
                    ID["EPSG",4230]]],
            METHOD["Inverse of HORIZONTAL_SHIFT_GTIFF"],
            PARAMETERFILE["Latitude and longitude difference file","es_ign_SPED2ETV2.tif"],
            OPERATIONACCURACY[1.0],
            ID["INVERSE(DERIVED_FROM(EPSG))",15933],
            REMARK["Parameter values from ED50 to ETRS89 (12) (code 15932). Assumes ETRS89 and WGS 84 can be considered the same to within the accuracy of the transformation. Replaces ED50 to WGS 84 (40) - see supersession record."]]],
    STEP[
        CONVERSION["UTM zone 30N",
            METHOD["Transverse Mercator",
                ID["EPSG",9807]],
            PARAMETER["Latitude of natural origin",0,
                ANGLEUNIT["degree",0.0174532925199433],
                ID["EPSG",8801]],
            PARAMETER["Longitude of natural origin",-3,
                ANGLEUNIT["degree",0.0174532925199433],
                ID["EPSG",8802]],
            PARAMETER["Scale factor at natural origin",0.9996,
                SCALEUNIT["unity",1],
                ID["EPSG",8805]],
            PARAMETER["False easting",500000,
                LENGTHUNIT["metre",1],
                ID["EPSG",8806]],
            PARAMETER["False northing",0,
                LENGTHUNIT["metre",1],
                ID["EPSG",8807]],
            ID["EPSG",16030]]],
    USAGE[
        SCOPE["unknown"],
        AREA["Spain - mainland, Balearic Islands, Ceuta and Melila - onshore."],
        BBOX[35.26,-9.37,43.82,4.39]]]

-Jukka Rahkonen-

Lähettäjä: gdal-dev <gdal-dev-bounces at lists.osgeo.org> Puolesta Elena Ruiz
Lähetetty: torstai 2. helmikuuta 2023 11.44
Vastaanottaja: Javier Jimenez Shaw <j1 at jimenezshaw.com>
Kopio: gdal-dev at lists.osgeo.org
Aihe: Re: [gdal-dev] Image reprojection with GDALWARP

Hi Javier,
It is true that the ED50 is a very old system, and that it is hardly used, but I came across this case and the lack of precision shocked me (they are more than 9 m), I know the grid that you suggest ((es_ign_SPED2ETV2.tif)) and it is the one that is currently used with everything related to ED50, but I assumed that this was already taken into account by GDAL when indicating the EPSG:23030 system. Is there an option in the GDAL configuration where you can indicate where to get those grids from? If it happens with ED50, it will also happen with other systems, right? How to tell GDAl or the gdalwarp command to choose the most correct or precise coordinate transformation? Greetings and thank you




________________________________

Elena Ruiz

[https://www.aplitop.com/descargas/logo.png]
Sofware Development & Technical Support
Tel.: +34 952 43 97 71
eruiz at aplitop.com<mailto:eruiz at aplitop.com>
www.aplitop.com<http://www.aplitop.com>










De: Javier Jimenez Shaw <j1 at jimenezshaw.com<mailto:j1 at jimenezshaw.com>>
Enviado el: miércoles, 1 de febrero de 2023 18:57
Para: Elena Ruiz <eruiz at aplitop.com<mailto:eruiz at aplitop.com>>
CC: Even Rouault <even.rouault at spatialys.com<mailto:even.rouault at spatialys.com>>; gdal-dev at lists.osgeo.org<mailto:gdal-dev at lists.osgeo.org>
Asunto: Re: [gdal-dev] Image reprojection with GDALWARP

Hi Elena

I see that you are converting to ED50, a very old system. I assume that you have to do it, instead of ETRS89 (I'm sorry about that).
How big is the reprojection error? I do not find that in your previous email. To get a proper transformation you need a grid file (es_ign_SPED2ETV2.tif<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fcdn.proj.org%2Fes_ign_SPED2ETV2.tif&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C2da75d97da044846b7fb08db050212f9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638109278711733097%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000%7C%7C%7C&sdata=%2FeauKHkk8tMTyV5FLvTMqSj7KHmDEvjQRR9QNjc3l24%3D&reserved=0>). I do not know if it is supported in gdal 2.4.4
The transformation from ED50 to ETRS89 is defined in different ways in different locations.

Cheers,
Javier
.___ ._ ..._ .. . ._.  .___ .. __ . _. . __..  ... .... ._ .__
Entre dos pensamientos racionales
hay infinitos pensamientos irracionales.


On Wed, 1 Feb 2023 at 18:06, Elena Ruiz <eruiz at aplitop.com<mailto:eruiz at aplitop.com>> wrote:
Hello, thank you very much for answering, I have added the "-r bilinear" option to the gdalwarp command and this solves the distortion problem, but the reprojection is still wrong with the coordinates of the drawing matching those of the image, and the drawing does not fit the image as it does the original image. How do I have to call gdalwarp so that it does the reprojection correctly, that is, in addition to modifying the image that changes the coordinates and the image projection system? Regards




________________________________

Elena Ruiz


Sofware Development & Technical Support
Tel.: +34 952 43 97 71
eruiz at aplitop.com<mailto:eruiz at aplitop.com>
www.aplitop.com<https://eur06.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.aplitop.com%2F&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C2da75d97da044846b7fb08db050212f9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638109278711733097%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000%7C%7C%7C&sdata=dPSweDSS2VuoGGcPQjP4nNEI3gL0AItdPZMUsPJeF%2Bw%3D&reserved=0>










De: Even Rouault <even.rouault at spatialys.com<mailto:even.rouault at spatialys.com>>
Enviado el: miércoles, 1 de febrero de 2023 16:58
Para: Elena Ruiz <eruiz at aplitop.com<mailto:eruiz at aplitop.com>>; gdal-dev at lists.osgeo.org<mailto:gdal-dev at lists.osgeo.org>
Asunto: Re: [gdal-dev] Image reprojection with GDALWARP


Try adding -r bilinear or -r bicubic. By default gdalwarp performs nearest resampling which is rarely what you need for high quality results
Le 01/02/2023 à 16:44, Elena Ruiz a écrit :

Hi, I'm trying to reproject images from WGS 84 / Pseudo-Mercator EPSG:3857 to ED50 / UTM zone 30N EPSG:23030, but the resulting image doesn't adapt well to drawing. I will explain it in more detail so that you understand:

I have an image in EPSG:3857, the result of a WMTS service, and a drawing in the same SRS, as you can see in the image the drawing matches the image perfectly.
Image source with drawing: https://aplitop.sharefile.com/d-s4d6dcf6d8ff34023a8b6e95277519c56<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Faplitop.sharefile.com%2Fd-s4d6dcf6d8ff34023a8b6e95277519c56&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C2da75d97da044846b7fb08db050212f9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638109278711733097%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000%7C%7C%7C&sdata=9HlbKtVRg8lkgeXRjWF2yrQMFeVCYlxY1Ktdg1h164A%3D&reserved=0>

I want to convert the previous image to EPSG:23030 and I have used gdalwarp command with different parameters and I get the same result :

1.gdalwarp -s_srs EPSG:3857 -t_srs EPSG:23030 -co "TFW=YES" -overwrite "C:\Users\elena\Downloads\Nueva carpeta\imagen_t.tif" "C:\Users\elena\Downloads\Nueva carpeta\nueva_6.tif"
2.gdalwarp -s_srs EPSG:3857 -t_srs EPSG:23030 -te 335000 4084501.25 335550 4084848.75 -co "TFW=YES" -overwrite "C:\Users\elena\Downloads\Nueva carpeta\imagen_t.tif" "C:\Users\elena\Downloads\Nueva carpeta\nueva_6.tif"

And the resulting image, in addition to noticing that the road is distorted (like a staircase effect), does not match the drawing.
Image destination with drawing: https://aplitop.sharefile.com/d-s98f76513feb34c268764c0492ccb40d3<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Faplitop.sharefile.com%2Fd-s98f76513feb34c268764c0492ccb40d3&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C2da75d97da044846b7fb08db050212f9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638109278711733097%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000%7C%7C%7C&sdata=0SEHZyxLt3BCju9dsOGb1Wkf0Zg8UPRB1RH1q1bDGj0%3D&reserved=0>
Zoom of image destination: https://aplitop.sharefile.com/d-sf9a270a9300e4ac796eddeac7131a058<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Faplitop.sharefile.com%2Fd-sf9a270a9300e4ac796eddeac7131a058&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C2da75d97da044846b7fb08db050212f9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638109278711889331%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000%7C%7C%7C&sdata=KWIdNuSi8aKeWfZm4p6%2B8ZDNd3I9xc1zj%2FIqI9nA%2B2c%3D&reserved=0>

I have also done the test from EPSG:23030 to EPSG:25830 and it does not do well either.
I am using the version of gdal 2.4.4 for x64, this is the link to download the source and destination image: https://aplitop.sharefile.com/d-sdae2552bd2054fdf8cc29ce5b559e278<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Faplitop.sharefile.com%2Fd-sdae2552bd2054fdf8cc29ce5b559e278&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C2da75d97da044846b7fb08db050212f9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638109278711889331%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000%7C%7C%7C&sdata=aYb5Wig9wrJ13soJaCSfCDR0Kaz4nT5h3zzGTrw%2F%2F%2Bs%3D&reserved=0> , I hope someone can help me and tell me if there is something wrong. Regards
________________________________
Elena Ruiz
Sofware Development & Technical Support
Tel. +34 952 43 97 71
eruiz at aplitop.com<mailto:eruiz at aplitop.com>
Sumatra, 9 - 29190 Málaga (Spain)
www.aplitop.com<https://eur06.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.aplitop.com%2F&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C2da75d97da044846b7fb08db050212f9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638109278711889331%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000%7C%7C%7C&sdata=i%2BPF%2F84lMX0pYl06w0BXcpP%2BTQhh4diWuSqbGXPnmUA%3D&reserved=0>

[cid:image001.png at 01D936FF.CB24B080]<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.facebook.com%2FAplitopSL%2F&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C2da75d97da044846b7fb08db050212f9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638109278711889331%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000%7C%7C%7C&sdata=IQ4ftCJ2MBnpe9Qf1jPthGWLVEmp9fEIeoF2uxaaAe4%3D&reserved=0>   [cid:image002.png at 01D936FF.CB24B080] <https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Ftwitter.com%2Faplitop&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C2da75d97da044846b7fb08db050212f9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638109278711889331%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000%7C%7C%7C&sdata=Ad%2F3FdKYOcbNN8Uz4zNl3VOGDhwPyV4AEy0Ata2HORw%3D&reserved=0>    [cid:image003.png at 01D936FF.CB24B080] <https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.linkedin.com%2Fcompany%2Faplitop&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C2da75d97da044846b7fb08db050212f9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638109278711889331%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000%7C%7C%7C&sdata=MUf7Gvq8nC23izMzbD2U1tSnL10ygIAiOrG7hcyjxz8%3D&reserved=0>    [cid:image004.png at 01D936FF.CB24B080] <https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.youtube.com%2Fuser%2Faplitop&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C2da75d97da044846b7fb08db050212f9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638109278711889331%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000%7C%7C%7C&sdata=xLPilUsHQJOXoIy1Im54sBE2QV5DSKyZH9MAlybbFng%3D&reserved=0>

In accordance with the provisions of the European Regulation of Data Protection 2016/679 (Reglamento Europeo de Protección de Datos 2016/679), we inform you that the data and the information you provide us through this medium will be used by APLITOP, S.L., with C.I.F. B-92543396 and with address at C / Sumatra, 9, Malaga, 29190, in order to answer your questions and inform you about our products. The data provided will be kept as long as it does not request its cessation and will not be transferred to third parties except in cases where there is a legal obligation. You have the right to access your personal data, correct inaccurate data or request its deletion when the data is no longer necessary for the purposes that were collected, as well as any rights recognized in the RGPD 2016/679




_______________________________________________

gdal-dev mailing list

gdal-dev at lists.osgeo.org<mailto:gdal-dev at lists.osgeo.org>

https://lists.osgeo.org/mailman/listinfo/gdal-dev<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.osgeo.org%2Fmailman%2Flistinfo%2Fgdal-dev&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C2da75d97da044846b7fb08db050212f9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638109278711889331%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000%7C%7C%7C&sdata=zW8FiL6InBqWhIGEnWNxe99276NlbG8CT0uPCLs8r10%3D&reserved=0>

--

http://www.spatialys.com<https://eur06.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.spatialys.com%2F&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C2da75d97da044846b7fb08db050212f9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638109278711889331%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000%7C%7C%7C&sdata=wvxE%2FRf21QBgmBUWlKxRf3%2BMOS3%2FaNqL2d1XuVhcRkU%3D&reserved=0>

My software is free, but my time generally not.
_______________________________________________
gdal-dev mailing list
gdal-dev at lists.osgeo.org<mailto:gdal-dev at lists.osgeo.org>
https://lists.osgeo.org/mailman/listinfo/gdal-dev<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.osgeo.org%2Fmailman%2Flistinfo%2Fgdal-dev&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7C2da75d97da044846b7fb08db050212f9%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C638109278711889331%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000%7C%7C%7C&sdata=zW8FiL6InBqWhIGEnWNxe99276NlbG8CT0uPCLs8r10%3D&reserved=0>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20230202/01a7cea9/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.png
Type: image/png
Size: 365 bytes
Desc: image001.png
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20230202/01a7cea9/attachment-0004.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image002.png
Type: image/png
Size: 619 bytes
Desc: image002.png
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20230202/01a7cea9/attachment-0005.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image003.png
Type: image/png
Size: 468 bytes
Desc: image003.png
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20230202/01a7cea9/attachment-0006.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image004.png
Type: image/png
Size: 817 bytes
Desc: image004.png
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20230202/01a7cea9/attachment-0007.png>


More information about the gdal-dev mailing list