Dear all,
I'm working on the new projection to handle Meteosat Second Generation images. It is some sort of a "Vertical Perspective" (Snyder) as I understand. A reference C/Fortran code (pixel row-col to lat-lon) is available from Eumetsat website. My draft implementation of that code computes right results, however only under assumption that Earth radius parameter is equal to 1: set _MSG_ = PROJCS["MSG", GEOGCS["WGS 84", DATUM["World Geodetic System 1984", SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree", 0.017453292519943295], AUTHORITY["EPSG","4326"]], PROJECTION["MeteosatSG"], PARAMETER["semi_major", 1.0], PARAMETER["semi_minor", 1.0], PARAMETER["latitude_of_origin", 0.0], PARAMETER["central_meridian", 0.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["m", 1.0]] That looks strange for me. I'd appreciate any hints how to do it in more straightforward (proper?) way. My code is here: https://github.com/mfilocha/geotools/blob/b42be2ffc30784b61bfd114bc273bd534996713b/modules/library/referencing/src/main/java/org/geotools/referencing/operation/projection/MeteosatSG.java Test data: https://github.com/mfilocha/geotools/blob/b42be2ffc30784b61bfd114bc273bd534996713b/modules/library/referencing/src/test/resources/org/geotools/referencing/test-data/scripts/MeteosatSG.txt Regards Maciej ------------------------------------------------------------------------------ _______________________________________________ GeoTools-Devel mailing list GeoTools-Devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/geotools-devel
