Global Mapper Pro

Downloading from WMS Geoportail.fr server

falqnfalqn Global Mapper UserTrusted User Posts: 123
edited April 2015 in Downloading Imagery
Hi,
I have problems with downloading data from French geoportail.fr server http://wxs.ign.fr/geoportail/wmsc - list of layers appear but GM can not download any portion of data...
When I try with http://wxs.ign.fr/geoportail/wms it's working with one available layer - imagery. Is WMS-C service available in GM? What is the difference between WMS and WMS-C?
more info www.geoportail.fr and https://api.ign.fr/geoportail/api/doc/index.html

Best regards
Wojtek, Poland
«1

Comments

  • global_mapperglobal_mapper Administrator Posts: 17,238
    edited March 2010
    Wojtek,

    I just tried and am getting an error return of 403 from the server, which means that access is forbidden. Is this server perhaps locked down so that your IP address has to be approved by them before they will server up the data?

    Thanks,

    Mike
    Global Mapper Support
    support@globalmapper.com
  • GeoGeo Global Mapper User Trusted User Posts: 92
    edited March 2010
    falqn wrote: »

    The 2D API allows you to create maps by getting connected to Geoportal's services.
    Geoportal's WMS-C, WMS, WFS services can be plugged to desktop GIS applications. One needs to develop a connector to our GeoRM implementation along with the above services connectors. For more information, please contact us at 'contact dot api at ign dot fr'.

    A proof of concept of using Geoportal's 2D streams into 3D application has been to use Nasa World Wind
    This application is based on the NASA World Wind Java Development Kit to allow Geoportal's layers overlaying in this virtual globe.
    These Geoportal's layers are served by the WMS-C service as explained there.
    Click here to launch the application.


    Geoportal API Documentation in PDF
  • papaverpapaver Global Mapper User Trusted User Posts: 25
    edited January 2011
    Hi,

    Some informations about Geoportail and WMS-C : https://api.ign.fr/geoportail/api/doc/developpeur/geodrm.html

    Geoportail is a great data source for France... but impossible to import any data in GM at this time. Any solution in the future ?
  • global_mapperglobal_mapper Administrator Posts: 17,238
    edited January 2011
    I'll take a look when I get a chance and see if it would be simple to add support for these layers in Global Mapper.

    Thanks,

    Mike
    Global Mapper Support
    support@globalmapper.com
  • fritzpeterfritzpeter Global Mapper User Posts: 4
    edited June 2011
    Hi,
    There are news about the problem of this thread?
    The http://wxs.ign.fr/geoportail/wms in my case don't work, list of layers appear but I can't download from FXX.ORTHOIMAGERY.ORTHOPHOTOS
    Suggestions?

    Best regards, Fritz Peter
  • global_mapperglobal_mapper Administrator Posts: 17,238
    edited June 2011
    Fritz,

    Sorry I haven't had a chance to look at adding WMS-C support yet to know how much work it is.

    Thanks,

    Mike
    Global Mapper Support
    support@globalmapper.com
  • fritzpeterfritzpeter Global Mapper User Posts: 4
    edited July 2011
    Hi Mike,
    thanks a lot for your reply. From my point of view, would be really important to add this feature to Globalmapper, that is a complete and wonderful software, perhaps the best in its kind.

    Best regards. Fritz.
  • fritzpeterfritzpeter Global Mapper User Posts: 4
    edited September 2011
    Hi Mike,
    with the new version of Globalmapper the WMS-C support is been added?

    Best regards. Fritz
  • global_mapperglobal_mapper Administrator Posts: 17,238
    edited September 2011
    There isn't any support for WMS-C yet.

    Thanks,

    Mike
    Global Mapper Support
    support@globalmapper.com
  • fritzpeterfritzpeter Global Mapper User Posts: 4
    edited September 2011
    Thanks.

    Fritz.
  • GeoGeo Global Mapper User Trusted User Posts: 92
    edited September 2011
    hello,

    with the informations given at this url
    http://api.ign.fr/geoportail/document.do?doc=qgis
    and
    https://api.ign.fr/geoportail/document.do?doc=6593220
    
    "INSPIRE WMS" is working with Global Mappper, but it's slow, and there are less layers available.

    "INSPIRE WMS-c" could work but the parameters "&WIDTH=512&HEIGHT=512" are always sent by Global Mapper with "&REQUEST=GetMap" so "&BBOX=" is always wrong because the server accepts only "&WIDTH=256&HEIGHT=256".

    according to
    GET http://wxs.ign.fr/inspire/wmsc?ServiceName=WMS&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetCapabilities
    

    the <TileSet> parameter describe for each available layer this parameters :
    • <SRS>
    • <BoundingBox>
    • <Resolutions>
    • <Width>
    • <Height>
    • <Format>
    but Global Mapper seems to not use all of them, especially <Width>, <Height> and <Format>



    how to reproduce the errors :
    1. create an account at
      https://api.ign.fr/geoportail/myAccount.do
      
      and get a "service de consultation INSPIRE" key
    2. download in an unique directory, the GeoDRM wrapper, the configuration file and the batch file.
      https://api.ign.fr/geoportail/api/doc/sources/geodrm-wrapper.jar
      http://www.ign.fr/partage/api/tutoriel/wrapper-config-wmsc.xml
      http://www.ign.fr/partage/api/tutoriel/qgis_proxy_geoportail.bat
      
    3. edit "wrapper-config-wmsc.xml" and replace XXXXXXXXXXXXXXXXXXX with your API Key
      <?xml version="1.0"?>
      <wrapper>
        <!--
        <proxy>
          <host></host>
          <port></port>
        </proxy> -->
        <remoteServices>
          <remoteService>
            <localHost>localhost</localHost>
            <localPort>10001</localPort>
            <localPath>/wmsc</localPath>
            <remoteHost>wxs.ign.fr</remoteHost>
            <remotePort>80</remotePort>
            <remotePath>/inspire/wmsc</remotePath>
            <tokenProviderUrl>http://jeton-api.ign.fr/getToken</tokenProviderUrl>
            <apiKey>XXXXXXXXXXXXXXXXXXX</apiKey>
            <tokenLifeTime>600</tokenLifeTime>
          </remoteService>
        </remoteServices>
          <remoteServices>
          <remoteService>
            <localHost>localhost</localHost>
            <localPort>10002</localPort>
            <localPath>/wms</localPath>
            <remoteHost>wxs.ign.fr</remoteHost>
            <remotePort>80</remotePort>
            <remotePath>/inspire/wms</remotePath>
            <tokenProviderUrl>http://jeton-api.ign.fr/getToken</tokenProviderUrl>
            <apiKey>XXXXXXXXXXXXXXXXXXX</apiKey>
            <tokenLifeTime>600</tokenLifeTime>
          </remoteService>
        </remoteServices>
      </wrapper>
      
    4. for debugging purpose, i used Paros Proxy (free of charge man-in-the-middle proxy which captures all traffic between clients and servers)
      http://www.parosproxy.org/
      
      and updated the proxy section of "wrapper-config-wmsc.xml" with
        <proxy>
          <host>localhost</host>
          <port>8080</port>
        </proxy>
      
    5. if debugging, launch "paros" then "qgis_proxy_geoportail.bat"
    6. if not, "qgis_proxy_geoportail.bat" only from the console (if configuration errors) or with a double clic.
    7. edit [drive:]\[installation_folder_of_Global Mapper]\epsg_codes.txt and add this lines at he bottom
      //
      //IGNF:GEOPORTAL
      310024802,Equirectangular,RGF-93,meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=46.5,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
      
      there should be other EPSG code to be added, but this one is for France Metropolitan (FXX)
    8. open Global Mapper
    9. add a new WMS source
      http://127.0.0.1:10001/wmsc
      WMS
      
      get the list of available layers and choose one tagged with (FXX)
    10. an error message will popup at the end
      if debugging, Paros will show all the requests and Responses
    11. NB : with Paros / Tools / Manual Request Editor, you can build your own request
    12. add a new WMS source
      http://127.0.0.1:10002/wms
      WMS
      
      get the list of available layers and choose one tagged with FXX
    13. images will be loaded in Global Mapper


    it would be great if you could add support for "INSPIRE WMS_c" in Global Mapper.
    with the GeoDRM wrapper, i think it's very close to the end, with your support of course ;)

    thank you for your attention,
    Richard.
  • global_mapperglobal_mapper Administrator Posts: 17,238
    edited September 2011
    Richard,

    I will take another look and see if support for that particular source can easily be added.

    Thanks,

    Mike
    Global Mapper Support
    support@globalmapper.com
  • global_mapperglobal_mapper Administrator Posts: 17,238
    edited September 2011
    Richard,

    It's not quite a simple matter of switching to use to use a tile size of 256. Global Mapper needs to read the <TileSet> tags and setup its tiling of the WMS data based on the tile layout described in that rather than using its own default access zoom levels and tile sizes. Not trivial but certainly not rocket science either. I'll see if I can take a look in the next week or so and make this happen.

    Thanks,

    Mike
    Global Mapper Support
    support@globalmapper.com
  • GeoGeo Global Mapper User Trusted User Posts: 92
    edited September 2011
    Mike,

    you're right, WMS_c needs that all request have to be made with pre-defined parameters and ask for a pre-existing tile.
    i've found a WMS Tiling Client Recommendation documentation here :
    http://wiki.osgeo.org/wiki/WMS_Tiling_Client_Recommendation
    

    there is also a small documentation in english on how to build a WMS-C request for the Geoportail.
    see : API Geoportail / Documentation / Advanced WMS-C use
    https://api.ign.fr/geoportail/api/doc/developpeur/wmsc.html
    

    the full english Geoportail API Documentation is also available in PDF
    https://api.ign.fr/geoportail/api/doc/GeoportalAPI.pdf
    
    see chapter 17, Advanced WMS-C layers and the internal link for layer resolutions.

    hope this help,
    thanks for your support,
    Richard.
  • GeoGeo Global Mapper User Trusted User Posts: 92
    edited September 2011
    hello,

    i made a mistake in the epsg_codes.txt
    //
    //IGNF:GEOPORTAL
    310024802,Equirectangular,RGF-93,meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=46.5,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    
    with this parameters, it produces a vertical shift and a horizontal shift.

    i can't figure out the good parameters for the epsg code in Global Mapper.

    test procedure :
    i retrieved from the Geoportail wmsc server an area with pixel size=0.25m ; i used Gdal and a file called cad025.xml
    <GDAL_WMS>
    <!--cad025.xml-->
    	<Service name="WMS">
    		<Version>1.1.1</Version>
    		<ServerUrl>http://127.0.0.1:10001/wmsc?SERVICE=WMS&TILED=true</ServerUrl>
    		<SRS>EPSG:310024802</SRS>
    		<Layers>CADASTRALPARCELS.PARCELS</Layers>
    		<ImageFormat>image/png</ImageFormat>
    	</Service>
    	<DataWindow>
    		<UpperLeftX>541952</UpperLeftX>
    		<LowerRightX>542976</LowerRightX>
    		<UpperLeftY>5448192</UpperLeftY>
    		<LowerRightY>5447168</LowerRightY>
    		<TileLevel>4</TileLevel>
    	</DataWindow>
    	<BlockSizeX>256</BlockSizeX>
    	<BlockSizeY>256</BlockSizeY>
    	<BandsCount>1</BandsCount>
    	<MaxConnections>5</MaxConnections>
    	<Timeout>300</Timeout>
    </GDAL_WMS>
    


    command : gdalinfo cad025.xml
    Driver: WMS/OGC Web Map Service
    Files: cad025.xml
    Size is 4096, 4096
    Coordinate System is:
    PROJCS["unnamed",
        GEOGCS["unnamed ellipse",
            DATUM["unknown",
                SPHEROID["unnamed",6378137,298.257222101],
                EXTENSION["PROJ4_GRIDS","null"]],
            PRIMEM["Greenwich",0],
            UNIT["degree",0.0174532925199433]],
        PROJECTION["Equirectangular"],
        PARAMETER["latitude_of_origin",0],
        PARAMETER["central_meridian",0],
        PARAMETER["standard_parallel_1",46.5],
        PARAMETER["false_easting",0],
        PARAMETER["false_northing",0],
        UNIT["Meter",1],
        AUTHORITY["EPSG","310024802"]]
    Origin = (541952.000000000000000,5448192.000000000000000)
    Pixel Size = (0.250000000000000,-0.250000000000000)
    Corner Coordinates:
    Upper Left  (  541952.000, 5448192.000) (  7d 4'21.26"E, 48d56'30.99"N)
    Lower Left  (  541952.000, 5447168.000) (  7d 4'21.26"E, 48d55'57.87"N)
    Upper Right (  542976.000, 5448192.000) (  7d 5' 9.37"E, 48d56'30.99"N)
    Lower Right (  542976.000, 5447168.000) (  7d 5' 9.37"E, 48d55'57.87"N)
    Center      (  542464.000, 5447680.000) (  7d 4'45.32"E, 48d56'14.43"N)
    Band 1 Block=256x256 Type=Byte, ColorInterp=Gray
      Overviews: 2048x2048, 1024x1024, 512x512, 256x256
    

    retrieve the image
    command : gdal_translate -of GTiff cad025.xml cad025.tiff
    Input file size is 4096, 4096
    0...10...20...30...40...50...60...70...80...90...100 - done.
    

    get image info
    command : gdalinfo cad025.tiff
    Driver: GTiff/GeoTIFF
    Files: cad025.tiff
    Size is 4096, 4096
    Coordinate System is:
    PROJCS["unnamed",
        GEOGCS["unnamed ellipse",
            DATUM["unknown",
                SPHEROID["unnamed",6378137,298.257222101]],
            PRIMEM["Greenwich",0],
            UNIT["degree",0.0174532925199433]],
        PROJECTION["Equirectangular"],
        PARAMETER["latitude_of_origin",0],
        PARAMETER["central_meridian",0],
        PARAMETER["standard_parallel_1",46.5],
        PARAMETER["false_easting",0],
        PARAMETER["false_northing",0],
        UNIT["metre",1,
            AUTHORITY["EPSG","9001"]]]
    Origin = (541952.000000000000000,5448192.000000000000000)
    Pixel Size = (0.250000000000000,-0.250000000000000)
    Metadata:
      AREA_OR_POINT=Area
    Image Structure Metadata:
      INTERLEAVE=BAND
    Corner Coordinates:
    Upper Left  (  541952.000, 5448192.000) (  7d 4'21.26"E, [COLOR="#FF0000"]48d56'30.99"N[/COLOR])
    Lower Left  (  541952.000, 5447168.000) (  [COLOR="#FF0000"]7d 4'21.26"E[/COLOR], 48d55'57.87"N)
    Upper Right (  542976.000, 5448192.000) (  [COLOR="#FF0000"]7d 5' 9.37"E[/COLOR], 48d56'30.99"N)
    Lower Right (  542976.000, 5447168.000) (  7d 5' 9.37"E, [COLOR="#FF0000"]48d55'57.87"N[/COLOR])
    Center      (  542464.000, 5447680.000) (  7d 4'45.32"E, 48d56'14.43"N)
    Band 1 Block=4096x2 Type=Byte, ColorInterp=Gray
    

    command : open cad025.tiff in Global Mapper and get metadata
    FILENAME=cad025.tiff
    DESCRIPTION=cad025.tiff
    UPPER LEFT X=541952.000
    UPPER LEFT Y=5448192.000
    LOWER RIGHT X=542976.000
    LOWER RIGHT Y=5447168.000
    WEST LONGITUDE=[COLOR="#FF0000"]7° 05' 4.0161" E[/COLOR]
    NORTH LATITUDE=[COLOR="#FF0000"]49° 01' 26.8527" N[/COLOR]
    EAST LONGITUDE=[COLOR="#FF0000"]7° 05' 52.2051" E[/COLOR]
    SOUTH LATITUDE=[COLOR="#FF0000"]49° 00' 53.6816" N[/COLOR]
    PROJ_DESC=Equirectangular / RGF-93 / meters
    PROJ_DATUM=RGF-93
    PROJ_UNITS=meters
    EPSG_CODE=310024802
    COVERED AREA=1048576.2 sq m
    NUM COLUMNS=4096
    NUM ROWS=4096
    NUM_BANDS=1
    PIXEL WIDTH=0.25 meters
    PIXEL HEIGHT=0.25 meters
    BIT_DEPTH=8
    PHOTOMETRIC=Greyscale (Min is Black)
    BIT_DEPTH=8
    SAMPLE_FORMAT=Unsigned Integer
    ROWS_PER_STRIP=2
    COMPRESSION=None
    PIXEL_SCALE=( 0.25, 0.25, 1 )
    TIEPOINTS=( 0.00, 0.00, 0.00 ) --> ( 541952.000, 5448192.000, 0.000 )
    MODEL_TYPE=Projection Coordinate System
    RASTER_TYPE=Pixel is Area
    

    => LONGITUDE and LATITUDE are shifted

    reproject cad025.tiff to RGF93CC49 with gdal
    command : gdalwarp -t_srs EPSG:3949 -ot Byte -rb -of GTiff cad025.tiff cad025cc49.tiff
    Creating output file that is 4224P x 4394L.
    Processing input file cad025.tiff.
    0...10...20...30...40...50...60...70...80...90...100 - done.
    

    get image info
    command : gdalinfo cad025cc49.tiff
    Driver: GTiff/GeoTIFF
    Files: cad025cc49.tiff
    Size is 4224, 4394
    Coordinate System is:
    PROJCS["RGF93 / CC49",
        GEOGCS["RGF93",
            DATUM["Reseau_Geodesique_Francais_1993",
                SPHEROID["GRS 1980",6378137,298.2572221010002,
                    AUTHORITY["EPSG","7019"]],
                AUTHORITY["EPSG","6171"]],
            PRIMEM["Greenwich",0],
            UNIT["degree",0.0174532925199433],
            AUTHORITY["EPSG","4171"]],
        PROJECTION["Lambert_Conformal_Conic_2SP"],
        PARAMETER["standard_parallel_1",48.25],
        PARAMETER["standard_parallel_2",49.75],
        PARAMETER["latitude_of_origin",49],
        PARAMETER["central_meridian",3],
        PARAMETER["false_easting",1700000],
        PARAMETER["false_northing",8200000],
        UNIT["metre",1,
            AUTHORITY["EPSG","9001"]],
        AUTHORITY["EPSG","3949"]]
    Origin = (1998175.253576367200000,8201596.376676536200000)
    Pixel Size = (0.244426308951942,-0.244426308951942)
    Metadata:
      AREA_OR_POINT=Area
    Image Structure Metadata:
      INTERLEAVE=BAND
    Corner Coordinates:
    Upper Left  ( 1998175.254, 8201596.377) (  7d 4'21.40"E, [COLOR="#008000"]48d56'32.69"N[/COLOR])
    Lower Left  ( 1998175.254, 8200522.367) (  [COLOR="#008000"]7d 4'18.57"E[/COLOR], 48d55'57.97"N)
    Upper Right ( 1999207.710, 8201596.377) ( [COLOR="#008000"] 7d 5'12.07"E[/COLOR], 48d56'30.89"N)
    Lower Right ( 1999207.710, 8200522.367) (  7d 5' 9.23"E, [COLOR="#008000"]48d55'56.17"N[/COLOR])
    Center      ( 1998691.482, 8201059.372) (  7d 4'45.32"E, 48d56'14.43"N)
    Band 1 Block=4224x1 Type=Byte, ColorInterp=Gray
    

    command : load cad025cc49.tiff in Global Mapper and get metadata
    FILENAME=cad025cc49.tiff
    DESCRIPTION=cad025cc49.tiff
    UPPER LEFT X=1998175.254
    UPPER LEFT Y=8201596.377
    LOWER RIGHT X=1999207.710
    LOWER RIGHT Y=8200522.367
    WEST LONGITUDE=[COLOR="#008000"]7° 04' 18.5698" E[/COLOR]
    NORTH LATITUDE=[COLOR="#008000"]48° 56' 32.6888" N[/COLOR]
    EAST LONGITUDE=[COLOR="#008000"]7° 05' 12.0680" E[/COLOR]
    SOUTH LATITUDE=[COLOR="#008000"]48° 55' 56.1734" N[/COLOR]
    PROJ_DESC=Lambert-93 / RGF-93 / meters
    PROJ_DATUM=RGF-93
    PROJ_UNITS=meters
    EPSG_CODE=3949
    COVERED AREA=1108868.2 sq m
    NUM COLUMNS=4224
    NUM ROWS=4394
    NUM_BANDS=1
    PIXEL WIDTH=0.2444263 meters
    PIXEL HEIGHT=0.2444263 meters
    BIT_DEPTH=8
    PHOTOMETRIC=Greyscale (Min is Black)
    BIT_DEPTH=8
    SAMPLE_FORMAT=Unsigned Integer
    ROWS_PER_STRIP=1
    COMPRESSION=None
    PIXEL_SCALE=( 0.244426, 0.244426, 1 )
    TIEPOINTS=( 0.00, 0.00, 0.00 ) --> ( 1998175.254, 8201596.377, 0.000 )
    MODEL_TYPE=Projection Coordinate System
    RASTER_TYPE=Pixel is Area
    

    => LONGITUDE and LATITUDE are identical
  • global_mapperglobal_mapper Administrator Posts: 17,238
    edited September 2011
    I think the datum is probably the issue. Are you sure that it should be RGF-93?

    Thanks,

    Mike
    Global Mapper Support
    support@globalmapper.com
  • GeoGeo Global Mapper User Trusted User Posts: 92
    edited September 2011
    hello Mike,

    I checked again, but I'm certainly missing something.

    according to Geoportal's API / Coordinates Reference Systems / Mapping between IGNF and EPSG codes for WMS-C services,
    https://api.ign.fr/geoportail/api/doc/webmaster/rig.html
    
    the WKT definition for IGNF:GEOPORTALFXX / EPSG:310024802 is :
    PROJCS["Géoportail - France métropolitaine",GEOGCS["Réseau géodésique français 1993",DATUM["Réseau géodésique français 1993",SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG024"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","RGF93G"]],PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9001"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6378137.0000],PARAMETER["latitude_of_origin",0.000000000],PARAMETER["central_meridian",0.000000000],PARAMETER["standard_parallel_1",46.500000000],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","GEOPORTALFXX"]]
    


    in Global Mapper / Tools / Configure / Projection
    -> "Init from EPSG" : 310024802
    content of [drive:]\[installation_folder_of_Global Mapper]\epsg_codes.txt is :
    //
    //IGNF:GEOPORTAL
    310024802,Equirectangular,RGF-93,meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=46.5,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    
    then, "Send To File" and save.


    comparison (m=match) :
    PROJCS["Géoportail - France métropolitaine",                                                    |  |  PROJCS["Equidistant_Cylindrical",
        GEOGCS["Réseau géodésique français 1993",                                                   |  |      GEOGCS["GCS_Geographic Coordinate System",
            DATUM["Réseau géodésique français 1993",                                                |  |          DATUM["D_RGF-93",
                SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]], |m |              SPHEROID["GRS_1980",6378137,298.257222101]],
                TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG024"]],                  |  |  
            PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],                              |m |          PRIMEM["Greenwich",0],
            UNIT["degree",0.01745329251994330],                                                     |m~|          UNIT["Degree",0.017453292519943295]],
            AXIS["Longitude",EAST],                                                                 |  |  
            AXIS["Latitude",NORTH],                                                                 |  |  
        AUTHORITY["IGNF","RGF93G"]],                                                                |  |  
        PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9001"]],                                  |m?|      PROJECTION["Equidistant_Cylindrical"],
        PARAMETER["semi_major",6378137.0000],                                                       |  |  
        PARAMETER["semi_minor",6378137.0000],                                                       |  |  
        PARAMETER["latitude_of_origin",0.000000000],                                                |  |  
        PARAMETER["central_meridian",0.000000000],                                                  |m |      PARAMETER["central_meridian",0],
        PARAMETER["standard_parallel_1",46.500000000],                                              |m |      PARAMETER["standard_parallel_1",46.5],
        PARAMETER["false_easting",0.000],                                                           |m |      PARAMETER["false_easting",0],
        PARAMETER["false_northing",0.000],                                                          |m |      PARAMETER["false_northing",0],
        UNIT["metre",1],                                                                            |  |      UNIT["Meter",1]]
        AXIS["Easting",EAST],                                                                       |  |  
        AXIS["Northing",NORTH],                                                                     |  |  
    AUTHORITY["IGNF","GEOPORTALFXX"]]                                                               |  |  
    

    except the naming of the projection,
    Equirectangular (in Global Mapper epsg_codes.txt)
    Equirectangular (in IGNF:GEOPORTALFXX)
    Equidistant_Cylindrical (in Global Mapper prj file)

    everything looks identical.
    I don't understand why there is a shift when I load the same area in IGNF: GEOPORTALFXX and RGFCC49.

    thank you for your attention,
    Richard.

    PS: i've uploaded the 2 files made with gdal (cad025.tiff and cad025cc49.tiff) to your ftp (folder IGNF:GEOPORTAL)
  • global_mapperglobal_mapper Administrator Posts: 17,238
    edited September 2011
    Richard,

    I see one thing in the datum definition that is odd. The ellipsoid defined with the datum is an ellipsoid, but then in the projection parameters the semi-major and semi-minor axes of the ellipsoid are explicitly defined and both are the same, which looks like a sphere. This looks suspiciously like the situation with the Google Maps datum that looks like WGS84 but behaves slightly differently. Try changing the datum name in your EPSG code definition to 'GOOGLE MAPS (SPHERE RADIUS 6378137)' and see if it then lines up.

    Let me know if I can be of further assistance.

    Thanks,

    Mike
    Global Mapper Support
    support@globalmapper.com
  • GeoGeo Global Mapper User Trusted User Posts: 92
    edited September 2011
    Bravo !

    this epsg_codes.txt works.
    //
    //IGNF:GEOPORTALFXX
    310024802,Equirectangular,GOOGLE MAPS (SPHERE RADIUS 6378137),meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=46.5,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    


    i tried to update epsg_codes.txt, the before last, IGNF:TERA50STEREO, isn't complete, i think a custom datum is needed.
    //PROJCS["Géoportail - France métropolitaine",GEOGCS["Réseau géodésique français 1993",DATUM["Réseau géodésique français 1993",SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG024"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","RGF93G"]],PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9001"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6378137.0000],PARAMETER["latitude_of_origin",0.000000000],PARAMETER["central_meridian",0.000000000],PARAMETER["standard_parallel_1",46.500000000],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","GEOPORTALFXX"]]
    //IGNF:GEOPORTALFXX Géoportail - France métropolitaine - Réseau géodésique français 1993
    310024802,Equirectangular,GOOGLE MAPS (SPHERE RADIUS 6378137),meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=46.5,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    //
    //PROJCS["Géoportail - Antilles françaises",GEOGCS["Réseau de référence des Antilles françaises (1988-1991)",DATUM["Réseau de référence des Antilles françaises (1988-1991)",SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG495"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","WGS84RRAFGEO"]],PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9002"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6378137.0000],PARAMETER["latitude_of_origin",0.000000000],PARAMETER["central_meridian",0.000000000],PARAMETER["standard_parallel_1",15.000000000],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","GEOPORTALANF"]]
    //IGNF:GEOPORTALANF - Géoportail Antilles françaises - Réseau de référence des Antilles françaises (1988-1991)
    310915814,Equirectangular,GOOGLE MAPS (SPHERE RADIUS 6378137),meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=15,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    //
    //PROJCS["Géoportail - Guyane",GEOGCS["Réseau géodésique français de Guyane 1995",DATUM["Réseau géodésique français de Guyane 1995",SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG486"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","RGFG95GEO"]],PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9003"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6378137.0000],PARAMETER["latitude_of_origin",0.000000000],PARAMETER["central_meridian",0.000000000],PARAMETER["standard_parallel_1",4.000000000],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","GEOPORTALGUF"]]
    //IGNF:GEOPORTALGUF Géoportail Guyane - Réseau géodésique français de Guyane 1995
    310486805,Equirectangular,GOOGLE MAPS (SPHERE RADIUS 6378137),meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=4,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    //
    //PROJCS["Géoportail - Réunion et dépendances",GEOGCS["Réseau géodésique de la Réunion 1992",DATUM["Réseau géodésique de la Réunion 1992 (RGR92)",SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG700"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","RGR92GEO"]],PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9004"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6378137.0000],PARAMETER["latitude_of_origin",0.000000000],PARAMETER["central_meridian",0.000000000],PARAMETER["standard_parallel_1",-21.000000000],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","GEOPORTALREU"]]
    //IGNF:GEOPORTALREU - Géoportail Réunion et dépendances - Réseau géodésique de la Réunion 1992 (RGR92)
    310700806,Equirectangular,GOOGLE MAPS (SPHERE RADIUS 6378137),meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=-21,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    //
    //PROJCS["Géoportail - Mayotte",GEOGCS["RGM04 (Réseau Géodésique de Mayotte 2004)",DATUM["RGM04 (RESEAU GEODESIQUE DE MAYOTTE 2004)",SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG702"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","RGM04GEO"]],PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9005"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6378137.0000],PARAMETER["latitude_of_origin",0.000000000],PARAMETER["central_meridian",0.000000000],PARAMETER["standard_parallel_1",-12.000000000],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","GEOPORTALMYT"]]
    //IGNF:GEOPORTALMYT - Géoportail Mayotte - RESEAU GEODESIQUE DE MAYOTTE 2004
    310702807,Equirectangular,GOOGLE MAPS (SPHERE RADIUS 6378137),meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=-12,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    //
    //PROJCS["Géoportail - Saint-Pierre et Miquelon",GEOGCS["Saint-Pierre-et-Miquelon (2006)",DATUM["St Pierre et Miquelon 2006",SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG706"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","RGSPM06GEO"]],PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9006"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6378137.0000],PARAMETER["latitude_of_origin",0.000000000],PARAMETER["central_meridian",0.000000000],PARAMETER["standard_parallel_1",47.000000000],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","GEOPORTALSPM"]]
    //IGNF:GEOPORTALSPM - Géoportail Saint-Pierre et Miquelon - St Pierre et Miquelon 2006
    310706808,Equirectangular,GOOGLE MAPS (SPHERE RADIUS 6378137),meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=47,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    //
    //PROJCS["Géoportail - Nouvelle-Calédonie",GEOGCS["Réseau Géodésique de Nouvelle-Calédonie",DATUM["RESEAU GEODESIQUE DE NOUVELLE-CALEDONIE (RGNC 1991)",SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG547"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","RGNCGEO"]],PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9007"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6378137.0000],PARAMETER["latitude_of_origin",0.000000000],PARAMETER["central_meridian",0.000000000],PARAMETER["standard_parallel_1",-22.000000000],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","GEOPORTALNCL"]]
    //IGNF:GEOPORTALNCL - Géoportail Nouvelle-Calédonie - RESEAU GEODESIQUE DE NOUVELLE-CALEDONIE (RGNC 1991)
    310547809,Equirectangular,GOOGLE MAPS (SPHERE RADIUS 6378137),meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=-22,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    //
    //PROJCS["Géoportail - Wallis et Futuna",GEOGCS["World Geodetic System 1984",DATUM["WGS84",SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.000000],AUTHORITY["IGNF","REG642"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","WGS84G"]],PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9008"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6378137.0000],PARAMETER["latitude_of_origin",0.000000000],PARAMETER["central_meridian",0.000000000],PARAMETER["standard_parallel_1",-14.000000000],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","GEOPORTALWLF"]]
    //IGNF:GEOPORTALWLF - Géoportail Wallis et Futuna - WGS84
    310642810,Equirectangular,GOOGLE MAPS (SPHERE RADIUS 6378137),meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=-14,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    

    part.1
    The following errors occurred with your submission
    * The text that you have entered is too long (16384 characters). Please shorten it to 10000 characters long.
  • GeoGeo Global Mapper User Trusted User Posts: 92
    edited September 2011
    part.2
    //
    //PROJCS["Géoportail - Polynésie française",GEOGCS["RGPF (Réseau Géodésique de Polynésie Française)",DATUM["RGPF (RESEAU GEODESIQUE DE POLYNESIE FRANCAISE)",SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG032"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","RGPFGEO"]],PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9009"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6378137.0000],PARAMETER["latitude_of_origin",0.000000000],PARAMETER["central_meridian",0.000000000],PARAMETER["standard_parallel_1",-15.000000000],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","GEOPORTALPYF"]]
    //IGNF:GEOPORTALPYF - Géoportail Polynésie française - RGPF (RESEAU GEODESIQUE DE POLYNESIE FRANCAISE)
    310032811,Equirectangular,GOOGLE MAPS (SPHERE RADIUS 6378137),meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=-15,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    //
    //PROJCS["Géoportail - Kerguelen",GEOGCS["Réseau Géodésique des TAAF (2007) (dms)",DATUM["RESEAU GEODESIQUE DES TERRES AUSTRALES ET ANTARCTIQUES FRANCAISES 2007",SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG036"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","RGTAAF07G"]],PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9010"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6378137.0000],PARAMETER["latitude_of_origin",0.000000000],PARAMETER["central_meridian",0.000000000],PARAMETER["standard_parallel_1",-49.500000000],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","GEOPORTALKER"]]
    //IGNF:GEOPORTALKER - Géoportail Kerguelen - RESEAU GEODESIQUE DES TERRES AUSTRALES ET ANTARCTIQUES FRANCAISES 2007
    310642812,Equirectangular,GOOGLE MAPS (SPHERE RADIUS 6378137),meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=-49.5,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    //
    /PROJCS["Géoportail - Crozet",GEOGCS["Réseau Géodésique des TAAF (2007) (dms)",DATUM["RESEAU GEODESIQUE DES TERRES AUSTRALES ET ANTARCTIQUES FRANCAISES 2007",SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG036"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","RGTAAF07G"]],PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9011"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6378137.0000],PARAMETER["latitude_of_origin",0.000000000],PARAMETER["central_meridian",0.000000000],PARAMETER["standard_parallel_1",-46.000000000],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","GEOPORTALCRZ"]]
    //IGNF:GEOPORTALCRZ - Géoportail Crozet - RESEAU GEODESIQUE DES TERRES AUSTRALES ET ANTARCTIQUES FRANCAISES 2007
    310642801,Equirectangular,GOOGLE MAPS (SPHERE RADIUS 6378137),meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=-46,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    //
    //PROJCS["Géoportail - Amsterdam et Saint-Paul",GEOGCS["Réseau Géodésique des TAAF (2007) (dms)",DATUM["RESEAU GEODESIQUE DES TERRES AUSTRALES ET ANTARCTIQUES FRANCAISES 2007",SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0,0,0,0],AUTHORITY["IGNF","REG036"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","RGTAAF07G"]],PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9012"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6378137.0000],PARAMETER["latitude_of_origin",0.000000000],PARAMETER["central_meridian",0.000000000],PARAMETER["standard_parallel_1",-38.000000000],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","GEOPORTALASP"]]
    //IGNF:GEOPORTALASP - Géoportail Amsterdam et Saint-Paul - RESEAU GEODESIQUE DES TERRES AUSTRALES ET ANTARCTIQUES FRANCAISES 2007
    310642813,Equirectangular,GOOGLE MAPS (SPHERE RADIUS 6378137),meters,CENTRAL MERIDIAN=0,TRUE SCALE LATITUDE=-38,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    //
    //PROJCS["Terre Adélie 1950",GEOGCS["Pointe Géologie - Perroud 1950",DATUM["Pointe géologie-Perroud 1950",SPHEROID["International-Hayford 1909",6378388.0000,297.0000000000000,AUTHORITY["IGNF","ELG001"]],TOWGS84[324.9120,153.2820,172.0260,0,0,0,0],AUTHORITY["IGNF","REG009"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","TERA50G"]],PROJECTION["Oblique_Stereographic",AUTHORITY["IGNF","PRC0901"]],PARAMETER["semi_major",6378388.0000],PARAMETER["semi_minor",6356911.9461],PARAMETER["latitude_of_origin",-90.000000000],PARAMETER["central_meridian",140.000000000],PARAMETER["scale_factor",0.96027295],PARAMETER["false_easting",300000.000],PARAMETER["false_northing",-2299363.482],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","TERA50STEREO"]]
    //IGNF:TERA50STEREO - Terre Adélie 1950 - Pointe géologie-Perroud 1950
    //datum name = XXXXXXXXXXXXXXXX
    //TOWGS84[324.9120,153.2820,172.0260,0,0,0,0]
    //PARAMETER["semi_major",6378388.0000], PARAMETER["semi_minor",6356911.9461]
    2986,Oblique Stereographic,XXXXXXXXXXXXXXXX,meters,CENTRAL MERIDIAN SCALE FACTOR=0.96027295,CENTRAL MERIDIAN=140,ORIGIN LATITUDE=-90,FALSE EASTING (m)=300000,FALSE NORTHING (m)=-2299363.482
    //
    //PROJCS["Géoportail - Monde",GEOGCS["World Geodetic System 1984",DATUM["WGS84",SPHEROID["IAG GRS 1980",6378137.0000,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0.0000,0.0000,0.0000,0.0000,0.0000,0.0000,0.000000],AUTHORITY["IGNF","REG642"]],PRIMEM["Greenwich",0.000000000,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","WGS84G"]],PROJECTION["Miller_Cylindrical",AUTHORITY["IGNF","PRC9901"]],PARAMETER["semi_major",6378137.0000],PARAMETER["semi_minor",6378137.0000],PARAMETER["central_meridian",0.000000000],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","MILLER"]]
    //IGNF:MILLER - Géoportail Monde - WGS84
    310642901,Miller Cylindrical,GOOGLE MAPS (SPHERE RADIUS 6378137),meters,CENTRAL MERIDIAN=0,FALSE EASTING (m)=0,FALSE NORTHING (m)=0
    
    i didn't already made any tests with the epsg custom codes.

    can you tell me how to write the EPSG IGNF:TERA50STEREO for Global Mapper ?

    thanks for your great support,
    Richard.
  • global_mapperglobal_mapper Administrator Posts: 17,238
    edited September 2011
    Richard,

    Create a custom datum based on the International 1909/1924 ellipsoid with the following X, Y, and Z transformation parameters:

    dx = 324.9120
    dy = 153.2820
    dz = 172.0260

    That should give you what you need. Make the name you use for the datum match what is in your epsg_codes.txt file.

    Thanks,

    Mike
    Global Mapper Support
    support@globalmapper.com
  • GeoGeo Global Mapper User Trusted User Posts: 92
    edited October 2011
    Hello,

    With gdal_translate, it is rather complex to retrieve data from the GeoPortail or via an INSPIRE API key to show them in GlobalMapper.

    IGN (Institut Géographique National) released a new version of IGNMap :
    Consult INSPIRE data with IGNMap.
    As part of the implementation of the INSPIRE directive, IGN provides the data flow. To access these feeds, you need to have a key, the latter being available free on the site api.ign.fr .

    see full thread in french
    http://ignmap.ign.fr/spip.php?article9
    
    or in english via translator
    http://translate.google.com/translate?sl=fr&tl=en&u=http://ignmap.ign.fr/spip.php/article9
    


    I tried this new version, exporting layer "CADASTRALPARCELS.PARCELS (EPSG:310024802)" at different resolutions (1m and 0.50m) because all tiles are not available at all resolutions. WMS-C is more performant than WMS for this job.
    Loading all images in GlobalMappper using EPSG code: 310024802 described above as native projection works fine.

    It is currently the fastest way I found to use the Geoportail data with GlobalMappper.
  • global_mapperglobal_mapper Administrator Posts: 17,238
    edited October 2011
    Note I am still working my way through my todo list to get back to WMS-C support, but it is not off my radar.

    Thanks,

    Mike
    Global Mapper Support
    support@globalmapper.com
  • GeoGeo Global Mapper User Trusted User Posts: 92
    edited October 2011
    hello Mike,

    You said that adding WMS-C support is "Not trivial but certainly not rocket science either."
    But it is often necessary to use several tools to achieve our goals.
    So I quoted the Gdal and IGNMap software because they are free and allow to use the INSPIRE Geoportail data with GlobalMappper, not the direct way, but a rather simple workaround.

    Of course, I'd rather do it directly with GlobalMappper and I thank you to make all changes and improvements that users are constantly asking for.
  • papaverpapaver Global Mapper User Trusted User Posts: 25
    edited April 2012
    Hi,
    Mike et Geo : I just want to thank you for this post, VERY usefull for me today.
    Papaver.
  • GeoGeo Global Mapper User Trusted User Posts: 92
    edited April 2012
    some changes are coming soon
    see : Géoportail 2012 changes


    automatic translation (bold and underline by me) :
    Changes to affect how the application accesses the services.

    The main developments in web services in the new version Geoportal are:
    • replacement for the standard flow of tiled WMS-C by the standard WMTS,
    • choice for flows tiled projection EPSG: 3857 and a pyramid OGC compliant (parameters identical to those of Google Maps)
    • the removal of the control system access token (GeoDRM), and the transition to a filtering method compatible with the GIS clients.
    • geocoding on land parcels in addition to place names and addresses.


    Innovations will also be phased in after the opening:
    • the possibility of access in HTTPS
    • dissemination of data vector WFS,
    • flux vector curled (kml) and a 3D API,
    • a service of coordinate transformations,
    • a calculation service elevation.

    From late March, will be issued a test platform services Geoportal three and documentation of new products.
  • global_mapperglobal_mapper Administrator Posts: 17,238
    edited April 2012
    Actually some of those sounds encouraging to get it working with GM right now, particularly the tiled system with parameters identical to Google Maps. Sorry I keep delaying the tiled WMS implementation, seems my todo list is never ending!

    Thanks,

    Mike
    Global Mapper Guru
    gmsupport@bluemarblegeo.com
    http://www.globalmapper.com
  • global_mapperglobal_mapper Administrator Posts: 17,238
    edited May 2012
    Do you know if the switch-over to WMTS (standard tiled WMS) has been completed yet? We now have support for WMTS under the covers in Global Mapper and would like to add support for adding general WMTS sources to Global Mapper.

    Thanks,

    Mike
    Global Mapper Guru
    gmsupport@bluemarblegeo.com
    http://www.globalmapper.com
  • GeoGeo Global Mapper User Trusted User Posts: 92
    edited May 2012
    Un nouveau Géoportail en juillet
    A new Géoportail in july
    L'infrastructure Géoportail entame sa seconde mutation avec une nouvelle version prévue pour juillet 2012.
    The Geoportal infrastructure begins its second mutation with a new release scheduled for July 2012.


    Comment modifier l'application dont vous êtes éditeur pour utiliser les nouveaux services du Géoportail (SIG, application mobile, CMS) ?
    How to change your application to use Géoportail new services (GIS, mobile application, CMS)?
    Depuis début avril, est publiée une plateforme de test des services du Géoportail 3 ainsi que la documentation des nouveautés WMTS, OpenLS, et de gestion de droit d'accès (GeoDRM).
    Il vous est possible de tester en avant-première les web services du Géoportail 2012. Pour cela, il vous suffit de contacter votre correspondant IGN habituel dont les coordonnées figurent sur http://pro.ign.fr/5/contacts-ign.htm
    Since early April, a test platform for Géoportail 3 services has been published and documentation of new WMTS, OpenLS, and management of access rights (GeoDRM).
    You can test a preview web services Géoportail 2012. To do this, simply contact your usual corresponding IGN whose address appears on http://pro.ign.fr/5/contacts-ign.htm

    developer documentation in english :
    WMTS
    Advanced WMTS use

    OpenLS
    Search engines

    GeoDRM
    New upcoming GeoDRM (in middle of the page)



    that's all the news.
    i don't have an api key for this beta test.

    glad to hear that you have support for WMTS !
    i hope that the new GeoDRM will not be an insurmountable difficulty.
  • global_mapperglobal_mapper Administrator Posts: 17,238
    edited May 2012
    I'm finishing up the generic WMTS support now, so it will be in the v14.0 release later this summer. There should be builds available to try this with by the time GeoPortal goes live with the WMTS support.

    Thanks,

    Mike
    Global Mapper Guru
    gmsupport@bluemarblegeo.com
    http://www.globalmapper.com
Sign In or Register to comment.